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PREFACE 


The importance of Heaviside’s contributions to electrical 
theory is now generally recognized and appreciated. His 
teachings nevertheless are available to only a comparatively 
few; to the many engineers and physicists who could profit much 
by it, the work of Heaviside is more or less a sealed book. 
This may be accounted for largely as due to the novel and original 
mathematical processes he has introduced and applied with 
such extraordinary skill in the solution of many problems; and 
also to some extent to the lack of any attempt to correlate and 
present his teachings in a systematic manner suitable for one 
approaching the subject for the first time. 

Heaviside was not only a great mathematician, he was also a 
great physicist; and it is the knowledge of the physics of the 
problems which guided him correctly in many instances to the 
development of suitable mathematical processes. He concerned 
himself little with formal proofs or rigorous demonstrations. 
As he remarked: ‘‘In working out physical problems, there 
should be, in the first place, no pretense to rigorous formalism. 
The physics will guide the physicist along somehow to useful 
and important results, by the constant union of physical and 
geometrical or analytical ideas.” 

There is really nothing intrinsically difficult about Heaviside’s 
mathematical processes. Once the underlying principles of 
the operational methods are understood, and for that, a 
knowledge of comparatively elementary mathematics only, is 
required, a way is immediately opened to the student for the 
utilization of a considerable part of Heaviside’s work, particu- 
larly in so far as it relates to engineering problems. As an 
illustration, we may refer to the now celebrated Expansion 
Theorem. Electrical engineers could be easily taught its mean- 
ing and use in the solution of problems, and there is no reason 
why it should not have by this time found its way into engineer- 
ing text-books. 

This book is an attempt to present some of Heaviside’s work 
in a manner to appeal to the engineer. It is put forth in the hope 

Vv 
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that it will serve as a suitable introduction acquainting the stu- 
dent with the principles of operational calculus and their applica- 
tions to engineering problems. Great care was taken in the 
choice of material and its arrangement with a view to unfolding 
the subject in a gradual way, and at the same time emphasizing 
throughout, the utility of these mathematical processes in facili- 
tating the solution of engineering problems. All the way through, 
illustrations are introduced at every step in the development of 
the subject to make it more readily comprehensible, and at the 
same time serving to emphasize the utility of the mathematical 
methods. The Expansion Theorem is stressed throughout, 
showing its applicability to the solution of such varied problems 
as Transmission Lines, Electric Filter Circuits, Artificial Lines, 
and others. 

In the choice of illustrations, I have not confined myself to 
Heaviside’s writings, since he concerned himself mostly with 
problems relating to transmission lines and cables. It was 
rather my aim to show the wider applicability of Heaviside’s 
mathematical processes, and for this reason, I sought to include 
varied problems. It is hoped that the discussions of some 
of the problems given here, such as filter circuits, artificial lines, 
and others, have some intrinsic merit of their own, and will be of 
interest to engineers. 

The author wishes to acknowledge his deep obligations to 
Heaviside for the benefits derived from the study of his writings. 
This little book is put forward as a personal tribute to him and 
also with the expectation that it may serve as a stimulus for 
others to study his works. 


THE AUTHOR. 
WasuHInaTon, D. C. 
August 2, 1928. 
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INTRODUCTION 


Doctor Louis Cohen’s studies of Heaviside’s work are well 
known from his publications, and they are highly appreciated. 
Electrical engineering is to be congratulated upon the fact 
that Dr. Cohen in this book gives us a splendid summary of 
that part of Heaviside’s mathematical analysis which bears upon 
the theory of the electrical circuit. It will be found that Dr. 
Cohen’s presentation of this subject is so clear and simple that 
many of the difficulties which the student finds in Heaviside’s 
original work have disappeared. MHeaviside’s work is epoch 
making and every electrical engineer should familiarize himself 
with it, to that extent, at least, which is given in the judiciously 
selected parts discussed in this excellent book. 
M. I. Purr. 


NorFouk, Conn. 
August 20, 1928. 
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HEAVISIDE’S 
ELECTRICAL CIRCUIT 
THEORY 


CHAPTER I 
OPERATIONAL CALCULUS 


The Application of Operational Calculus to the Solution of Circuit Problems 
of Concentrated Inductance, Capacity, and Resistance. 


Any investigation in electric-circuit theory, to be at all compre- 
hensive, must necessarily concern itself with the study of the 
voltage and current distribution in the branches of any circuit 
network during the transient state as well as in the permanent 
state. The problems in electric-circuit theory relating to the 
steady-state condition do not generally offer any mathematical 
difficulties, even for the case of circuits of distributed inductance, 
capacity, and resistance; such problems are now fully discussed in 
engineering textbooks. The mathematical equipment of the 
‘ engineering student is quite sufficient to enable him to follow 
intelligently any discussion relating to such problems. In the 
case, however, of the study of transient phenomena, problems 
arise which are much more difficult, and special mathematical 
methods must be provided to facilitate and, in some cases, to 
make at all possible a solution. 

Heaviside was chiefly interested in telegraphic problems, deal- 
ing with matters relating to the propagation of electric impulses 
in conductors, and as such had to deal exclusively with transient 
effects. This led him to the development of the operational 
calculus, which he used with remarkable skill and brilliance in 
the solution of many problems which would otherwise have been 
either altogether impossible of solution or at least would have 
involved mathematics of great complexity. The mathematical 
methods which Heaviside evolved are, of course, applicable to 
the solution of problems in other branches of engineering, and 
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emphasis should be placed on the methods rather than on the 
specific problems used to illustrate the applicability of the 
methods. 

To begin with, we shall consider a few simple problems utiliz- 
ing direct operational methods to obtain the solutions, which 
may serve as an introduction to familiarize the student with the 
general concept of operational solutions. 

Circuit of Inductance and Resistance.—Consider the case of 
the current rise in an inductive circuit. Suppose we have a 
circuit of inductance LZ and resistance R, and a steady voltage 
E applied, what is the currect in the circuit at any time after 
closing the circuit? We have, of course, the 
well-known differential equation for this cir- 
3 cuit condition, which is as follows: 


I= +h = H#. (1) 


The usual method for solving this equa- 
tion, to obtain an expression for the current, 


Fi@. 1. 
is to assume a solution of the form: 
t= Ae + B. (2) 
which, on substitution, gives 
(LD) + R)Ae‘ + BR = E. (3) 


Since this is to hold for all values of ¢, the following relations 


must be satisfied: 
DBXV+ R = 0; BR = £, 


which give 
R E 
= a B= R 
Hence, 
re 
i=Ae 4+ = 


To satisfy the condition that the current in the circuit is to have 


zero value when ¢ = 0, the constant A must have the value eS 


R 


We finally arrive at the following equation: 


Bove a 
s=Ri-e*) (4) 


the well-known expression for the current rise in an inductive 
circuit. 

We shall now obtain the solution to the same problem by the 
operational method. Use the symbol p to designate differ- 
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entiation to time—that is, p = d/dt—and equation (1) takes 


this form: 
(Lp + R)i = E, 


E 
| "Ip +R . 
This is a symbolic solution, giving the current in terms of the 
constants of the circuit and the symbol p, which is to be alge- 
brized, as Heaviside would say, to obtain the real solution. We 
proceed as follows: Write equation (5) in this form: 


and 


a 


: E E 
eae RY (8) 
P Bais 
by division, we have 
; iy R R? ie 
Be IG a oe 
iy) ig tap te 


The operation on LE by 1/p= p7' signifies integration of # from 
Otot. This is a consequence of the following consideration: 


Given 
pu = 2, 
then 
u= p's, 
and 
pu = ppv = v. 
Hence, 


pp i=1. 


Thus, p~! represents an operation on any quantity that if the 
operation by p be subsequently performed, the quantity is 
unaltered. An operation by p~' is, therefore, equivalent to an 
integration. Operating by p—! on a unity function (that is, one 
which has zero value for t < 0 and is equal to 1 for all values of 


t > 0), we get 
t 
= [aa 
0 


1 
P 
1 t t2 
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Using these integration values in (7), it transforms into the 
following: 
E R WEY feed EE 


“Tal Dee leae Aes 


1 oe) oa Mes (9) 
The bracket expression in (9) is the expanded form of eae 
hence, 7 
Loe 


= To (10) 


a 


R 
Operating now on ek by ; that is, integrating the exponential 


factor, we obtain the complete solution for the current in the 
circuit as follows: 


; EH ’ 2, - | mal 
(= ak to TR € 


E ==, 
which is the same as (4). 


In this simple problem, the operational method does not appear 
to offer any particular advantage over the older method, but it 
serves to illustrate, by a simple example, the general process of 
the operational method. In more complex problems, the advan- 
tage of the operational method will become evident. We note, 
however, that in the development of the solution to the preceding 
problem, we have established the relation 


i ” eats 
a 
eae Te 
or, more generally, 
Z Pp ; 
=——_ =e, (12) 
ieee 
Pp 
Operating on unity function by = r z yields -“. This formula 


turns up frequently, and we shall have occasion to use it in 
connection with other problems. In fact, with the relation 
(12) given, which could have been established without regard to 
any particular problem, the solution to the preceding problem 
would have been much simplified. 
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Circuit of Capacity and Resistance.—Take the case of a 
circuit of capacity C and resistance R, and a steady voltage L 
applied, to obtain the expression for the charging current. The 
circuit equation for this case is 


Ri + ff ea = H, 
symbolically 


Le\ 
(n+ a) = E, 


i i fe (13) 


a 1 1 
ote RI fe Rp) 
This is of the same form as (12), a = 1/RC, and the solution can 
be written down at sight. 


liteby Yo, 


and 


1 
= he Fe (14) 

Circuit of Inductance, Capacity, and Resistance.—As another 
example illustrating the applicability of direct operational 
method, we may consider a circuit comprising inductance, 
capacity, and resistance. The circuit equation expressed 
symbolically is as follows: 


(Ipt+Rt+p\qz 


ye a eee (15) 


i i 
Ui Doan lear aps 


Eig. 3: 


This equation can be put in the following form: 
Ep 


= ) 16 
*= L@— p)@ — pe») oe 
where pi and pz are the roots of the quadratic equation 
LCp? + RCp +1 =0, 
that is, ‘ 


Pp, = —a— j6 
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where 


Beas, 2 AE ee 
a= oF and B= [A 


Equation (16) can be resolved by partial fractions into the 
following: 


eee ee Dot e a 
L(p2— 1) |D— D2 P— Pi 
EH 1 1 
: — el 18 
L(p2 — pr) po Pe 4 oe (18) 
P Pp 


By (12), each of the bracket terms operating on unity function 
is converted into an exponential term; thus: 


steed = €P2t: = eit, 
a hesad 2 
Pp Pp 
Hence, 
: E 
ee ee DOC Se et 
1 Wie ral ePit} (19) 
Substituting the values of p. and p; from (17), we get 
; He-# ; 
= — er Bt _— Bt 
ieee org ae 
_ EKe~* sin Bt 
hes Se 


In the problems considered above, the application of a steady 
voltage was assumed. The operational method, however, is not 
at all restricted to that particular type of voltage; it is equally 
effective in the solution of problems for alternating voltages. 
Since, however, a periodic alternating voltage of any form may 
be resolved into simple harmonic components, it will be sufficient 
to consider only the case of applied simple-harmonic voltage. 
Inductive Circuit, Applied Alternating Voltage.—In the case 
of an inductive circuit, applied voltage simple harmonic, 
E cos wt, the circuit equation is 
(Lp + R)t = E cos ot, 

and 

ee E cos wt. 

hp Re? 
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or 


; D E(ei* a ea) 


2L0( 1 ats i) 


To algebrize the above equation involves operating on an expo- 


(21) 


—1 
nential function by the operator (1 + =) . It is desirable 


to establish first a general formula for this operation. Consider 
the following equation: 


u= d ee (22) 
1+¢ 
er P 
By division, 
(GM RUE ies | 
4 = ie — ef ue sh 23 
(OS OT ie! ee) 


which implies successive integrations of the exponential term 
from 0 tot. Thus: 


t 
pio = | out = Me) 
0 


t 
z 1 Py ; 
pe = fine — 1)dt = OO —1)- < 


t 
1 t| 1 fee 
—3 At — Bas (eh | ee (ee) eee 
PS if | Le ) r| xa ) 2 2! 
Repeating the process over and over, we get for the nth 
integration 
1 t ? fe 
me = eee aoe 1) ath yr-l is, n-29! = X”2-33 ! ee 


Substituting these integration values in (23) and combining 
terms of the same powers of t, we obtain the following: 


w= (1-248 a+. jes 


hee Pe 
(1-8 4% tas ) 
mee a 
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Each of the bracket terms is the expansion of the fraction 


2 hence, 
1+ x 
2 3 
pager aes PS 1-—at (at? _ (at) a) lay te 
eee a nN 2! 3! 
oN 
a 
Rea cea ae eat 
= a é a a5 i. (24) 
Lo x Loe x L-F > 


By the aid of this formula, (21) can be algebrized at sight, 
replacing \ by jw or —jw, andaby R/L. This gives the following 
result: 


eee aE 
EB ejot e let € L € L ; (25) 
tp la ie al eh hee ae 
Lio Lj IR TR 


Operating on the bracket terms of (25) by 1/p, that is, integrat- 
ing between the limits 0 and ¢, we obtain: 


Apa E clot 2 Pal oe ee 1 
~~ 2(R4+Lj0o R+1ljo' R-Ljo R—-Ljo 
R R 
2 el or 1 x gk He 1 
~ R+jle° R+Ljo R—jlw  R—jlw 
R 

nee ¢ jot e—dot - Re | (26) 
~ 21|R+Ljw R — Ljw R?2 + Lo? 


Replacing e' by cos wt + 7 sin wt, and e~/“ by cos wt — 7 sin wi and 
combining, the above simplifies to the following: 


R 
pa | Reswttlwsinwt Re ¥ 
3 R+ D2 - B+ De? 
R 
_ E cos (wt — ¢) ERe © 
VR + ie | RP Da? Ce 
tan ¢ = Ee 
ier 


The first right-hand term of (26) is the steady component 
current, and the second term the transient component, which 
decreases to zero as ¢ increases. 
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The derivation of the above expression was carried through to 
how the straightforward process in operational solutions. The 
same result can be arrived at by a simpler method, utilizing the 
“shift principle’; that is, shifting the operand, the exponential 
factor, to the left, suitably modifying the operator, and 
~ operating on unity function. This will be clear from the follow- 
ing mathematical considerations: 

Shifting of Operand.—If X(t) denotes an algebraical rational 
function of ¢ which can be expanded in ascending or descending 
powers (or both) of the variable, the following relations obtain: 


S(p)ew* = flaje (28) 
For since p stands for d/dt, we have 
pet = aew 
Operating on both sides by p7'!, we get 
pope = pact 
But p-!p = 1, hence, 
aaleat = pleat 
Repeating the operation, we get 
poppet = pact = are, 
also, 
pa p-tet = pla tett = a Pert, 
By repeated operations, we obtain the following relations: 
meat = meat, 
ie ee (29) 
Now as f(p) is an algebraical function which can be expanded 
in powers, we may write, 
Hepet i= (Ay + Arp + Aop? + Asp? +--+ 
B (pee o3 pas bre: eae etree 
(Ao Ase + Asa? + Asa® > + + = 
Bia +- Boa? + Bsa~* + + + + fe™ 
= fae”. (30) 
which establishes the relation (27). 
If T denotes any function whatever of t, then 
Hig) Vertl ja ers(p +a) 7. (31) 
that is, changing p to p +a, shifts the factor e% to the left, 
and the operand is reduced to 7. The proof of this is as follows: 
A single operation with p gives 


petT = ae*T + pT = e*(p + a)T. 


lI 
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Repeat the operation, 


ppe"T = pe*(p + a)T 
ae*(p + a)T + e*p(p + a)T 
e“(p + a)(p + a)T. 


prety = e“(p + aT’. 
By repeated operations, we get the result 


pretyT = ep +a) T. 
Obviously, (32) 
pte eT = & (pp — a)*T.| 


Hence, 


Consider now the case of negative indices; write, 
pPtayT=T1;T =@+4)"?%, 
introduce this value of 7 in the first equation (32), 
pre*(p + a)*T, = #71. 
Operate on each side by p~”, and the result is 
e*(p + a) *Ty = pre T:. 
Now, no limitations were assigned to the form of 7, and there is, 


therefore, none to 7;, which can represent any function of 7; 
replacing it, therefore, by T, we get 


pret T = e%"(p = a)—"T. | 

Similarly, (33 ) 

pe tT = e~#(p ee. a)-"T.| 
Since f(p) is expressible as a series in ascending and descending 
powers of p, we can for each term separately shift the factor 
e~ to the left and change p to p + aineachterm. There results 
the relation 

Spe? = eto aT. (34) 
That is, in operating on a time function which has a factor e*, 
this factor may be shifted to the left by changing every p in the 
operator to p + a. If the entire operand is e*, then shifting it 
to the left and changing p to p + changes the operand to 
unity time function. 

Solution of Inductive Circuit Problem Utilizing the Shift 
Principle.—The utilization of the shift principle established in 
the preceding section is frequently very useful in operational 
solutions. As an illustration of the application of the shift 
principle, we may consider the same problem, that of the cur- 
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rent rise in an inductive circuit under the application of an 
alternating voltage. By (21), we have 
ee (i ciwt ¢— Jot 
a ip Pe Palin oR | 
Consider each right-hand term of (34) separately. For the 
first term, shift the exponential factor to the left, and change 
p to p + jw in accordance with (34); we get 
eiat ee 1 a R+ oe 
= ¢/ = — 1 —— ) 36 
Meer Lemp = Ip\ “Le oe 
and, in this case, the operand is unity. Operating by the bracket 
term on unity, we have, by (12), 


g\7t RF Liw)t 
(Geet re bea, 


and this, again, is to be operated on by 1/p; that is, integrating 
from 0 to t, which gives 


(35) 


t _(R+Ljw)t - iE | _R+Djut 
S50 Feet) eres ae Tow, Say ih. 
ik dt puaiee is 1 
Introducing this in (36) gives 
cjot got if _R+Ljot 
fs Se SE = € L = jl 
Ip +k R + Ljw 
R 
at, clot 37 
= Eayecy Bg ah oR ery (37) 


Carrying through the same operation in connection with the 
second right-hand term of (35), the same result will follow except 
for the change in sign before jw; that is, 
R 
e— Jot ae e— sot 
pie i ie Rie. 
Introducing the values from (37) and (88) in (35) gives the 
expression for the current. 
Se a Geary es PN 
~—6h2NR+A jo R-Ljo R4+ De 
which is the same as (25). 

It is seen that, in this problem, at least, utilizing the shift 
principle has introduced a considerable simplification in the 
solution. It is a useful principle, which can, in some cases, be 
taken advantage of to facilitate operational solutions. 


(38) 


_, 
e ts (39) 
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Circuit of Capacity and Resistance Applied Alternating 
Voltage.—The circuit equation is as follows: 


i(n - | = E cos oat, 
Cp 
sa E cos wt _ ECp(e+t + «~%t) 


1 2(1 + RCp) ou 
Us 
Cp 
Shifting the operands ¢ and e—7' to the left, we get 
._ Bet} Cp + jw) Ee-it| C(p — jw) | 
_ : 2 - (41 
i= “3 [T+ RO@ ca +°5" re Ro@ cae f OD 


Each of the bracket terms is to operate on unity time function. 
Consider each bracket term separately. 


C(p + jw) Hh” Rp 


1+RC(p+Jjo) | 1+ RCje 


operating on unity by feng Ps gives, by (12) a (gotie)! 
: ; : . 
1 + Lt RCje 


RCp 


and this, again, is to be operated on by G — = 


Hence, 
ACOH IO al aye alk - (zetia)! 
L-bRC@ jo)  R a " 
1 t 1 : 
= Je tote) oat ae a3 — (41) 
R+ oe ona 


In a similar way, we get for the second right-hand bracket term 
C(p — jw) _ 1 -(ge-#): lL -(ge-): 
1+ RC(p —jo) BR aay. M4 
Cj 


eels aaa 
Cis 
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Introducing the values from (41) and (42) into (40) gives 


ef wate t Ay es 
i= Size Ro 1 Re + 1 et + Je Bot — 
OSG a sae 
1 = 
ees | e RC —_é* : (43) 
Chew ~ Cje 
This transforms by simple algebra into the following: 
t R La. 
«RG cos wt — a sin wt 
t= Ki— ~ + - — ) 
R11 + R?C%w?) : 1 
| ci | 
z « RC cos (wt + y) | 
yf (at 
1 
tan yp = RCo: 


The first term is the transient component, and the second term 
the permanent-current component. For w = 0, that is, steady 
voltage, (44) reduces to 

ae az 
R 

The derivation of the solutions to the above problems by the 
operational methods illustrates the general process of operational 
solutions. Simple problems were selected with a view to empha- 
sizing the process rather than showing the application to the 
solution of more difficult problems. These will also serve as an 
introduction to the discussion of the expansion theorem, taken 
up in the next chapter, which is a development of the operational 
methods but offers a more direct way of obtaining a solution, 
avoiding many intermediate steps and saving much work. 
Other illustrations of direct operational methods will be given 
in the discussion of circuits of distributed electrical constants. 


CHAPTER II 
EXPANSION THEOREM 


The importance of the expansion theorem in the solution of 
problems in electric-circuit theory cannot be overemphasized. 
It applies to all kinds of circuits; concentrated inductances and 
capacities, or distributed inductances and capacities, giving 
directly the solution for the current and voltage distribution 
in the circuits; the steady as well as the transient state. Of 
course, its greatest utility is in the facility it affords for obtaining 
solutions for transient effects, which may be, in some cases, very 
difficult if not altogether impossible to obtain by the usual 
methods. 

Why this tremendously useful theorem should have remained 
so long unappreciated and not made more extensive use of can 
be explained, perhaps, as due to the fact that Heaviside’s 
mathematical methods, because of their originality and novelty, 
were not fully comprehended for some time and were, therefore, 
neglected. Within recent years, however, more attention has 
been given by physicists and engineers to Heaviside’s work, and 
quite a number of papers were written to elucidate and amplify 
his writings. 

The mathematical formulation of the expansion theorem is as 
follows: 


. EB ePnt 
ae cays err en . 
Pn dp P = Pn 
The meaning of it is this: Given any network of circuits, and a 
steady voltage applied at any point in the system, the current is 
given by the above expression; the first right-hand term repre- 
sents the steady-state component, and the summation term 
gives the transient component. Z(p) is the generalized 
impedance of the system. Z(p)p-o is the value of Z(p) when 
p = 0, and the summation term is to extend to all the roots of 
the equation Z(p) = 0. The meaning and method of using it 
will be clearer as we go through with the process of derivation 
and the application to the solution of various problems. 
14 


EXPANSION THEOREM 15 


The method and reasoning by which Heaviside arrived at this 
theorem is rather difficult to follow, involving, as it does, closely 
theoretic reasoning of a dynamical system, this is given in Chap. 
VIII. It is not at all necessary, however, that we follow his 
method; the derivation may be obtained by a simpler method 
involving only algebraic processes, which the author has given in 
a paper published in the Journal of the Franklin Institute, 
December, 1922. 


DERIVATION OF EXPANSION FORMULA 


Given any network of circuits, the combined impedance of the 
circuit system being designated by Z(p), then, obviously, the 
current is given by 


E 
= Z(p) (2) 


Z(p) always contains p in various powers, and, in general, Z(p) 
can be written, 


LAD ee -lesGa part tid apie teens, (3) 


Let the roots corresponding to Z(p) = 0 be pi, po, ps, ete., then 
(3) may be written in this form: 


Zp a= (P= Ps) (Piz Da) (Pp — Ds) eo sce, 


a 


and 
B | 
Pa) Dy Ps) ar 


The denominator of (4) may be resolved by partial fractions 
into the following: 


(4) 


: 1 Roa! B 
Pip — pap — ps) = =p — pi “Pp — pe f 
@ 
= Pie Rese ’ 5 
pe ee DP — Pm ©) 
where 
the i 
~ (pi — p2)(pi — ps)(pi — pa) + * + (Pi — Dm) 
1 
ON 900s 7): =p) a — Dw) (6) 
E f a | 
~ (pm — Pi)\(Pm — P2)(Pm — ps) * + * (Pm — Dm—) 
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For the explanation of the above transformation, see any algebra. 
The following notation 1s convenient: 
1 1 1 
A=,—;B=,—;50=s5 :: :: 
Z(p1) Z(p2) Z(ps) 
Introducing these in (2), we get 
1 1 


“Gp — p1)Z(p1) as (p — p2)Z (po) Bi 


fi 
1 

(Dp — Din)Z(pm) 

: I ; 1 

Now write each of the factors ———— in the form ————_ = 

es) 


an m Lm 
Pp 


; (7) 


A ae Be) and expanding in a series, we obtain 
P Pp 


il a 114 Pe Bey ey ay 
P—DPm ?p\ 
teens + Pa es | 
1 
= — €Pm! 
p 


Operating on ¢?m‘ by 1/p, which means integration to ¢t from 
0 to t, we get 


P — Pm 
Substituting the relation given by (8) for the factors in each of 
the terms of (7), we get the following: 


t 
pee J, Pntdt = = (nt — 1) (8) 
0 Dm 


: E ePit eP2t ne os €P mt 
anes piZ (pr) DoZ (pa) ae DmZ (Pm) 
LY ee 
PiZ(pi) — p2Z (po) 


(9) 


Now suppose we write, 
Z(p) = (p — pi)F(p) 
F(p) = (p — po)(p — ps)(p — ps) - - - 
By differentiation, we obtain 


IP) _ (p — 
Op 


where 


oF ae 


=P) + PG. 
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For 
P = pi; p — pi = 0; and F(p) = Z(pi). 
Hence, 
OZ . 
ee p= pi= 4p). 
Similarly, (10) 
dZ(p) 


dp Pp = p2= Z(pr2), 


and so on for all other values of p. 
Also, for p = 0 by (5) and (6), 
ete 1 1 1 
Z@)m0 PZ) pis) pls) °°° OY 
Substituting the results from (10) and (11) in (9), we finally 
obtain 


‘ E ( eP it 
t=>——_ +# 
Z(P)pa0 | Own F 
P1 ap foie ay U2 
eP2t ae (12) 
Cae ae 
2 op 2 
We may write the above in a more convenient form, 
é E . sy €Pnt 
(re 13 
Zea + © Dar) oe) 
n=1Pn op n 


which is the form given by equation (1). 

The physical interpretation of the formula is this: Given any 
electrical circuit and a steady electromotive force impressed on 
it, the current in the system is the algebraic sum of two com- 
ponents, one of which is the steady-current component, the 
first term of the formula in which p = 0; that is, the reactance 
terms are zero and only the resistances are effective; and the 
summation term represents the transient current, the number of 
terms occurring in the summation term depending on the number 
of degrees of freedom of the system, the number of roots of the 
equation Z(p) = 0. Also, the transient terms may be in the 
form to show simple decaying currents or decaying oscillatory 
currents, depending on the character of the circuit system, 
Further, the summation term represents, also, the subsidence 
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of the current in a circuit system when the electromotive force 
is suddenly removed. 

It is to be observed that, in the derivation of the expansion- 
theorem formula, it was assumed that all the roots of the equa- 
tion Z(p) = 0 are unequal. For the case where some of the 
roots are equal, the problem is more complex, and we can refer 
here only to an interesting discussion of it by Goto.! Formula 
(13), however, is quite sufficient for most practical cases. 

Extension of Expansion Theorem for Alternating Electromotive 
Forces.—Any alternating e.m.f. may be resolved into.sinusoidal 
components, and it is, therefore, sufficient to consider only the 
case of asinusoidal e.m.f. Assume the e.m.f. impressed on the 
circuit system to be represented by the real part of He, and 
put for brevity \ = jw, equation (2) assumes the form 

Pe Gee 

1 Z(p) € | (14) 
Resolving Z(p) by partial fractions, as in the preceding case, we 
get, as before, 


1 1 
| (= pZp.) * (p—p)Z@) 7 °° * | 
1 
- (15 
. (p — Pm)Z (Dm) | ee 
Operating by 1/p on e“ gives 
i INA : dt = 1 a 
Bs alte dt = y(¢ 18 
a Une BL At — 1 At t 
a ae Idt = 3(e" — 1) - § 
Repeating the process over and over, we get 
1 i il 1 Hert? ee. 
eh a aa os — a 
pr € nn oN hnz-1 t m2 2) ym-3 3! aaa (16) 
Now expanding the term , as previously and operating on 
ery 1 
e“, using (16), we obtain 
1 1 Pi Pe 
eS = NAN i fl Fl 25 
p— pi. >| (at wae yc ) 


Pr Pr, Pi 
il = Shite © 
n Tato ae 
*Goro Macuinort, “Extension of the Heaviside Expansion Theorem,”’ 
Researches of the Electrotechnical Laboratory, Tokyo, Japan, 
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a 
aePa ? 


eee 
eee) 


But 
DIP Dia nD; 1 
Ct ’ 
Xr 2 NG 1 = P1 
r 
hence, 
1 ane Ley ee P1 Be 
ier (1 + pit + Be + 
IN 
pit 
i) 
— 1 1 @ = P1 ePit i 
Pip oN { 
hee 
IN 


Operating now on the bracket term by 1/p, we finally obtain 
oat 1 


Gaik fo arn maa aD 


In a similar way, for the other terms of (15), 
m= fet — or, 


D—DPr — pe 


A = : { et — Pmt} 
DP—-— Pm X— Du 
Substituting these values, equation (15) transforms to 


1 1 
(= p20) 7 — ps2) 


eat ae ts) 


NK pZpy) oo 


— Bet 


ePil eP2t 


(A — p1)Z(pi) as (X — p2)Z (po) q 


ePat 


(\ — ps)Z(ps) qyatnt 


- 2 


Ree gee 
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The first bracket term is ; Zl Aa » which is obvious on com- 
p=xr 
paring with (5) Se (6). In the second bracket term, Z(pn) 


a 


may be replaced by — ) = Pn by equation (10). Hence, (18) 


reduces to the eck 


; Ee — ePat | 
fee : 
Tie aZ 19 
ae Or Ciel 29) wee) P = Pn ue 
Pp 
Replacing \ by jw, we have, finally, 
Heit ec EPrt 
. =H) é ~ 
Dering aZ(p) (20) 


n=1 (jw eae p= Dn 


This we may consider the modified expansion formula applicable 
for alternating voltage, and which we shall refer to hereafter as 
the second expansion formula. The first right-hand term of (20) 
is the permanent-current component, and the summation term 
the transient component. It is clear, of course, that only the 
real parts are to be used in the solution of any problem by the 
application of this formula. In the summation term, however, 
the imaginary part is readily eliminated. We may write (20) 
in this form: 


Heit ie 
ER : 
eee aN ss aZ(p) (21) 
n=1 Op P= Pn 


Thus, the steady-current component is expressed in complex 
quantities of which the real part only is to be used in the final 
solution, and the summation term giving the real solution for the 
transient components. 

If the applied voltage is at some phase angle 6, that is E cos 
(wt + 6), then in place of Ee, we have Hef@it+® = Hei%iot, and 
in the final formula, E is to be multiplied by the factor e®, which 
introduces a modification in (20) as follows: 


; Eeiet+o ’ sanity ePnt 
t= 57. — Ee? oh 
P)p=i AG ame ee (22) 


op 
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Taking only the real part of the summation term gives 
n=m 
2 Beer) E S COS (gn + O)e?n! 


— cae : 
Z(D) p= 0 a pe a TAO) ey (23) 
a zs yp ; 


tan Qn = me 
Pr 
We shall now consider a number of problems to illustrate the 
application of the expansion formula to the solution of electric- 
circuit problems. We shall confine the discussion in this chapter 
to problems relating to circuits of concentrated inductance and 
capacity. Other illustrations of the application of the expansion 
formulae to the solution of more difficult problems, those relating 
to circuits of distributed electrical constants, are given in the 
chapters following. 
Circuit of Inductance L and Resistance R.—Assume a steady 
voltage applied to an inductive circuit, the current equation is 


di ; 
Li + Ri = E, 
and 
E 
“Ip tk ee 


The determinantal equation for this case is 
Z(p) =Ilp +k =0 


a one-degree equation, only one root, p = —R/L; 


ne = J; and, for p = 0, Z(p)p=0 = BR. 


Substituting these values in (13) gives the solution of the problem, 


eS BE a ht 
-f- ee 
L 
= (1 7 ett), (25) 


the well-known formula for the current rise in an inductive 
circuit, obtained in the preceding chapter by direct operational 
process. 
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For alternating e.m.f., apply the second expansion formula 
(21), which gives 


Tee ee 
ree Hee 10 cabo 
Tjw + R @ \L 
TL? 
Heit ER ue 
TUiseeR: Le ae ee) 
Taking only the real part of the first right-hand term, 
ee eecas t+ Lw sin Peay Po 
,= R? + List Ww [63] Ww ) 
EH ae t| 
St COS (wt — ¢) — COS ge L = (27) 
a/R? Do? | 
tan ¢ = =e 


The first term in the above equation is the permanent alternating- 
current component, and the second term the transient compo- 
nent. For ¢ = 0,7 = 0 as we should expect, and fort = « 


1 = ———— 0s (wt — (2) 
Wir 
only the permanent-current component active. 
Circuit of Capacity C and Resistance R.—The voltage across a 
condenser for a variable current in the circuit is V = 1/Cfidt, 
and the circuit equation, therefore, is 


(orgies 


ie (28) 


and 


The determinantal equation is 
1 
Zp) = R + G = 0 


which gives only one value of p. 


eee 
P RC 
aZ(p) _ _ 1 _ _me 
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For steady voltage, 


1 
Z(p)p=0 = © because ae o for p = 0. 


The steady-current component is zero, in this case, as we should 
expect. Substituting the values of p and dZ(p)/dp in the sum- 
mation term of (13), we obtain the following expression for the 
current in the circuit: 


t 

: He Fe Beat 

ay ee OF ey Ce 
RC 


the well-known expression for the charging current of a condenser 
circuit. 

For an applied alternating voltage, He’«‘, use the same values 
of p and 0Z(p)/dp in formula (21). 


pe =o 
Eeiot RC‘ 
as 1 1 
xe poeceeas 4 2 
aerarr bee sae \r : 


Taking only the real part of the first right-hand term, and 
rearranging slightly the second term, gives the following: 


1 
. B(R cos wt — G sin ot) He FG! 
i= : f+ = nae 
2 + Ca RC*a*( Fe + Tis) 
which may be written in this form: 
agin E sin y eee , 
1= Seely {eos (wt + yw) + RCw* RC (30) 
DN chy 
tan y = we 
RCw 


the complete solution for the current in a condenser circuit; the 
first term the permanent current component, and the second term 
the transient component. 

Circuit of Inductance L, Capacity C, and Resistance R.— 
For a circuit of inductance, capacity, and resistance in series, 
applied steady voltage, the circuit equation is 


di ee be 
La + ht + G it = E, 
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and 
t= Bene eT : (31) 


( 
DB iol Cp 
The determinantal equation is 
i) 
Ap) =oL R+— =0, 
(p) js ad de Cp 


or 


Ley? + RCp + 1 = 0, 


a second-degree equation, the roots of which are: 


Rigs tee sf] 
Pao Oe INI Ce Kw tee 
CLs eee 32 
R 5 OL Re? oo 
BE OT, wit NO ate 
We also have 
OZ(p) _ 5 _ ae. 
ergo L Cp (33) 
Introducing the values of p; and ps2 from (32) into (33), we get 
ACU ie tei NL 
OLD eee, 
ee ag Cea S35" 


For p = 0; 1/Cp = ~ and Z(p) = o, the steady-current com- 
ponent is zero, as is to be expected. 

For the transient current we have, on substituting the values 
of p and dZ(p)/dp from (33) and (34) into (13), the expansion 
formula 

eit 


Gat iL — oa ey) 


e— JBt \. 


i 
Cla +58) 
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By a simple algebraic transformation, the above reduces to 


RB — ot, | ye eae hs 
SOT PIAA 
Ee—“ sin Bt 
a damped oscillatory current, the damping factor is 
— R , 
aoe 
and the frequency is given by 
- 5 Nee 4L? (36) 


If we neglect the resistance term in the expression for the fre- 
quency, (35) simplifies to 


—ataq) Tal 
_ He“ sin Bt _ B, [Ce-etsin Bt (37) 


1 
—— 
1/LC 


The above results hold also for the case h?/4L? > 1/LC. When, 
however, 1/LC = R?/4L?, the critical case, the two roots p; and 
p2 are equal and the expansion formula in the form given is no 
longer applicable. We may arrive, however, at the expression 
for the current for this case in the following way: When p; and 
2 approach equality, 8 approaches zero value. For very small 
difference between p; and pe, B is very small, sin Bt = Bt, and (35) 
reduces to 


C= (38) 


independent of 8, and holds also for the limiting value 
Pi = p23 B= 0. 


The expression for 7 for the critical case given by (38) can be 
arrived at in another way. Put the circuit equation (31) in 
this form: 

ae 


1 a ) 


1c ( 7 ip 25 ro) 
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and for 
1 | ee 
Le ca 4]? a’, 
: Ep 
‘=To +a) ee 
We note, however, that 
aN ee are 
(pta)? da&pta/)’ 
and 
- f “ie l aS; e—* (see equation (12), Chap. I). 
1+° 
Pp 
Hence, 
d 
©: oe a? = mr es emo, (40) 


Introducing this value in (39), we obtain 


age 
e ) 
which is the same as (38). 

For the case of an alternating e.m.f. impressed on the circuit, 
the modified expansion formula (20) or (21) is to be applied. 
It is somewhat simpler to use formula (20) in this case. 
aZ(p) 

Op 
given by (82) and (34). The expressions for 0Z(p)/dp may be put 
in a simpler form: 


The values of pi, po, and 


are, of course, the same, 


SAN Slaps taxon, | ae 1 _ LCC? — 8 — 2jaB) — 1 
Opes Gla E98)" Ce jey 
I 1 : 
auc ars mete jab) 
~~ C(=a + 58)? 
ai Ah & 2bB 
=e yo. | Parge 
Similarly, 
OZ (p) 2Lj8 _ 2L8 


Op eer a typ Bb ja 
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Introducing these values in (20), we obtain 
Eet ne es” eer (8 1 Ja) 


t= ; a+ 
jij far aig SO re ea 
Chw 
(8 — ja)e~i 
jotatge, “OY 


The first right-hand term in the above expression represents the 
steady-state, periodic component; and the other two terms, the 
transient component. Combining the two terms of the transient 
component, we get 


He 
2LB 


[98a + j6? — wa + jan ee ot (7B SR OO a Je eae 
6? + a? — w? + Qjee 


Di 


This reduces to 

He~~ {jBw cos Bt — (8? + a” + jwa) sin Bt} 
Lg 6B? + a® — w? + Qjwa 

The real aa of (42) is 

ee caper cos Bt — [ (ate 6")? lee — B)] sin Bt} 


(42) 


= 


ie 


1B (B? + a? — w?)? + 4atw 
which may be put in the simpler form 
os lee eo ce em (By) : (43) 
ole LB V/ (8? ot — wt)? + Oe: 
en 2aBbw 


(a? + 6)? + w(a? — 6) 

Taking the real part of the first right-hand term of (41) and com- 
bining with (43), we have the complete expression for the 
current: 


E cos (wt — ¢) Ee-* (a? + 6) sin (Bt — y) | (44) 
Rr EP ee ort aae 
(tym 


For w = 0, the steady-current component is zero, and the 

transient component reduces to 
Ee-# 
Lg 


= 


sin £t. 


= 


the same as (85). 
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We have so far considered the application of the expansion 
theorem for the determination of the current in a circuit. The 
method, however, is not at all limited to that particular variable. 
We could have, for instance, in the preceding illustration, 
developed the solution for the voltage across the condenser, by 
the application of the expansion theorem to the expression for the 
voltage across the condenser, 


E 
~ LCp? + RCp +1 


V 


In this simple case, there is no particular advantage in develop- 
ing the solution for the voltage in preference to the current. 
In the case, however, of problems relating to circuits of 
distributed electrical constants, it is sometimes an advantage to 
apply the expansion theorem to determine the voltage instead 
of the current at some point in the circuit 
and derive the solution for the current from 
the voltage by the circuital relation between 
them. 

Divided Circuits.—A circuit comprising 
an inductance and a capacity in parallel, 
and this parallel circuit in series with a 
resistance, the circuit arrangement is shown Fie. 4. 
in Fig. 4. 

If we designate the impedance of the inductance branch by 
z, and the capacity branch by 22, we have the circuit equations 


Roto a 2111 = E, | 


Roto + Zale = E,, (45) 


and the auxiliary relation 
10 = 14 a 19. 


From these, the expressions for 7; and 72 are readily derived, thus: 


11 == Eas 
Ro(2i + 22) Fez. 
te ae | a 


~ Ro(e: +. 22) a= 2122 
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In one case, the determinantal equation is 


Ro(21 -- 22) + 212 


Z(p) = - =0, 
and in the other case, (47) 
ees ee IEG) 


The roots of the determinantal equation are the same in both 
cases, determined from the equation 


Ro(Z1 + 22)+ z122 = 0, 
1 1 
Fa Lip + Rit Ga) =i Cup (Lip - Rk) = 0, 


which yields the following values of p: 


PTO 8 
where 
See ae 
Sa a2Ty Oh Cs (48) 
_— (Roti (R: aye 
os \ ReEAG: (sr. a siae.) 
For 
ee ol 
() As 2 Cop ae 
Hence, 


for One ACD) 56 — Ry + Ry 
for 72; Z(p)po0 = ©. 


We have also for 7, 


dZ(p) Pre ( 1 ) ele ale ae 
Op = ap | Ro Lp Ry Cop “EAC Ri) Cop 
R Rk 
= Cw | Roba - ft, - P|: (49) 


remembering that the bracket term is zero for the values of p 
which are the roots of the determinant equation, and these are 
the values of p to be used in the expansion formula. For 72, 


| 


1 1 
sain) _ a |RALD+ Rt ap) toghe + FO) 
dp ~—sSOp Iyp + Ri 
= RoLiC op? =a ort Ri) (50) 


Cop*(Lip + Rx) 


30 HEAVISIDE’S ELECTRICAL CIRCUIT THEORY 


Introducing these values in the expansion formula, we obtain the 
following expressions for the currents in the two branches: 


E Ee(-@+ia)t 
io as Ro+ Rk, ae RolaCop? -- (Ro + Ri)’ (51) 
is EC.p(Lip + Ry) e(— ot JB)E, 


~ Reln Cap — Che eR) 


Substituting the value of p given by (48), the complete solutions 
for the current distributions in the circuits are obtained. 
The current in the main circuit is 


a es zit + 1,Cop? + RiCop} 
Role” hdc gee re 


fo=titie= e(—a+i8)t (52) 
When C = 0, that is, the condenser branch open, the expression 


for the current 7) should reduce to that of the current in an induc- 
tive circuit whose resistance is the sum of the resistances Ry and 


R,;. This it actually does, as is seen from the following considera- 
tion: For C very small, 


ee eL, 
MOR Cn Mae, 


= 1 1 R,\ 
a ee > (oom: a5 a 


= eras 2 ae fy (ne lectin oS 
SANTO. ACR one, ae 


oe 1- 2C ok -- RiC2Ro\, approximatel 
= IOO2Ro 6 je | ppro? Wy 
and 
ae Begs eee 1. _ Bo ls 
PSST de ER ee ies 2k ere, 
ie ieee 
+. Ty ee 
Putting C = 0 in equation (52), it reduces to 
RotR1 
do di ee a (54) 


ae eR ye Saree ee 


the expression for the current in an inductive circuit. 


EXPANSION THEOREM ol 


For Li = ©, R; = ~, the inductive branch open, 
a) 
21 ~ 2C2Ro 


ee | aa 
TCs LORY a ee) Ben ie 


= i it - #) USter 
IN GYGRI FE OY Pa fn Te am 


a 


Nine A oe eee Ry 
Mabe = ole SCR, CR. | OF 
1 d 
Pens Om 
For 
E 
fie tare Sree 
and 


Hel a IC 2p? ++ R,Cop) irs E(1iC op? ao R, Cop) 


“RolCop? — (Ry ~ Ri) RoLiCsp? — Ry 
Substituting the value of p from (55), it reduces to 


Introducing these values in equation (52), we have, finally, for 
the condition L; = Ri = ~, 

; Yh ee 

oo = Re RoC2, (56) 
the expression for the charging current in a circuit of resistance 
and capacity. 

Coupled Inductive Circuits.—Two circuits each consisting of an 

inductance and a resistance, coupled magnetically, a steady 
voltage applied to one circuit. The circuit equations are 


Kia. 5, 
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(Lip + Ry)ts ak M pie = E, (57) 
(Lop + Ro)i2 + Mpir = 0, 
from which the expressions for 7; and 72 are readily obtained. 
E(Lep sia Re) 
(LiL, — M?)p? aii (Riles be RoLi)p + ft Re 
= —EMp ; 
*2 = (ile — M?)p? + (Rile + Rola)p + Bike | 
To develop the solutions for the currents in the circuits from the 
above expressions, it is necessary, of course, to obtain the roots of 
the determinantal equation, 


Z(p) = (LiL = M?)p? == (RiLe a RLi)p == RR, = 


The roots are readily determined as follows: 


11 — 
(58) 


pi| _ — (Rule + Rls) + V(Rila = Rel? + 4B iRMs (59) 
Do 2(L,L2, — M?) 
This may be put in a more convenient form, 
Pi a —(ay a a2) + Va; — Qe)? + dora.” (60) 
Pe 1 — Kk? _ 
where 
x Pee _ R,. Po pee be 
ROT a aS, Tai 
We also have 
OZ 
a = 2p(Libe — M?) + (Rule + Rela) 
=> 20,L2{ (1 = K?)p a a) _ ay} 
= +20 iLeV/ (ar an a2)? + 4k?eajag. (61) 
For 
p = 0: 
Z(p)p—0 = R, for the primary circuit, 
and 


Z(p)p—0 = 0 for the secondary circuit. 
which is obvious from equations (58). 
Substituting these values in the expansion formula (13), we 


obtain the following expressions for the currents in the two 
circuits: 


11 = E oe E {Lopi-t Re . p 
R, DT Lente — a2)? ae TERT ae Pi 
(Lope + Re) ne 
SS SS = ae ee 

Pe 


[erasers (62) 


—EM 
A Gel bg eVe Coe —— a2)? + 4aja,k? 


i 
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The values of pi and pe are given by (60). 
For t = 0 the bracket term in the expression for 7; reduces to the 
following: 


Tep1 of Re 7 Lepe as Rk, = Re(pe2 = Pi) 
P1 P2 Pipe 
Introducing the values of p; and pe from (60) gives us 
Lep1 + Re a Lope + Re = 


P1 Pe 

—2ReV (a1 02) oF 4ajack? ; 
_ x2 | (+ a2)? = (a1 = a2)? — 4eraek? 
a | eee 

— 2ReV(a1 — a2)? + doraeK? 

4aja2 

Substituting this in the first equation (62), we have, for t = 0 
oe QRoE Hi oak 


Pie R, = Ranenlile ae R, x R, 
The current in the secondary circuit is obviously zero for ¢ = 0, 
the bracket term reduces to zero. 

Coupled Oscillatory Circuits.—For two circuits each consisting 
of an inductance, capacity, and resistance, the circuit equations are 


1 
(Ln + Ri + ae + Mopre = 


1 
(Le seb) eae on) + Mpii 


| 
& 


(63) 


I 
© 


Fig. 6. 


Solving for 7; and 72, 
1 
allan + Re ¢5) | 
1 
(Ln + hi+ opi? selva oa) me huikg ys (64) 
—EMp 
ne 


= ah — 2 " 
(Ln +R,+ Cap) +e + Ga) Mp 


ja = 


12 = 
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It is observed that in the above expressions for either 71 or 72, 
the determinantal equation is a fourth-degree equation. To 
develop the complete solutions from the above expressions by 
the application of the expansion theorem requires the determina- 
tion of the four roots of a fourth-degree equation, which would 
introduce considerable mathematical complexity. This problem 
is considered fully in Chap. III, Filter Circuits, where questions 
relating to multiperiodic circuit systems are considered. We 
shall limit the discussion here to circuits of negligible resistance, 
in which case the determinantal equation is reduced to a biquad- 
ratic equation, the roots of which are readily determined. 

For R = 0, the circuit equations (64) reduce to the following: 


(Lap ae Cp Ore 

OS C:Csiilig = Mp Cale Ep 
—EMC,C2p3 

CiCo(LiLe — M*)p* + (L2C2 + 11C3)p? + i 


12 = 


For either circuit, the detexminantal equation is 
Z(p) = (Le — M*)CiCopt + (L2C2 + LiCi)p? + 1 = 0. (66) 
The roots of this equation are given by 


0, .= + {th EV CLG? FAO (67) 
2(Lile = M?*)C1C2 


The four values of p, given by (67) may be put in this form: 
Pi = jB1; pe = —JB3; 


ps = jb2; Pps = —JB2, (68) 
where 
ae a, + L2C2) +-V EiC1 — L2C2)? + 4,0, 
2(LiL, — M)C,C, 
Bo =< ome + D2C2) ak has — L2C2)? + 4M?C,C. fy 
Tale — MCCS 
For L,C, = LC, the above reduces to 
B i! 
1 SSE 
VLC + K 
Sia (70) 
B 


NIN LOC Tan 
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where 
M 
K SS 
VLiLs 
The values of 0Z(p)/dp for the different values of p are readily 
obtained from (66). 


hee = 2p{2(LiL, — M?)p® + (LC. + LiCy)}) 


and substituting for p? the value given by (67), it simplifies to 


Se = +2pV/(LiC1 — L2C2)? + 4M C2. (71) 


For p= 0, Z(p) = © obviously. The steady-current com- 
ponent in either the primary or the secondary circuit is zero, 
which we should expect. By substituting the values of p and 
0Z(p)/dp in the expansion formula (13), we obtain the complete 
solutions for the currents in either the primary or secondary 
circuit. Consider the secondary circuit: 


pier 
= —EMC.C Se ee 
; De +2p? /(I1C, — LaC2)? + 4M°C,C, 


~EMC1C 
RC 10. rr eR aa tue 
4 ~EMC,C, 

MON). — L2C2)? + 4M?C C2 


{7B ie! — jBi1¢—763# = 
JB2e82! + 7 Boe} 
{B, sin Bit — By sin Bot}: (72) 


EMC,C2 
eon) (0, = 1,0.) + 40°C, C, 
The current consists of two components, oscillatory currents of 


different frequencies; f; = 61/27 and fe = B2/27. When the 
circuits are syntonized, that is, LiC; = L2C2, the above reduces to 


rs 5 VEC: sin Bit — B2 sin Bet}; (73) 


and the values of 6; and £2 are given by (70). 

Subsidence of Current in Circuits——The problems we have 
considered thus far all relate to this question: A voltage, steady 
or alternating, is suddenly applied to a circuit system; what is 
the character of the current at any time after the application of 
the voltage? We found that, in every case, the formula for the 
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current, which is, indeed, the expansion formula, consists of two 
parts, one of which is the permanent-current component, and 
the other the transient-current component. In some cases, the 
permanent current component is zero, as in the case of applied 
steady voltage to a condenser circuit. The transient current is 
expressed by the summation term in the expansion formula, the 
number of terms occurring in the summation depending on the 
character of the circuit system. 

The expansion theorem, however, is not limited in its applica- 
tion to the determination of the transient-current rise in a 
circuit after the closing of the switch or the application of the 
voltage; it also gives correctly, except for a reversal of sign, the 
subsidence of the current in a circuit system when the voltage 
source is suddenly short-circuited. This can be demonstrated 
in this way: Suppose the circuit was closed for a sufficiently 
long time to insure the establishment of the steady-state condi- 
tion, the transient component completely wiped out, which, in 
most practical cases, obtains within a short time after the 
closing of the switch; the current in the circuit is then expressed 
by the first term of the expansion formula, thus: 


Piaee Bute 
Z(P)p=0 


eee 
Z(D)p =e 


depending on whether the voltage is steady or alternating. Now 
suppose that, instead of short-circuiting the voltage source, we 
introduce in the circuit another voltage of exactly the same 
character and magnitude but of opposite sign; the effect is the 
same as if the original voltage were removed or short-circuited. 
We have then, in this case, two similar voltages of opposite sign 
acting simultaneously on the circuit. The original voltage pro- 
ducing a current in the circuit 


or 


ata ee 
Z(P) p=0 


and the suddenly applied additional voltage producing a current 
in the circuit in accordance with the expansion formula 


; E ePt 
te S7Gieaan >» aZ(p) (74) 
Pap 


D 
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The resultant current in the circuit is the sum of the two 
currents, that is, 
e?! 
aZ(p). (75) 
Op 
which is the formula for the subsidence of the current in a circuit 
on short-circuiting the voltage. The summation term, therefore, 
of the expansion formula gives also the current subsidence in a 
circuit on short-circuiting the voltage. 


CHAPTER III 
ELECTRIC FILTER CIRCUITS 


The term electric filter is used to designate a circuit system 
consisting of a number of recurrent sections connected in tandem, 
each section comprising inductance and capacity in various com- 
binations depending upon the requirements for which the system 
is designed. Actually, a circuit system of this kind is nothing 
more, of course, than a circuit combination comprising several 
circuits coupled electrically and is characterized by a multiplicity 
of degrees of freedom depending on the number of sections or 
individual circuits in the system. That is, the number of free 
vibrations of the circuit system is determined by the number of 
individual circuits. Obviously, the system will respond reson- 
antly to the frequencies corresponding to the frequencies of the 
free vibrations of the system. It permits to pass freely several 
frequencies within a predetermined range governed by the 


number of sections and the choice of the electrical constants of 
each section. <A better term to designate a circuit system of this 
kind would have been multiple periodic circuit system. Since, 
however, it is now widely known under the designation filter 
circuit, the term is retained here. 

Consider a circuit system diagrammatically shown in Fig, 
7. The series impedances are designated by 21, and the shunt 
impedances by 22. 

Either 2: or 2. may be any combination of inductance and 
capacity. An e.m.f. is applied to the first section through an 
impedance z;/2, and the last section is closed through a series 
impedance 2;/2. The mathematical analysis is much simplified 
by this choice of terminations. We shall also designate the 
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currents in the series elements by subscripts corresponding to 
the number of the section. For any section, say the qth, except 
the first and last sections, the voltage drop in the circuits in 
accordance with Kirchoff’s law is given by the following equation: 
Zitg + 22(tq — tq41) — 22(tg-1 — tq) = 0, 
or (1) 
(21 aia 222)tq ar 22(tg44 + Way) = 0. 

We assume a solution of the form 

ty = Aet¥ + Be-u, (2) 
where A and B are independent of g. On substitution and after 
a simple rearrangement, we get 


(2+ _ — ey — e~7)(Aety + Be-27) = 0. (3) 
2 


Clearly, (2) is a solution of (1), provided the value of y is chosen 
to satisfy the relation 


ey +e—t= 2+ ae 
&2 
or 
M Lae 
coshy = 1+ op (4) 


The constants A and B are determined from the termina. con- 
ditions, that is, the circuit equations of the first and last sections. 
For the first, that is, the zero section, we have 


gio + ealto — th) = E, (5) 
and for the last, the nth section, 
- rad 23nd Sy 1) = 0. (6) 


Using the values of 7 given by (2), we obtain the following two 
equations from which the constants A and B are determined: 


e + 22 —- ze") + G + 2 — zac”) B = H, 
21 a &1 — pal 
@ + 2. — 2€ *)Am + é + 2 —-— 220”) OY ts (0): 


Replacing . + 2. by 5 (er + ¢—7), the relation given by (4), and 


(7) 


rearranging, we obtain the following: 


—2H 
meets ae | (8) 
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From these, we obtain 


: 2he 
A = Zo(ev ae e~7)(e" er) | ¥ 
B QHev (9) 
Ga e~7) ("7 — e—”7) | 


Substituting these values in (2), we finally get the complete 
expression for the current in any section, as follows: 


QH{e"—-D7 + Ee Ort 


a Z2(e% — e—V) (7 — e-"7) ce) 
or 
. _ Eecosh (n — g)y 
“¢ * 3) sinh y sinh ny (it) 
For the last, the nth section, g = n, 
in = Z (12) 


2. sinh y sinh ny 

The complete solution can be readily developed from (11) and 
(12) by the application of the expansion theorem. In practical 
work, the main concern is the output current, the current in the 
last section given by (12). The solution and result will depend 
on the particular combination of inductance and capacity in the 
series and shunt elements z; and 22. We shall consider several 
typical arrangements for steady and alternating voltages. 

Series Element z; = Lp + R and Shunt Element z. = 1/Cp 
Steady Voltage E.—The determinantal equation is 


Z(p) = 22 sinh y sinh ny = 0, (13) 


Fie. 8. 
which gives 
ny Sse se= OF 1-2, 3, 2 a, 
ST 
Be (14) 


The zero value of s is to be disregarded, for the reason that this 
would not satisfy the relation given by (4), leading to the result 
21/2, = 0, an impossible condition. Also, all values of s greater 
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than n will give only repetitions of preceding values. Hence, 
values of s from 1 to n are to be used only. The values of p 
corresponding to the roots of the determinantal equation (14) 
can be determined from the relation given by (4). In some 
cases, 22, which is a factor in the determinantal equation, may be 
a combination of inductance and capacity to yield additional 
roots which should be taken into consideration in the complete 
solution. We shall have occasion to discuss this farther on. 

For the case under consideration, 2: = Lp + Rand 2, = 1/Cp, 
we have, by (4), 


cosh y = cos - =} + 5Cp(Lp + R), (15) 
which, on solving for p, gives 
(03. = 3 sie alee (16) 
st hie 
; Be 17 
a =258.= aoe cos ) are (17) 
also 
OLD) mo a: d(ny). 
aes Z2 sinh y cosh ny ap (18) 
By (15), we have 
sinh 2% = LCp + 
dp 
and 
LCp + ne 
oes me (19) 
Op sinh y 
Introducing the value of 0y/dp from (19) into (18) gives 
OZ(p) _ RC 
op na( 2p + 3) cos (sr) 
= ni -+ >) cos (sr), (20) 


for p = 0; Z(p) as given by (13) is indeterminate for sinh y 
and sinh ny are separately equal to zero, and 2 is infinity. We 
can arrive, however, at the value of Z(p)p-0 = 0 by taking its 
value as p approaches zero. 

For small values of p, 


2e 
cosh y=1+% = 1+ Pap +R), 
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and 
vy’ = Cpo(Lp + R). 
; _ ne? _ nCp(Lp +R) _ = 
Z2 sinh y sinh ny = Cr Cs = nk for-p.= 0. (21) 


Introducing the values of p,, OZ(p)/dp, and Z(p)p—o from (16), 
(20), and (21) in the expansion formula (13), Chap. II, we 
arrive at the solution for the current in the nth section as follows: 


=. (—a+jBs)t 
Pa eS i80(1 + ae jp, nt cos (sr) 
E | Be~ttwy taal 
~ aR 2 nL + 78, cos (sr) ‘ (22) 


s=1 


Taking proper account of the double-sign terms, including both 
terms under each double sign, the above reduces to 


EH 2He—* ts sin Bgt 


ay an nL B, cos (sr) (23) 
s=1 
The value of 8, given by (17) is 
B, = Vr? — cos) — a’. 
If we neglect a? compared with Tae — cos =) the above 


reduces to 


2 sin pa 
| 2 st 2n (24) 
Bs = 4/>>{ 1 — cos— ) = ———) 
LC n / LC 


and introducing this in (23), we get the following: 


i Be Ries eae sin B.t 
Un a nR + n VED ee d 


an ©8 (sr) 


_E , Ben (6N(=1)' sin gt 
~ nk ay n o> sin 37 ; (25) 
sel 2n 
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The first right-hand term is the steady-current component 
given by the impressed voltage divided by the total resistance 
of all the series elements. The summation term gives the transi- 
ent current consisting of several components, the number 
determined by the number of sections or meshes in the circuit 
system; each component an independent damped oscillatory 
current, all having the same damping but of different frequencies, 
which are given by 


fs Bs a 2n (26) 
; 2a rv/LC 


The highest frequency of the oscillations is when s = n 


fn = 


1 

rere } 

r/ LC 

the upper limit of the frequency range is independent of the 
number of sections used. The lower frequency limit is decreased 
as the number of sections used is increased. When the number 
of sections is very large, the lower frequency limit approaches 
zero value. This type of circuit is commonly called low pass 
filter. Actually, it is nothing more than a circuit system capable 
of responding resonantly to a number of frequencies whose upper 


1 
a/LC 
desired by increasing the number of sections and suitably select- 
ing the values of the inductance and capacity in the sections. 
The separation between adjacent frequencies is obviously less 
and less as the number of sections is increased. For very many 
sections, the system will be resonantly responsive to many 
frequencies and thus approach the condition of what is generally 
designated as band tuning. Theoretically, the condition of 
so-called band tuning is realized only when the number of sections 
is infinitely great; but practically, the condition is approximated 
with a reasonably large number of sections. 

Alternating E.m.f.—For applied alternating voltage, of sine 
wave form of frequency f = w/2z, the modified expansion formula 
(20) of Chap. II is to be applied, 


cy EES per 27 
AGN =jo ce yp 02) CQ 
ap 


» and the lower limit may be made as small as 


limit is 


44 HEAVISIDE’S ELECTRICAL CIRCUIT THEORY 


the first right-hand term giving the permanent-current, and the 
summation term the transient-current components. 

Consider first the summation term. The values of p and 
dZ(p)/dp have already been obtained (equations (16) and (20)); 
hence, 


= S co 


ed Ca! + jBs \e(—a 98, yt 
ae fe. erent jnt(1 + —— ja) (sr) 


Syn 


a aS (—a 4 9Bs )2¢ + jpt > 
+ j8.nL{(—a + jB;)? + w?} cos (sz) 
Taking both terms under each double sign, we have 


pee eee (—a + j8.)(cos Bat + 7 sin B.t) 
ae wis pee) jbsnL{(—a + jBs)? + w?} cos (sz) 


_ (a — jBs)?(cos Bt — j sin Bat) 
jbnL{(— a— jBs)? + w?} cos (sr) 
By simple algebraic transformation, it reduces to the following: 


eee eee 

ie ae ae 

e—at 2{ (a? ois Bs)? + w(a? — B2)} sin Bt — 4aB.w cos Bet 
Sa nLB,{ (a + 62)? + Qw(a? — B2) + w 4} cos (sz) 


For a = 0, it reduces to 


per B. sin Bt 
> 32 (p) > Le SPO on 
(pt + ot) 4D) 


Comparing (28) and (29) with (25), it is seen that the oscillation 
frequencies of the transient current are the same in both cases, 
independent of the character of the applied voltage. The ampli- 
tudes of the oscillations, however, are modified by and depend 
upon the frequency of the applied voltage. If the frequency of 
the applied voltage is the same as one of the frequencies of the 
free oscillations of the system; that is, w is equal to one of the 
values of 8, the amplitude of the transient component of this 


- (28) 
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frequency will be infinite when the resistance is of zero value, 
which is evident from (29). This is what we should expect, the 
system responding resonantly to this frequency. 
It remains yet to investigate the steady-state component: 
; Ee 
ji aeyee 
Z(D)p =iu 
2 Ee 
~ 22 8inh y sinh nyp=jo 


(30) 


We have the expression for cosh y 


Cp(Lp + ) 

Z 
Cjw(Ljw + R) 

ae 
From this it should be possible to express the denominator of (30) 
in terms of w. It is clear, however, that for n circuits, this would 
lead to an expression in w of the nth degree, a complicated expres- 
sion which would not throw any light on the problem. Of more 
interest is to consider for which values of w the system would be in 
resonance. The preceding discussion relating to the transient 
components would lead us to expect a resonance effect at the 
frequencies of the free oscillations of the system. This is actually 
the case, and it can be shown in this way. If we neglect the 
resistance, the denominator in (30) must have zero value for the 
resonance frequencies; we shall determine these frequencies and 
show that they are precisely the frequencies of the free oscillations 
of the system. 

By a well-known trigonometrical transformation, we may 

express the sine functions in a product series, thus: 


cosh y = 1+ 


in2 in2z in2 
sin v 1 sin” x sln* x 502 (31) 


sin nx = K sina/ 1— 
Pa 4 9 20 cor 
sin? — sin? — sin? — 

n n n 


n is an odd integer 


Put . . . . . 
x = jy; sin nz = sin jny = j sinh ny, 


| 


sin + =(sin. 7y = 7 sink’7,; 
and (31) transforms to 


sinh ny = K sinh y i a ls UG es (32) 


e 71 0 . 
sin? = sin? — 
n n 
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By (15), putting p = », and neglecting resistances, 
2\2 LCow? 
sinh? y = cosh? y — 1 = (1 + ie) em | = 10ur( as 1); 


Substituting in the first factor of the power series (32), we obtain 


YH APNge7e 
pe tee 7 tice 
1+ T= — 
sin? — sin? * 


The value of w for which this is zero is readily obtained by solving 
the above quadratic equation, which gives 


eae eee 
| n 


PN iaay Si ee 
n 


| | 


PO = 
as LC LG 
and 
T cma a: 
2(1 — cos = 2 sin (33) 
Qo = = eS 
J LC 1/LC 
In a similar way, the other factors of the power series give 
2 sin Ze 
J@2 = —< 
VLC (34) 
5. BN 
2 sin — 
Ws == ae 
VLC 


These are precisely the same as the values of 8,, equation (26), 
which establishes the fact that the frequencies for which the 
system is in resonance are the frequencies of the free oscillations 
of the system, as we should expect. This is but another way of 
arriving at the same result, determining the frequencies of the 
free vibrations of the circuit system. 

Series-elements Capacities; Shunt-elements Inductances.— 
A circuit system of several recurrent sections in which the series 
elements are capacities, and the shunt elements are inductances, 
is commonly designated high pass filter. A system of this kind 
responds resonantly to frequencies above a fixed lower limit 
determined by the choice of values of the inductances and capac- 
ities in each section, the upper frequency limit depending on the 
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Fig. 9. 
number of sections used. For this case, 


Cp (35) 
Equation (12), the expression for the current in the nth section, 
is perfectly general, applicable for all cases, irrespective of the 
choice of inductance-capacity combination in either 2; or 22. 
In developing, however, the full solution from (12) by the applica- 
tion of the expansion theorem, cognizance must be taken of the 
changed conditions. We have, then, for this case also, 

E 


‘n = 2 sinh sinh ny BEAD 
the determinantal equation as before, 
Z(p) = 22 sinh + sinh ny = 0, (13 bis) 
and 
ny = [srs s = 15 2.3.7 8H 
eee a (14 bis) 


The frequencies corresponding to the roots of the determinantal 
equation are obtained from (4), 


2 wal 
cosh y = 1+ on (4 bis) 
which for this circuit system gives 
Sr il 
cosh Va nCOS (*) =i1+ 2Cp(Lp +R) (36) 


The values of p which satisfy the above equation are readily 
obtained, as follows: 


ji) a ae eR 
where 


R B 1 R? 
On eae ae aA ye OT: 
2h aLc(1 — cos =) ae ep) 


n 
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We have, also, 


a = = z, sinh y cosh ny—— om). (38) 
By (4), 
: Oy _ 2Lp + R 
aren dp 2Cp(Lp + RK)?’ 
and 
) (2Lp + R) ( 3 5) 
Ag P 
z2 sinh y. a 
2 16, 
ap 2Cp*(Lp + R) cpi( 1 +4) 


Introducing this in (38) gives 


1 oe 
20) n( : ‘) 


wee 


a FnjBs cos (st) 
Cla 3 Joya Be jBs)* 


For p = 0 in this case, zp sinh y sinh ny = ©. The steady- 
current component is, therefore, zero, as we should expect, 
because the series elements are condensers. 

Substituting the values of p, and dZ(p)/dp from (37) and (40) 
in the expansion formula, we obtain the developed expression 
for the current in the nth section, 


e(— @EBs)t 
igs E> (—a + jB;)n(+jB;) 
G(r Ge) = ee Be) 


€#180t O(a? — 62) 
+njBs cos (sr) 


cos (sz) 


(40) 


) 
cos (sz) 


== He-@ 


(41) 


Taking both terms under each double sign, the above reduces to 


s=n 


in = —2Ee-+% C(a® — 62) sin Bet. (42) 
nB; cos (sz) 
s=1 
For a = 0, it reduces to 
CB, sin Bt 
mee Bae cos (sr) (43) 
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If we neglect the effect of the resistance in the expression for B, 


in (87), it reduces to 
1 


Bs a 
\ze(t — cos a) 
n 


1 


Figen ney a Ra (44) 
2 s1n on 1G: 


and introducing this in (43) gives 


eae C sin Bat 
on ES 7 (45) 


ai n sin 
n 


Comparing (25) and (45), it is seen that the transient currents 
in the two cases, low pass filter and high pass filter, are of exactly 
the same character, expressed by equations of the same identical 
form; the only difference is in the frequency characteristics of 
the two circuit systems. In this case, by (44), 


Pegi ge 1 


fs or 1 
4a sin 9, VLC 
For s = n, 
1 

Ue Teenie 

ior @ = Al (46) 
1 6 
fi= 


. 1 — 
47 sin anv EG | 
The frequency range is between these limits, the lowest fre- 
: 1 bet ee 
uency is =——= and the upper frequency limit is 
q Wy is i CG pp q yy 
1 
= ) 
4r sin ae /LC 


depending, therefore, on n; as n is increased, the upper 
frequency limit is raised, which may be made to approach 
infinity by making n, the number of sections, infinitely large. 
In any case, a circuit system of this kind is necessarily responsive 
resonantly to a number of frequencies confined within the range 
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between f; and f, given by (46), and by making n sufficiently 
large, an approximation to a uniform response to all frequencies 
within that range is obtained; and the condition for so-called 
band tuning for all frequencies above a certain lower frequency 
limit is approximately realized. 

The investigation was limited, in this case, to steady voltage 
which is believed to be quite sufficient to explain the behavior 
and characteristics of the system. The derivation of the expres- 
sion for applied alternating voltage is readily obtained by the 
second expansion formula, following the method used in connec- 
tion with the preceding case, the low pass filter. 

We have so far considered the simplest cases, those of a single 
reactance element, inductance, or capacity in either the series or 
shunt branches of the circuit system. It is quite obvious, how- 
ever, that a large variety of series shunt combinations of induc- 
tance and capacity, in either the series or shunt branches, is 
possible, each combination affording a distinct circuit system of 
definite characteristics of its own, particularly in the matter of 
the frequencies of the free vibration periods of the circuit system. 
It is not within the scope of this book to consider fully the theory 
of filter circuits; the discussion is given here more for the purpose 
of illustrating the application of the expansion theorem to the 
solution of problems in electric-circuit theory. We shall, how- 
ever, discuss briefly a few typical cases, confining the discussion 
to the determination of the frequencies of the free vibration 
periods, the governing characteristic in differentiating one circuit 
system from another. The complete solution in any particular 
case can be readily derived by the application of the expansion 
theorem in a manner outlined in connection with the derivation 
of the developed solution in the preceding cases. 

Series-elements Inductances, Shunt-elements Inductances, 
and Capacities in Parallel.—Disregarding resistances, we have, 
for this case, 


a: = lip, 


pose te ANSE, (47) 
1+ L2Cop? 


The frequencies of the free vibrations of the system are deter- 
mined, as in the previous cases, from the conditions: 


sinh ny = 0; 
ny = jst;s = 1, 2,3, °°: 4, 
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and 
oh hes 
Yeo at yi 
cosh y = 1+ nae 
uf 229 
which give for this arrangement 
LL 
Ae e) = 14570 + LiCw). (48) 


From this, the values of p are readily determined, 


a : |2E4{ cos ™ — 1 a2 
Ee L2C2| Li n 
a _ 1 
Ee N aio 


This determines the frequencies of the system, 


4De . g ST 
EL ASG By 
7 sin a ui (49) 


1 sr 
2 
fr = EA WEP : + ~ sin Qn. (50) 


The extreme frequencies at either end are given by 


a pt 
fi san +7, sin’ 5, S28, 
if = mee ene 41 
eo: Dra) Dine N Dy 
When n is very large, many sections, sin 27/2n approaches zero 
value, and f; approaches the value 1/2r~/L2C2. All the fre- 
quencies are then confined within the limits 1/2r~/L2C2 and 


J! + ar Obviously, by making L2/L: small, the 
1 


QnV/ L202 
frequency range is made small, and thus a condition of a 
circuit system capable of responding resonantly to a number of 
frequencies, all confined within a narrow range, is obtained, and 
the condition of so called band tuning for a narrow band of 
frequencies is thus approximately realized. 

Series-elements Capacities, Shunt-elements Inductances, and 
Capacities in Parallel—For this case, 
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3S} = = = 3 
S G5 s) aos) 42 = 
Fia. 10. 
1 Le 
a= Aj = oe 52 
4 Cyp ; il = LC op? ( ) 


The frequencies are determined as before from the relation 


cos = 1+ ad, 
n 2 22 
oe 1 + LC op? 
ome 21.C ip? oy 
which gives 
1 
Ps = joe 
afi2] 264( cos ies 1) _ 0,| 
Sree = (54) 
. aU 1 : 9 Bul 
iyfbace( — Ce sin = 
and the frequencies of the system are 
1 
Is aa — i (55) 
ae |LsC fee 
3 (Os 2n 
The extreme frequencies at either end are given by 
it 
Qn ned + 2% sine in) 
aor (56) 
Te ad ‘ 


Paa/ 2041 = :) 

Ce 
When n is very large, many sections, sin s1/2n approaches zero 
value, and f; approaches the value 1/2r-/L2C2. By making 
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C,/C2 small, the difference between the extreme frequencies may 
be made as small as desired, and thus confine the frequencies of 
free vibrations of the system within a narrow range. This 
arrangement therefore offers another suitable method for band 
tuning, the width of the band determined by the difference 
between the two extreme frequencies f; and fy. 

Series Elements: Inductance and Capacity in Series; Shunt 
Elements: Inductance and Capacity in Parallel—This arrange- 
ment gives 


a ds Lap 
Ae Ip + Cip’ 22 cee au T+ LiCop”’ (57) 
and 
i 
(tn + \e aR 1.0") 
il #1 Cip 
es = 1 is =1+ (58) 


2L 2p 


This transforms to the following: 


—4Lp sin? = = Inp + LiL2C2p* “iW abe ite ee 


or 
LyL2C\Cep4 + (u Cs + LC, + ADC, sin? * 9 mn) + 1 — 0. 
Solving for p? gives 


—(Iats Be Ont 41707 Bint sr) 


2 
Nine + EG, + 4D oC, sin? a) = 40, ,L0C1C2 


= » (59 
Ps OIFTAC sO: oo 
and the frequencies are given by 
(ze 4+ L1C1 + 4L2C; sin? ) + 
. LC + LC, + 412C, sin? “a — 40 ,L201C2 
(mee a . (60) 
Or 2L,L2C C2 
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SEAS 
If we make n very large, then for s = 1, sin? ae = 0, and 


yl 
Oe (L2C'2 -{- L4C;) nen (ORO i LC)? 


1 — 


V 2L1L2C1C2 
Fors =n (61) 
‘ia + 11C, + 4L201) + 
1 V (LeC2 +110; + 41201)? — 4LiL2CiC2 
nae Oe PITA GAGE 


The free periods of vibrations of the system have two frequency 
ranges, one between the values of f: and f,, when the positive 
signs are used in the double-sign terms of (60) and (61), and the 
other when the negative signs are used. When 110; = L2C2, 
the frequencies are given by 


2 St eu 
pe aft ttm + teers) 4 


2LC 
(62) 
cays Sek eos - 9 ST 
ee an + 2 SMa (1 + sin a 
Qa Ke | 


Magnetically Coupled Filter-circuit System.—The various 
types of filter-circuit systems considered so far are all of the 


uni 


Fiey ad, 


direct-coupled types, the recurrent sections connected directly in 
tandem. There is no need, however, to impose any limitation 
on the type of couplings to be employed between the recurrent 
sections. Results, similar to those obtained for the direct- 
coupled filter circuit, are obtainable with magnetically coupled 
circuits or mixed coupling. We shall investigate now the proper- 
ties and characteristics of a magnetically coupled filter-circuit 
system, that is, one consisting of a number of successively mag- 
netically coupled circuits, as shown diagrammatically in Fig. 11. 
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If we designate by z the impedance of each circuit, and by M 
the mutual inductance between any two adjacent circuits, we 
have the following system of equations for the current distribu- 
tion in the circuits: 

211 +e M pire = E, 
a3 + Mpro + Myr, = 0, 
ee alas on ai (63) 
etn = M pin—1 =F M ptm+1 = 0, 
Assume a solution of the following form: 
lm = Ae™y + Bey, (64) 
Substituting this in any equation of (63) except the first and last 
gives 
2(Aemy + Be) + Mp(Ae™Yy + Be“@—Dy) + Mp(Acmby + 
Be) = 0, 
and rearranging, 
(g + Mpey + Mpe7)(Ae™ + Be”) = 0. (65) 
The solution given by (64) will satisfy the equations of (63), 
except the first and last, if we assign a value to y to satisfy the 
relation. 
z+ Mp(er +7) = 0, 
or 
y= (66) 
cosh y = 5 Mp 
For the first and last equations we have 


z(Aevy + Be-7) + Mp(Ae? + Be”) = E, 


2(Aery + Be) + Mp(Ace™ Y7 + Bea FY (). (67) 


Two equations from which the constants A and B are deter- 
mined, thus: 
Eev(z + M pe’) 
en 7(2 + M pe’)? 7 ee D(z + M pe?) 
— Hev(z + Mpe™) ; 
eO-D7(z + Mper)? — D(z + Mpe v7)? 


) 


A= 


B= 
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Taking into account the relation given by (66), the above reduce 
to 


— Hem Dr 
ha Mp(e- oD — et D7) 
B Eer)y f ee 


Mp ot 7 — ety) 


Introducing these values in (64), we obtain the expression for the 
current in any of the circuits, say the mth circuit, 


E{ "Dy — ea} 


Im = Mp{e*rby— eMEDy} (69) 
Expressed in hyperbolic functions 
pe Bsn Gin Te by 
CRS Ay si (nee Dy (70) 
For the last, the nth circuit, 
; —F sinh 
tn eee (71) 


7 Mp sinh (n + ly 


The complete solution is readily developed from the above 
expression by the application of the expansion theorem. The 
determinantal equation is 


Z(p) = Mp sinh (n + 1)y = 0. (72) 
Hence, 
(n ae 1)y = jst, 
seiner Ta (73) 
rene. be a, n+ 1. 


The values s = 0 and s = n + 1 are to be disregarded, because 
for these values, the numerator in (71) is also zero. We take 
therefore, only the values of s from 1 to n. 

The values of p corresponding to the roots of equation (72) 
are obtained from the relation given by (66). If each circuit 
consists of an inductance and capacity in series, z = Lp + R + 
1/Cp, and by (66) 


1 
cos ei = (te “ : c) 
“ntl 2M p 
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which gives a second-degree equation in p, 


2LC + MC cos 


S71 9 = 
wal? + 2RCp + 1 = 0. 


Solving, we obtain the values of p corresponding to the roots of 
equation (72) as follows: 


IO (oe = rere (74) 
where 
—R 
ag at a ae A 
2( 1 + 2M cos 4) 
212 
Nie + 2MC cos 7 o 
B. = —. (75) 
AOR OTC oor 
We also have 
aie DREN cosh Gio Pn, (76) 
By (66), 
ey 2 
sinh yo os = oMp (2 a Gp) 
Substituting this value of dy/dp in (76), gives 
dZ(p) _ (n +1) 2 \ cos (n + 1)y. 
dp =——s«ip oie sinh a) 


Introducing the values of p, and 0Z(p)/dp from (74) and (77) in 
the expansion formula, and remembering that it is to be multi- 
plied by sinh y, the numerator of (71), we finally obtain the 
completely developed solution for the current in the nth circuit. 


4 2 sinh? rye( tas tiBs)t 
aes NN 6 eee cosh (n + 1) 
4 ( aC + js) : 


= 2ES 


sin? — 


(n + n(- + Oe 5 ay cos (sr) 


é( tastjBs)t 
ri (78) 
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Including both terms in each double-sign term, the above reduces 
by simple algebraic transformation to the following: 


Qa Bs 3 : Sir 
ast Sun 2 
: (8+ oatm) cos Bt — eam tin Bat} in ee | 
4Ras z 
(ir 1) re [F Cd +B or nese cos (sr) 
(79) 
By (75), however, we have 
af + B= -_ 
o(t + 2M cos ele 7 
Gee ipe, ee 
Making these substitutions, (79) simplifies to 
i= 
ce 2k cost —29,(L+2 M cos St, ) sims sin? oe 
a(z +2M cos — 
(n+1); 3R?+ a ss cos ( sr) 
(80) 


This is the complete solution for the current in the nth circuit of 
n successively coupled circuits. The current consists of n 
separate oscillations, each of an independent frequency, damping, 
and amplitude. The frequencies and damping factors are given 
by (75). Thus: 

For two circuits, n = 2. 


(eee _ —-R 
2L(1 + Ky’ So a Cao). 
foe ‘a nee 1 
2a /LOLE RK) Ww A/LCG =m) 
M 
1G ia (81) 


Ql = 
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For three circuits, n = 3. 


a1=-3 za 5 ag= = a3 = a : 
2n(1 + 1.414K)’ ie CHER se ae 
ne 1 1 1 
"On ECOG 4141) AT TE 1—1.414K) 
For four circuits, n = 4, 
of aah ; mes —R 
Dp Cie-e holon): 2L(1 + 0.616K)’ 
=f : =f 
ee QL (10616 R)” Ap GlP = eL.O1OR): 
eee 1 1 ae 1 1 
Penn TO FT Gi6R) © 1 247700 0.610K)’ 
ae 1 1 eo 1 
aor TC 0.610R). 27 4/L0G —1.616K) 
If we neglect the resistances of the circuits, equation (80) sim- 


plifies to the following: 
DEC: ibs sin Bt sx 
n+1 cos (sr) apy n+1 


(82) 


. 
Ue 


For two-coupled circuits, n = 2, the current in the secondary 
circuit is given by 


yy = _ — 6, sin Bit sin? @ + 62 sin Bet sin? (r) 
= a {82 sin Bot — B, sin Bit} - (83) 


which is in exact accord with (62) of Chap. II, the formula arrived 
at for this condition by another method. 
The corresponding dampings of the two oscillations of the 
two-couple circuits are 
ot eae a i —R 
Peon eee. SOT sek 


For a three-coupled circuit system, n = 3, the current in the 
third circuit is given by 


HO 


i3 = oP sin Bit — Bz sin Bot + , 5B sin Bit (84) 
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The corresponding dampings are 


fs lia <P tel ees Tk 
"OL AL AlaRY "OTR ee 


a1 


For a four-coupled circuit system, n = 4, the current in the 
fourth circuit is given by 


— 0.34558; sin Bit + 0.90468. sin Bot | | 
— 0.90468; sin Bst + 0.345584 sin Bat | 


and the corresponding dampings, 


a4 = 0.4EC (85) 


—k res -R 
te OE + L.616K)? DE OG1GR 
=- =F 
a3 


; SLi — 0.616K) = okt elon 


The method followed here made it possible to obtain the 
complete solution for the current distribution in any number of 
coupled circuits; to determine the frequencies, dampings, and 
amplitudes of the several oscillatory components. The only 
limitation imposed is that the electrical constants of the individ- 
ual circuits L, C, and R, and the mutual inductance between 
any two circuits are the same for all the circuits. 

Resonance Effects of Multiperiodic Circuit System.—We have 
established, now, that a circuit system consisting of a number of 
coupled circuits has a number of periods of free oscillations, each 
oscillation characterized by an independent frequency, damping, 
and amplitude. This is perfectly general, applying to all types 
of couplings, magnetic, static, or direct. We may reasonably 
infer from this that if we should impress alternating voltages on a 
circuit system of this kind, it would respond resonantly to several 
frequencies, which are precisely the frequencies of the free 
oscillations of the system, and to those frequencies only. Hence, 
a resonance curve of a circuit system of this kind would, in effect, 
be the resultant of several resonance curves, each corresponding 
to one of the frequencies of the free oscillations of the system. If 
the circuit system is designed for small frequency separation 
between adjacent frequencies, then, obviously, parts of adjacent 
resonance curves would overlap, the curves would partly merge, 
and the effect would be as if the system were in partial resonance 
even for the frequencies which are between resonance frequencies. 
If the circuit system is designed to have a large number of 
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resonance frequencies, all confined within a narrow range of 
frequencies, the resonance curves would be crowded together 
and the effect approximate a uniform resonance for all the 
frequencies within that range of Be This is sometimes 
designated as band tuning. 

The author has emphasized this point, at the expense, perhaps, 
of repetition, because of the confusion in the literature of this 
subject regarding the meaning of so-called band tuning. Further 
to bring out this point, I have worked out a numerical 
problem to show the resonance effects for impressed voltages of 
different frequencies in a two, three, and four magnetically 
coupled circuit, and plotted the resonance curves which are 
given in Figs. 12, 13, and 14. 

To compute these curves, use is made of the basic formula 


E sinh y : 
Mp sinh (n + 1)y 


in = 


For impressed alternating voltages, this formula, with the change 
of p to jw, gives the steady-state current. We also have the 
auxiliary formula 


—z 
cosh = 2Mp 
and for p = jw, 
-(u-8) 
hy = e ) 
wad 2Mw 


neglecting resistances. 

It is clear that, since the frequencies of the free oscillations 
were determined from the determinant equation sinh (n + 1)y 
= 0, then obviously reversing the process, impressing voltages 
of these frequencies on the system would also give sinh (n + 1)y 
= 0, and the currents for these frequencies would be infinitely 
large if the resistances are neglected. This can be readily verified 
working backward. Substituting the expression for 8, defining 
the frequencies of the free oscillations, in place of w, we have 


sinh y = Vcosh? y — 1 


= | e = te) - (86) 
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Substituting for w the free oscillation frequency value 


ee aA (87) 


St 
s{ue( +2H cos rai :) 


the above transforms to 


sinh y = 
rs L 2M Si 4M? 
eee _ yaks bi) xc =) 
Ge; 7 098 n+l ea 
2M 
2M Sqr 
Nisa eee cos) (88) 
which reduces to the following: 
; SLE Sa RN nL 
sinh y = J- 1 + cos ae Pe aad 
Hence, 
enSE 
es eat 
and 
sinh (n + 1)y = sin (sr) = 0. (89) 


For these frequencies, the current in the circuit is infinite for 
zero resistances. ‘To determine the general form of the current 
curve in response to voltages of different frequencies, it will be 
sufficient to calculate the current values for frequencies in 
between resonance frequencies and for frequencies on either side 
of the two extreme resonance frequencies. The following values 
of the constants of the circuits were assumed: 


L = 0.5 mh; C = 0.0002 mf; M = 0.05 mh. 
A two-coupled circuit system: 


The resonance frequencies are 1/27 107/+/11, and 1/2 107/+/9. 
For the frequency in between these frequencies, say 1/27 


107/+/10, cosh w = 0, and sinh y = j; y = jr/4 and sinh (n + 
1)y = sinh he 9 


For yA = 


sinh y = 1.732 and sinh 3y = 25.96. 


mvt 


For fs = se sinh y = 1.732 and sinh 3y = 25.96. 
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A three-coupled circuit system: 
The resonance frequencies are 


Fg eee 100 Je ns, 107 
a a —— area 
Qn /11414"*” — Inr/10"** A \/8.586 
107 : : : : 
For js = am sinh y = 70.937 and sinh 4y = 70.49. 
107 é : : : 
fo= an x 3.05" sinh y = j0.937 and sinh 4y = 70.49. 
; 10’ : . 
fo= an X35) sinh y = 0.515 and sinh 4y = 2.094. 
107 : : = 
f= On X 2.75" sinh y = 0.698 and sinh 4y = 3.454. 
amet 
| \/ ' 
= 
< 


eee AD, Soe So S| 


Frequency 
Fie. 12. 


A four-coupled circuit system: 
The resonance frequencies are: 


ee Oe ee 
mo (Gis. | rN 10.618 9.882 
107 
Wrv/ 8.382 
7 
For fs = Ag Shih op = FURS e oH ey eee ee 
T 
vf 
For fs = Stes SA sinh y = 71.0; sinh 5y = —71.00 


Qa 
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7 
POR fp Se eo minty a OIE? wine oe 


Qn 
7 
For fs = es sinh y = 0.95; sinh 5y = 34.00. 
T 
a 
For fy = oar eee sinh y = 0.55; sinh 5y = 6.9. 


From these values, the graphs of Figs. 12, 13, and 14, are plotted 


Amplitude 


Frequency 
Fie. 13. 


Amplitude 


Frequency 
Fig, 14. 


CHAPTER IV 
OCEAN CABLES 


CIRCUITS OF DISTRIBUTED CAPACITY AND RESISTANCE 


A considerable portion of the second volume of Heaviside’s 
“Electromagnetic Theory” is devoted to a discussion of problems 
relating to wave propagation on cables. A large number of 
problems touching various phases of the subject are analyzed, 
in which direct operational solutions and the expansion theorem 
are applied. It is altogether beyond the scope of this book to 
give as comprehensive a discussion here; nevertheless, the topics 
selected are such as to cover the most important phases of the 
subject and at the same time show the application of the opera- 
tional solutions in a variety of forms to serve as a guide in the 
solution of practical problems that are apt to arise in connection 
with this subject. 

For a cable of negligible inductance and leakage, the relations 
between the current and voltage at any point on the cable, 
distance x from the transmission end, are given by the following 
equations: 


Ri = —2 
di (1) 
CpV = i 


R and C are the resistance and capacity per unit length of cable, 
and p = d/dt. 

From these equations, we readily derive the propagation 
equations for either the voltage or the current, as follows: 


PV 
RCpV = FEE 

d1 (2) 
RCpi = de 


It is sufficient if a solution is obtained for either one of the 


equation (2); the other is readily evaluated by the aid of the 
65 
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relations (1). Consider the voltage equation and assume a 


solution, 
V = A&* + Be-™. (3) 


This satisfies the differential equation if 
K = VRCp. (4) 


A and B are independent of x, but functions of time, and depend 
on the terminal conditions. By either of equations (1), we get 
the following expression for the current: 


pie Bee (5) 
R 
These are general solutions involving indeterminate time func- 
tions A and B, which must be separately evaluated in each 
special case. 


CABLE OF INFINITE LENGTH 


In the case of a uniform cable of infinite length, any electric 
impulse impressed on the cable travels along the cable without 
suffering any reflection at any point on the cable, and it is 
because of the absence of any 
reflections that the mathematical 
solution of the problem is con- 
siderably simplified. We shall, 
therefore, begin the investigation 
of this subject by considering some 
problems relating to cables of 
infinite length. 

Assume a steady voltage EH impressed at the transmitting 
end of the cable x = 0, to determine the voltage and current at 
any point on the cable distance x from the transmitting end. 
The terminal conditions for this case are 


Hie. 15. 


x 0; V=8#, 
z=w0;V=0, 
which give 
A'=Qand B= #. 
Hence, 


V = Be, 

ees (6) 
= tee fi tl) 

V R € . 
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At the transmission end, ¢ = 


Ko = ) 


te n,/°P. (7) 


What meaning are we to attach to the fractional derivative 

Vp = Vd/dt which appears in above equations? Heaviside 
asserts, and apparently he discovered it experimentally, that if 
we operate by Vp on unity function, that is, a function which 
has zero value for ¢ < 0 and is constant for all values of t > 0, the 
result is 1/1/7rt. That is, 


1 
Vl ae (8) 
There is no question about the correctness of this relation; 
many problems worked out by Heaviside on this assumption 
yield correct results. Farther on (see pages 73, 97), we shall 
give partial proof establishing the validity of this assumption. 
With the definition of ~/p1 as given by (8), the current at the 
transmitting end is given by 
2 CG 
io = EH =PE (9) 
Obviously, operating by ~/p on any constant quantity gives 
1/+/zxt times that constant. Hereafter, if not otherwise specified, 
unity function is to be understood. 
To obtain expressions for the voltage and current at any 
point on the cable, distance x from the transmitting end, we 
must algebrize equations (6) which involve operating by ¢«~*’. 


Expand e-*? in a series, and 


( K2y? K3y3 K4z4 
YS a | aah ae 
RCp)x2 (RC p)*x8 
= E}1—(RCp) x + | pe - Be = +--+}. (10) 


The terms of integral powers of p in the above expression are 
to be disregarded, because they imply complete differentiation 
of a constant, in this case unity, which gives, of course, zero. 
Hence, (10) reduces to 


343 COp)25 
“yas TEN eo UO oo ae | (11) 


3! 5! 
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Since 


Sar 
Vp = ar 


we get by successive differentiation 


- 1 
% = or/p = ———i-% 
ws oo Wa 


; 3 
Y= p/p = ———__1-¥% | (12) 
p 1 ae FDA 
a 3°5 7 
%— 3 Sey ee eet 
pt = pvp Ree 


Introducing these values in (11) gives 


t 7 BCNA 1 RC\* 
(gop itn el eee ey A 
| AT) *trawdt)? 


3 RC 2 5 ear ec 
Pa t ) pee 
or 


V RCx*\% L/RCx# 1 /RC2\2 
pra) esas) vec 4) a 


} 


which is the expression for the voltage at any point on the cable. 
It is seen that it is a function of RCz?/4t, and if we put this for 
brevity y?, we have 
Ve 2 | UE id An 
Hit rae aa a inera 


eo] 


which may be put in the form of a definite integral. 


V Oa 
—=1]———]| e-*dy. 
E wa y. . (14) 


This integral term is sometimes designated as the error function 
and denoted by erf. y for which tables have been worked out. 
A short table is given in the appendix. 

The expression for the current in the cable could be obtained in 
the same way, utilizing the same operations as in the derivation 
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of the voltage formula (13). It is simpler, however, to derive the 
expression for the current from (13) by the circuital relations (1). 


on 3 
~ Rade 
EHU/HO\s 1 /RO\# 88 (ROK, | 
Be l\ xt 2-2~/n\ t - OBI Faye f 
Bea ERG PRON, 
Bt) See tag ie 
7 a oe 
StS ae, (15) 


Equations (13) or (14) and (15) determine the voltage and current 
at any point on the cable for any time ¢. 


CABLE OF FINITE LENGTH 

In a cable of finite length, the character of the voltage and 
current depends not only on the electrical characteristics of the 
cable itself but also on the terminal conditions at the ends of the 
cable, the latter determining the character of the reflections at 
either end. The voltage and current waves on the cable are, in 
fact, composite waves, the resultant of practically an infinite 
series of waves arising from the repeated reflections at either end. 

We shall consider first the simplest cases, those of a free cable, 
the ends either grounded or open, no terminal impedances, in 
which case the reflected waves are of the same character as the 
incident waves with or without reversal of signs. For a grounded 
end, the voltage wave is reflected negatively; for an open end, 
the voltage wave is reflected positively. 

Whatever the terminal conditions are, the equations for the 
voltage and current in the cable are still given by (8) and (5) 


Vi Ale Be =", 


i ire + Be*), Ce 


The factors A and B are time functions, but independent of z, 
and are determinable from the terminal conditions. 
A Cable of Length 1 Grounded at Both Ends.—A constant volt- 
age E applied at the end x = 0. We have for this case, 
z=0;V=H, 
ap = Is 6 0p 
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which gives on substitution in (16), 
A+B=E&, 
Ack! + Be* = Q, 
from which A and B are determined as follows: 
— Bex! 
is e- Kl’ 
Hex! (17) 


eK! —_ eK! 


Fie. 16.° 


Introducing these values in (16) gives 


JB fekt—=) zie eK (l—2) } 
eK! = eK! ’ 


2 Ee fe) a8 pias (18) 


eX! — ¢-Kl 


VY = 


These equations may be put in this form: 

E sinh K(l — 2) 
sinh Kl 4 
Cp cosh UG eee 
R sinh Kl 


To develop the solution from the above equations, we may utilize 
either the direct operational process or apply the expansion 
theorem. We shall adopt the latter as being the simpler method. 
Heaviside carried through in detail the operational process. 

The expansion-theorem formula, of which the derivation is 
given in Chap. II, is 


lee 
(19) 


~=H 


Eat 


1 OF Va= Tor Si soe -aZ(0) . (20) 
n=1 dp P = Pn, 
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which may be applied for the determination of either V or 7. 


For the determination of V, we note that the determinantal 
equation is 


sinh Kl 
OD eran =in). 7S aN 
Hence, 
Kb = gna; n= 1, 2,.3,4 °°: 
and 


K = VRCp = ae 
from which we get 
nn? - 


Po TROP 22) 


There are an infinite number of values of p given by equation 
(22) which satisfy the determinantal equation. We also have 


Te = cosh xi, 
but 
CED IRC RCL” RCP 
Op 2N p 2/RCp gna’ 
therefore, 
dZ(p) _ RCP _ RCP 
mer = one cosh Kl = ie COS Nr. (23) 


For p = 0, we note that sinh Kl = 0, but so is also the numerator 
in (19), sinh K(l — x) = 0, which would make the expression 
indeterminate. By taking, however, the value of sinh K(/ — z) 
/ sinh Kl as p approaches zero, we get 


sinh K(I — x) a K(l — x) oy eee 
sinh Kl Kl t? 
and 
Z(p)p-0 = 1 — 7 (24) 
Also, 


sinh K(l — x) =jsin at —z) =jsin ne(1 = ") 


: PR Ree 
— j COs nr sin ——: 


i 


Il 
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Introducing these values in (20) we obtain 


= Ror! - Ne 
ae j cos nm sin ~~ 
£ 
i a (: re , +BD a ee 
AE ROR Qjinr 
n= 0 Sn 
z Oi Loe ae. Bee 
-2(1-7)-25,, a Rs (25) 


n=1 


which is the complete solution for the voltage at any point on the 
cable for any time ¢ after the application of the voltage. For 
infinite time, we get the steady-state condition which is, by 


(25), fort = o, 
V= B(1 — ?), (26) 


The expression for the current may be derived in the same way 
by the application of the expansion theorem to the second equa- 
tion of (19), or we may obtain it from (25) by the circuital 
relation 


; dY 
as dx 
which gives 
Nn=0 ny? 
: E 2E naz RCP 
i=a +p cos 7 € . (27) 


n=1 


If 7 is made infinitely long in above expression, the first right- 
hand term reduces to zero and the summation terms go by steps 
of x/l, and for l infinite, 7/l is infinitely small. If we put nz/l =« 
s, then 7/l = ds and the above becomes 


foo} 


g2 
Re 


oS cos (sx)ds. (28) 


= ae € 

For « = 0 at the transmitting end, this reduces to 
foo} $2 

BRO RG 


asd {| ds. 


OCEAN CABLES 73 


Put 


and 


This is a well-known definite integral whose value is ~/7/2 (see 
Osgood’s ‘‘Advanced Calculus.” 
Hence, 

_ 2EW/RC Va ee! 30 

ee) NR ce 

Comparing this with the second equation of (7), which is also 
an expression for the current at the transmitting end of an 
infinite cable, we get the relation 


Ne iz ee Vat 


ae 

Vp a: (30) 

which establishes, at least in an experimental way, the validity 

of the assumption concerning the fractional differential operator. 
Another proof is given in Chap. V. 

Cable Open at One End.—Another simple case is that of a 

cable open at the far end. The terminal conditions for this case 

are 


and 


x=0;V =H, 
Daa er =, 
which give, by (16), 
A B — EH; 
AeA: aes =a) (3b) 
Hence, 
Eek! 
Jt ar alae = 
— Kl 
a 2 
B= 
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Introducing these values in (16), we obtain the equations for 
the voltage and current as follows: 


cosh K(l — x) 
Lees cosh KI’ i 
op Cpsinh J K(l — a) co 
i R cosh Kl | 


Applying. the expansion theorem, we have, for this case, the 
determinantal equation 


Z(p) = cosh Kt = 0; (34) 
which gives 


Kl = jy5n = 1,3, 5,7, 
and 
K= VRCp =F 57° 


From this we get 


n? 2 
= — ROP (35) 
Also, 
OZ(p) _ om gees 
op Op 
but 
O(Kl) _ RCP 
Op so jnr 
hence, 
Pe aa 
OF) BOT sin Kt = BCE ain (36) 
Dp jinn * na40 2 
For 


p = 0; cosh K(l — x) = 1 = cosh Kl. 


Introducing these values in the expansion formula, we obtain 
the developed solution for the voltage as follows: 


(1-2) “akc 

cos. ae 

_n_ = 2 l 

Vik E> nen? RCP . nr : 
ARCE tage 2 


n%2 


i= em : ne "arc (37) 
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The expression for the current can be readily derived from (37) 
by the relation = —1/R dV /dz. 
Thus: 


N22 
t 


as ~ 4RCL ; (38) 
1 Sy by cos. ——F 


To show the application of the second formula of the expansion 
theorem, derived for the case of applied alternating voltage, 
we shall consider the preceding problem, assuming an alternating 
voltage of sine wave form and of frequency f = w/2zm applied at 
the transmitting end. For this case, we have the formula 


Heit t 


(5 ey sae ; (39) 
ZiDynie Di joa PP 


only the real part of the above to be taken. 

The determinantal equation is the same as in the preceding 
problem of applied steady voltage, and we have, therefore, the 
same values of p and 0Z(p)/dp. The summation term is to be 
modified only by replacing p, by p, — jw. For the first right- 
hand term, we have 


1 _ cosh K(U—2) _ cosh VjwRO(— 2) 
Z(P)p=iw cosh Kl cosh VW jwRCI 
cosh ee 417) (2) 
= cosh sade +] i 
Hence, 
cosh | YB 4 na 2) 
= Eeet 
pe cosh SAD iy iu : 


BN Pettcosh Ka) 40) 
= 2 1 2/02) _ 
(p2 + w’) Op Pp = Pn 
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Introducing the values of p, and 0Z(p)/dp from (35) and (36), 
equation (40) takes the form 


cosh a + 7)(l — x) 
V fast) EB jut —— — 
; Wh) ag 
cosh aa a (1+ ))l 
NEE (LE — {aint 
Rr n® sin —~e 
IRCP OG niet Soe 
16RCcR T° 


For w = 0 (41) reduces to (37), as it should. 

Voltage Applied at Intermediate Point on Cable.—We shall now 
consider a more comprehensive case, that of applied voltage HL 
at an intermediate point on the cable. We have, in this case, 
two series of waves, the primary wave is doubled, going both 


ways; each wave suffers reflection at both ends of the cable. 
Count distance from the end to the left of Z the voltage applied 
at the point x = 1,, and the total length of the cableis/. Desig- 
nate by Vi, 71 and Ve, 22 the voltage and current to the right and 
left of H, respectively. 

For the right of EZ, 


Vi = AyeKO OVD Bye KO 
Fae = [Ase 4 Bye 1D}, (42) 
For the left of Z, 
Vo = Agek? + Boe-*?, 
a pds Bye-). (43) 


The distances to the right and left of Z are counted in opposite 
directions, hence, the reversal of sign in 72. 
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For a cable open at both ends, we have the conditions, when 


r= i 11 — 0, 
Hence, 
—AreKO'D + Bye&') = 0, 
As — Bg — 0, 
and 


ee t—1 
B, = A 362K ( Dv, 


Be i Salo. cee 


At x = 1, the point of voltage application, the following condi- 
tions must hold: 


which give the following: 


—Ai + By = —Age®t + Boe*1, 
Ai + B, ==> A oe*41 = Boe = 10), 


Introducing the values of B, and Bz from (44), we get 


A ane Ot) Ae Ane A: 0, 


Ay + AyeXO ID — Agek1 — Age *1 = EH, 
or 
AyfeKC) — KCl Vet) $e Aa {eR — eK} = 0, 
AyfeRUID 4 eK) beK-) — Aafekli 4 eK} = E 


) 


which, on solving, give the following values for A; and Ag: 
Ai = 
HeK—p { eK pele e Ky } 


f eKOt) — Kip} f Kl f eK} 4 {KC FRED) fk — Rh} 


_ Be-Ky sinh Kl 

i; 2 sinh Kl 

4, = DH sinh KW = bi), 
eas sinh K1 


) 


(45) 
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Substituting these values in (42) and (43), we obtain, by simple 
transformation, the following expressions for Vi; and V2, the 
voltages to the right and left of the point of the applied voltage: 


Thus: 
Y,= eee cosh K(I — 2), 
sinh Kl (46) 
Vow —E sinh K(l — h) ah 
a=. sinh Kl Seg 


In a similar way, we can obtain expressions for the voltages on 
the cable when both ends are grounded or when one end is 
grounded and the other end open—two cases. 

The two ends grounded: 


_ EcoshKl, .. 
ase Ti Ci pag (47) 
“te —E cosh K(I —h) cahoie 
a= ah KE aan 


The end to the left grounded, and the end to the right open: 


) 
Vi= Hee ned — 2), 


cosh Kl re 
V _ —E sinh K(l —1|) . hK 48) 
an cosh Kl a 


The end to the right grounded, and the end to the left open: 


, ., 2 sinh Khe. 

Y= Tag, sinh K(l — 2), : 

Vee —Ecosh K(I — 1) hK ea) 
2= FO cosh Kz. 


Expressions for the currents can be readily derived from the 
voltage expressions by the relation 
1 dV 


{= ———— 


R dx 


For either of the above cases, the complete solution is readily 
developed by the application of the expansion theorem. It will 
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be sufficient to consider only one case, and we shall take that 
of cable open at both ends for which formulas (46) apply. The 
determinantal equation is 

Z(p) = sinh Kl = 0, 

Tm jiig a= 112535 

The values of p corresponding to the roots of this equation are 
given by 

K = VRCp = i 


and 
— 2772 
Po = RO (50) 
2 
EA osh (glee = ees cosh Kl, 
cs) Op 2ynr 
2 
- aos cos (nz). (51) 
For p = 0, 


cosh K(l — x) = 1 
sinh Kl, od Kl, ly 


sinh KL KL 1 


Introducing these values in the expansion-theorem formula, we 
get 


a =, 
Bl j sin ™™" cos ro(2 a i) es 
sy : 
Ue are mak E>, ue) Pee) (nm) 
= RCE 2ynr 


which simplifies to the following: 
4 OE 1. nal nme — Wade, 
Vi=E = n —7~ cos 1 € BOS, 


Put nr/l = s,, and the above takes the form, 


32 
wan 
Wire 


= 4 OR 2 (sl,) He (sa)e ; (52) 


For Vz we use the same values of p and 0Z(p)/dp, but in this 


cops raspy tl eet Ce a 
case, tor p = JY, Z(p) p= i Kl 1 
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ces ly nex 
1 j SID Nar 1 —- Ze COs Mae ae 
1 Cl 
a irae (2 is t) +E 2 a? ROP Ai we 
RCE 2jnr 


Hence, 


cos (nm) 


r 1 NTL 
s1n n> COS sy ine nur? 


I 
J ly l —Rot 
2 w(1 a ') is n>, = ¢ Ral, 


and for nz/l = s, 


Ve —H(1 A 2) ae pea (sl 2 — (sx) - act (53) 


For cable grounded at both ends, formula hi the determinantal 
equation is the same as in the previous case. Hence, for this 
case, also, 


—nr? dZ(p) _ RCP 
[Op “RCE ) dp = Qinm cos (nr). 
For p = 0, we have, for V1, 
eet Siege 
Z(P)p=0 l 
aele Ae 
Z(P) p=0 l 
Combining and simplifying, we obtain 


and for Ve, 


- cos site sin A __ next, 
Vie a(1 a *) - 22> —- ¢ Rat", 


c. x oe *) - Tipe oe = (sx) - a (54) 


a a SOS cos 2S ae (sx) - Fe (55) 


The complete solutions for the other two cases, (48) and (49), 
can be developed in the same way. 

Cable with Terminal Impedances.—We have so far confined 
the discussion to a free cable, that is, a cable with free ends, 
either grounded or open, in which case the reflected waves are 
copies of the incident waves, with or without reversal. The 
series of waves resulting from repeated reflections are all of the 


and 
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same character, which makes it possible to express the solutions 
in a comparatively simple form. When, however, impedances 
are introduced at the terminals, the reflections are no longer 
complete, and the problem assumes a more complicated char- 
acter. In some cases, however, it is possible, even for the condi- 
tion of impedances at the terminals, to obtain solutions either 
by direct operational process or by the application of the expan- 
sion theorem. 

We shall consider first the case of an infinite cable, an imped- 
ance Z, at the transmitting end. The terminal condition for this 
case 1s 
By (3) and (5), 


B+ Ds/PB = B, 


and 
ae 
1+ 2u|? 
Hence, 
V = ee ere (57) 
1+ Zon) P 


We shall investigate the input voltage, that is, the voltage at 
the transmitting end, for different types of impedances. 


Fia. 18. 


Resistance at Terminal; Z) = Ro.—For this case, 


By division, 


f Cone ee Dye (cy a fp 
¥e= B\1~ n°?) +H? pi SPY" 4 (58) 


82 HEAVISIDE’S ELECTRICAL CIRCUIT THEORY 


The terms of even powers of p imply complete differentiation of 
unity, which would give zero, of course, and are, therefore, to 
be disregarded, and the above reduces to 


Y 2 ) % 
Vy = B — HR] + SP) + nf Sa @) (59) 


Operating by p” on unity gives 1/~/rt by (8); hence, 


4 ( Y TANG 1 
Vo=HE — zr) io + me) + ri) + +--> ee (60) 


Operating on 1/+/t by p, p?, etc., involves repeated differentia- 
tion, which can be done at sight (see (12)) and we get the com- 
plete solution 


CN" ft Es) Ge ‘) 


This solution is in the form of a series in descending powers of ft, 
suitable for computation for large values of t. For very small 
values of t, the above formula is obviously not well adapted. 
It is possible, however, to develop (58) in the form of a series of 
ascending powers of ¢ convenient for computation when ¢ is 
small. 

Write equation (58) in this form: 


= ee (62) 


By division, 


RB \e R R \34 1 a 
Veer ef pa ee (ee one eee ee 
(2 +(e) (ier ») i, (race) 
This may be put in the following form: 
R Te eaNe R \3 1 ( RRCp\* 
ames eee eee 0 
: BI aes 3 (Face) st (aac) as i( p ) 
Qie Vig 
By) Crater ‘ eae 
"Rice t (rcp) i i 


In the first series of the above equation, the operation required 


is successive integrations of p’ = 1/+/zt; and in the second 
series, the operation required is successive integration of unity. 
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In either case, the operations can be performed at sight, bearing 
in mind that 


1 4 Id rp 

P 

1 WA 2a 2, 
A arcs 
| ees LAO Tay rg: 
= = ——_____}{” 
pe bois aor 


Performing these operations, we get the following: 


2R QR \? t% OR \’ v5 i RiC\ 
reap Raw ead Uap ee ee Pe ee Pte OS Aes 
Vo {Pa +(e) +R) pest Ges 


i 2 3 3 \ 
| Rt R P R t 
Bs ype LS) LT Ls oe ee ar 
Vac t a) as (ea) ez | 
and this simplifies to 
ee. BEY 3 QR t 1/2. 
Say on ei) 11 + 7203 T 3. (Re) ES a6 


AG _ civ), (63) 


a very convenient formula for calculating Vo for small values of 
t. Formulas (61) and (63) supplement each other, one to be 
used for large values of t and the other for small values of ¢. 

Terminal Impedance z, = 1/cop.—Substituting 1/Cop for Zo 
in (57), we get the expression for V, for this condition, 


Fig. 19. 
y.= E = E 
ae 1 [Cp g 

pel ER a oe 

Cop\N R C2Rp 
Put for brevity C/C?R = 1/a, and 

ee e 
1+ > 


Vap 
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By division, the above expands in the series 


i 1 1 bee 
vi = 2)1 ~ (apy * (ap) ~ Can) T 
1 
=F ts nwo * Ge * | 
1 
— Bit 4s ucheC meen: 


The first bracket term of (65) involves complete integration 
of unity-time function which can be effected at sight. Thus: 


a ae W chyna. 

ae: pe te dee eee 
t 

= Ee. (66) 


To evaluate the second bracket term of (65) involves integra- 
tions of the factor 1/(ap)”. 


Now, 
As 1 i Sle Is Ses 
pep PVrt Wr 
eam 2 gee 
eee Ey ees 
ppt /x J, ae 
ia 4 ad i: oe ey: 
z= = | dt = = 
pp Al 3° 5SW/z 
Hence, 
1 1 1 1 
Bute cas eae Hite te eee 
i oP * Gp? * (yt | (ap) 


t \7 1 fh 
are) [t+ae) teaa) + - | 
Combining the results of (66) and (67), we obtain the complete 
developed solution for the voltage at the transmitting end, 


t t \% 1 2¢ L/2i7 
Vo = Ee — 20(-,) oe ee (68) 
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a very convenient formula for computation when ¢ is small, the 
series rapidly convergent. Heaviside gave an alternate formula 
for this case, adapted for large values of t, which is as follows:! 


al ONG a Re BN 
rae) Bt sl (51) ae 


CABLE OF FINITE LENGTH; TERMINAL IMPEDANCES 


(69) 


When we depart from the ideal condition of a free cable, the 
ends either open or grounded, we run into considerable mathe- 
matical difficulties. The problem becomes more difficult because 
of the fact that the reflections at the ends are no longer complete; 
part of the wave energy is absorbed by the terminal apparatus, 
and only part of it is reflected, with the result that at each reflec- 
tion the entire character of the wave is changed. The final 
solution must, of course, sum up all the repeated reflections at 
each end, which introduces considerable complexity in any 
attempt to obtain a general solution. In some special cases, 
however, it is possible to obtain a solution in a form suitable for 
numerical calculations. 

It would be well worth while for the student who desires to 
obtain a mental picture of the physical process in the building 
up of wave solutions to read Heaviside’s ‘Electromagnetic 
Theory,” Vol. II, pp. 67 to 77. He sets forth clearly the whole 
process, starting out with a single wave generated at one end, 
taking proper account of its attenuation as it reaches the other 
end, multiplying by a coefficient of reflection, and transmitted 
back on the line; again attenuated and reflected, the process 
repeated over and over. Summing all up, he arrives in a physical 
way at the complete expression for the voltage on the cable. 

Let us assume a cable of length /; impedances z and 2, con- 
nected at the transmitting and receiving ends, respectively. 
We have the general expressions for the voltage and current, 
equations (3) and (5), which are applicable for all cases, irrespec- 
tive of the terminal conditions. 


Vo= Ac + Be-** 
7 = A (det |) Be), (70) 
K = VRCp. 


1See Hnavisipk, ‘Electromagnetic Theory,” Vol. 11, p. 42. 
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A and B are determined from the terminal conditions, which are 
as follows: 


Fors = 0; V, = #— 2 


71 
For’ =1:V; = 24r 


These give, on substitution in (69), 
ak 
R 

. iy 21S - 
Ac! + Be iL “p | —Ack -+ Be zy 


A(1 — 7) +8(1 +) ue 

12 
Ac! ae Be 2K a | ; 
€ (+ e+ € (1 - z=. | 


From these two equations, expressions for A and B are readily 


A+B=E-—**(-A +B), 


Rearranging, 


obtained: 
(1 - ae ] 
R 
as Cs K 
(: + 2,2 7) (eK! — ¢Kl) — R (2: + z,) (e 4+ e-*) 
(73) 
(1 z Ss ) kl 
ee KR? eK 
(1 + 2,2; =) (=! — €*) ari (2, + 2,) (eR! + e*’) 


Introducing these values of A and B into (69), we obtain the 
following expressions for the voltage and current in the cable: 


(eX¢-2) 3) ) 42K (exa-2 +_Kt-®)) 
R 
V=E ; ) 
i Kl —Kl K 1 1 
1+2,2; R? (eX! — e-Kl) + Ret 2,) (eX! +--*?) 


a (74) 
(ek U-2) +_eK-+)) +2rp (ekK¢—2) —¢K¢-2)) 


i=E— 
R Kk? 
ait t eret De ) (ee) +e +2,) (eX! -®") 
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Expressed in hyperbolic functions, 


: ics 
sinh K(l — x#) + 21 p cosh K(l — x) 
V i HK k2 ; K d 
1+ eee pe sinh Kl + R (a + 2,) cosh Kl 


K 
i cosh K(l — x) + Zp sinh K(l — x) 
Foye 


R TaN K 
(1 + ape) sinh Al + Jp (2 + 2,) cosh Kl 


For x = l, at the receiving end, 

; ” Ez, 
(& + 2a) sinh Kl + (2, + z2,) cosh Kl 
ee [ @) 

(& + 2,2, 7) sinh Kl + (z, + 2,) cosh Kl 


Vi= 


The problem of greatest interest in cable telegraphy is to deter- 
mine the current at the receiving end, usually designated arrival 
current. For this, we have the expression given by the second 
equation (76). To develop the complete solution from this 
expression by the application of the expansion theorem, it is 
necessary, in each case, to determine the roots of the determin- 
antal equation, 


R K 
Z(p) = (x ae a) sinh KI + (z,+2,)cosh Kl = 0, (77) 


or 


— (2 a= ey 


R K 
K = 220 p (78) 


tanh Kl = 


When the terminal apparatus z and z, are combinations of 
inductances and capacities, the determination of the roots of 
equation (77) may be a quite difficult matter. In special cases, 
however, the roots of the equation are readily determined, and 
the complete solution can be obtained. We shall consider here a 
few typical cases, 
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Resistance at Receiving End; Transmitting End Grounded.— 
The determinantal equation for this case simplifies to 


Ro Ro 
tanh Kl = ee I = — pl (79) 
Put 
Kl = jz, 
and the above transforms to 
tan x = —# x: (80) 


The roots of this equation can be found approximately by means 
of a diagram plotting two curves 
y; = tan 2z, 


Ri (81) 
lee Ae 1 ae 


the intersecting points of these curves giving the values of z, 
the roots of (80). In Fig. (20), tangent curves are plotted, 


Ai aes 
Vea ee eee Pee Soe F| 
FEC EERSEEEE EEE EEE 
cEEEEERSEP EES 


Fig. 20. 


yi = tan z, and tne two lines yz = —z and yz = —2/10, that 
is, for values of Ry = Rl and Ry = 4oRl. 


OCEAN CABLES 89 


The intersecting points for ye = —/10 are 
Pome L2.O% Onis Ost, 1162 14.8" 17583 
For 
Yo aad 
x = 2; 5; 8; 12. 


These values of x are the first few roots. By extending the 
graphs, additional roots of the equation can be obtained. 

To obtain the values of 0Z(p)/dp, we note that, for the condition 
of this problem, 


Z(p) = sinh Kl + R, cosh Kl, 


and 
OK 


ett Lis 
dp 


R 
ap Rl cosh Kl — R sinh Kl + Rol sinh Kl) 
DT LS Dea SS 
dp Ope 2% 
Hence, 
dZ(p) _ 1. 
dp =p 
Substituting Kl = jz, we have 


p ae) = (2 cos 7 —~ sin x — Rox sin (82) 


Rl nel Rie if Gn Ri BR ont kl | : 


By (80), however, 


Ro = —- se tan z, 
x 


and substituting this in (82), we get 


aZ(p) 1 Riese 
Dp = g( Bl cos 2 — x Dedtaaaa 
Rl 1 sin x , 
sag) & DS es ) se 


For p = 0, 


Kl 
Z(p) = BS + Ry = RL+ Ro. (84) 
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Substituting the values from (83) and (84) in the expansion 
formula, we obtain the complete solution 


area es : ate De ap é sin sy (85) 


the summation to extend for all the values of z, the roots of 
equation (80). When R, = 0, equation (80) gives tan 2 = 0 
or z = nz, and (85) reduces to 


= A =— a cos (nme REM, (86) 


which is the same as (27) for x = l. 

Condenser at Receiving End; z, = 1/cop. Transmitting 
End Grounded; z = 0.—By (76), the determinantal equation 
for this case is 


Doyen oe ie sinh KU + ta cae 


or 
>= Ke Cl 
tanh Kil = Oc & cake (87) 
Put Kl = 72 as in the previous case, and we have 
CLL : 
tan x = meee (88) 


The roots of this equation are obtained in the same way as in the 
previous case; plotting the two curves 


y1 = tan 2, 
; ee is 
Ya Co 


the intersecting points of these curves giving the values of zx 
required. The value of 0Z(p)/dp is readily obtained, thus: 


OZ(p) 1 1 Rl aes 
ae = Rl cosh Kl — : yj Sinh Kl +, (— sinh KI — —¢,,c08h Kl 
Pol elk? aaa : RCE 
= ofl cost — 5 sin x + 2 Cox sin 2 + Con? cos x. (89) 
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Substituting in (89) the values of 1/Cox = tan «/Cl by (87), and it 
simplifies to the following: 


VAG De! 
P Op cee 


COs & sine), Oy 


For p = 0, Z(p) = © no steady-state component, since there is a 
condenser in series with the cable. Substituting in the expansion 
formula, we obtain the solution for the current at the receiving 
end, 

2H < Rew" 

iL eaaphy Boiee 


i, = 


(91) 


COS & x 


the summation to extend for all values of x, the roots of the 
equation (88). 

Condenser at Transmitting End; z = 1/cyp. Receiving End 
Grounded; z, = 0.—The determinantal equation by (77) is 


Z(p) = dt x sinh Kl + @ os cosh Kl = 0 


or 


tanh KI = — aes =— oR (92) 
This is exactly the same form as (87), and, hence, the solution 
for this problem is the same as in the preceding case, except that 
1/co is replaced by 1/c. 
Condensers at Transmitting and Receiving Ends 7% = 1/cp; 
Z, = 1/c,p.—Substituting these values of z, and z, in the determi- 
nantal equation (77), we have 


| 


BPiay 21. 


R NG 1 
Z(p) = + am 2K) sinh K+ (or + ch, ) cosh KU = 0, (98) 
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or 
1 1 


he he Opa 
tanh KI = R Ne (94) 
Kk C,C yp? R 
Introducing the value of p = K?/RC, the above transforms to the 
following: 


RCP(C, + Ci) (95) 


2(72]4- 
RKPC,C, + athe 


Put K = jz, and we get 
RCP(C, + Ci) 


j tan «= Oop 
ra jac.c; +o) 
jx 


CG, + C; 
tan z = ONS — ay . i (96) 
a ZR ps PE ie Et 
PROM OF oe l Cop ; 


The roots of the equation can be obtained by graphical method 
if the ratios C,/Cl and C,/Cl are known; plotting the curves 


yi = tang, 
C, ase C. 
bel eect a (97) 
eC re 
"OP x 


the intersecting points giving the values of x which are the roots 
of equation (95). 
Differentiating (93) with respect to p, we get 


Ute) _ (_ Bek Loot OK eM pe Re 
= \ Kap Ct ee ope Cla ie be 


op 
Tee GENO ORD a Lf tent, a(KD 
+ ( +EG, »? * Jeosh Kl ok ae a. 7 sinh Kl 


i Op Op 
1/1 1 
= BG + a) cosh Kl, 
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Introducing the values p = K?/RC and 0K/dp = RC/2K, the 
above reduces to the following: 


2) = Ri|— bes i cewal vale ) ea 


x + 


op 2Kl 20,0, Keb + 2KING, CG: 
Wome lee 80) 9/11) 
THO +o CG, RR wale = >) cea 
Substituting Kl = jz, we get 
aZ(p) _ fe 1 3C# 1 ue oN 
Ua ra| Tra DON Ghee Shiv Wao anaes 
Vee Cl tee 
+|b-g00e tae, ta) | om] 8) 


Introducing these values in the expansion formula, we obtain 
the complete solution for the current at the receiving end. 


: E € 
=> Be orn +or(2 ft a 
2 ne Oy ae Ce 0; 
fA ae IPR Oc) EAN Teo 1H fe | 1 
De sal 2 Cres L 5 iG a =) 
the summation to extend for the values of x, the roots of equation 
(96). 
When the transmitting end is grounded, 1/C, = 0, equation 


(99) reduces to that of (91), which is readily seen on comparing 
(98) with that of (89). 


x2 
— Rot 


sin x 


ee cos x (99) 


CHAPTER V 
TRANSMISSION LINES 


Circuits of Distributed Inductances, Capacity, Resistance, and 
Leakage.—In the preceding chapter, the discussion was confined 
to problems relating to ideal cables, negligible inductance, and 
leakage. We shall take up here the more general problem in 
which all the electrical constants of cable or line, that is, R, L, C, 
g are taken into account. 

The relations between voltage and current at any point on the 
line are given by the following equations: 

(Lp + R)t = 7 
x 

«Sas (1) 
a gat) A a a 


From these we derive the equation of propagation 


PV 

(Lp + R)(Cp + g)V = 7a: (2) 
Put for brevity, 

Kes (lee Rg); (3) 
and use the following notation: 

a 

EG =— ye 

ihe neers (4) 

Ba oCe 


p=at+bjc=a-b. 
We have, then, 


1 
= ia(p + 2a)(p + 2b). (5) 
The solution of equation (1) is 
V = Ack* + Be-*=, (6) 


The expression for the current is readily obtained from (6) by 
the aid of either wae (1), thus: 


aps Ack + Be-**), (7) 
94 


i= 


inet 
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and with the notation given by (4), 


= ub p + 2b ay Wee —Ka 

ae Ser Ack* + Be-£*), (8) 
A and B are to be determined from the terminal conditions in 
any given problem. 


INFINITE LINE 


We shall consider first the case of an infinite line, and a constant 
voltage His applied at z = 0. A and B are readily determined 
as follows: 


Am—n (sande —wie 


Hence, 


ee p+ 2b x, 
al AP LPG (9) 


This is a symbolic solution involving the differential operator 
p from which the real solution is to be developed. For x = 0, 
the current at the transmitting end is given by 


weet p + 2b, 


ele ts p + 2a — 0) 


The operand, in this case, is unity-time function. Henceforth, 
the unity factor will be omitted from the equations, and if no 
other operand, unity-time function will be implied. 

To start with, we shall simplify the problem by neglecting the 
leakage, g = 0, b = 0, and equation (10) reduces to 


on dept eed Dee 
FE i p + 2a 
Now, 
% 4% 
P = coe? ) , (11) 
p + 2a p + 2a 


and e*t may be shifted to the right, provided we change simul- 
taneously p to p — a (see p. 9), which makes 


M % 
P ; = Ew Pie & ‘| “ea! 
p + 2a pta : 
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the operand changed from unity to e%. We have, however, 


apy and making this substitu- 


shown previously that «% = 5 


tion, we get 


ae SY ae ae p = Pe pa en: 
p + 2a pta/ p—a Tete ie 


Fo 
p? 


The operand is again unity. Expanding (12) by the binomial 
theorem gives 

Es, 1 = la 1-3 a, 1-3-5a° 

€ a ="€ Sea DEON | pit a | (13) 


ay 
operating on unity by each of the bracket terms involves suc- 


cessive integrations, that is, 1/p" = t”/n!. 
Hence, 


eu 1 = eat 4+. 


(=3) 
Pp 


The bracket term in the above expression is the expansion of the 
Bessel function of zero order, Jo(at), and, therefore, 


eeey = e—*] (at). (15) 


at? ate a’t® 
22 22.42 P. 42.62 7 a | 


(14) 


ae 
Oe 
Substituting the values from (15) in (11), we finally obtain the 
expression for the current at the transmitting end of the line, 


: E 
40 = Ipt 106): (16) 
A short table of the values of J,(at) for arguments from 0 to 6 
is given on page 164. 
For large values of the argument, the function approaches the 
value given by 
eat 


ll Ee 
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For negligible inductance, a is infinitely large, the above relation 
holds true, and, on substitution in (16), we get 


ee BV Le 
Lv 
Viomime Lalor XK, 


i Eulg Wa (17) 


On comparing this with equation (7), Chap. IV, the expression 
for the transmission current in a non-inductive cable, we obtain 
the relation 


hy = 


es 1 
Vl Ee A at 


which yields an additional proof of the validity of the operation 
of the fractional derivative, which was made use of extensively 
in the preceding chapter. 

In passing, it may be observed See we have incidentally 
established the equivalence of the following operations, all 
leading on development to a series which is the zero Bessel 
function, namely: 


Pee? ip us c = 0) bee: NE eeoe 
p + 2a pta (pa) 


1 


(3 
P 


Leakage Included, Applied Voltage Varying as e¢«°’.—In 
the more general case including leakage, we have, by (10). 


= Ii(at).. (13) 


ne ee p+ 2b 

home le p + 2a 
(Dae ee 19 
pip +o oy 


This may be transformed by multiplying by e~*‘e*, shifting e* to 
the right, and changing p to p — p, which would give 


iam B 4 [Pe (20) 
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If we should take the applied voltage to vary as e~**, the operand 
would be changed to unity and the above reduce to 


ye EL “of? a a 
ha: p +a 


7 o»\-% 
“fede Me-3)" 
Lv p Dp 


By (18), 


hence, 
ince Fer = * ie (21) 


To obtain the complete solution requires operating by 1/p on 
I,(ct), that is, integrating from 0 to ¢. If we use the expanded 
form of Io(ct) given by (14), the terms of which are readily 
integrated, we obtain the following: 


343 1 - 3e@5t? 
inlet): Pe ey te ee ae: |}. (22) 


eae 


For t = 0, the initial current is E/Lv, and the final current is 
zero, obviously, because the voltage falls to zero. 

All Electrical Constants, L, R, C, g Active; Steady Voltage.— 
To derive the expression for the initial current in a line, all 
electrical constants active, under the application of a steady 
voltage, we proceed in the following manner: 


Ae Ee + a 

een F p + 2a 
_E__(p+%) 
Lv (p + 2b)4(p + 2a)” 


g 
_z_ +8) 
~ Lo {+p = oe 


(23) 


TRANSMISSION LINES 99 


The development of p/{(p + p)? — o?}” is readily obtained by 
transforming it’to one of the known forms given by (18), thus: 


Aves See = erte—ot Pee ca Cue (ae pret 
\(p + p)? — 0} {(p + p)? — o}% {m2 — o2} 4 
pan th ee ee Revie 


C2 EA (ae 
—pt 
- — gw = eI (ot). (24) 
Co 72 
t-%) 


Introducing this in (23), we obtain 


iy = = (1 * fate (25) 
The initial current due to a steady voltage is thus expressed in 
terms of a known function and its time integral, meaning, of 
course, integration from 0 to ¢. Heaviside indicated several 
methods for effecting this integration, all leading, of course, to 
the same result. Shifting «~** to the left by changing p to 
p — p, we may write: 


of —pt 
ie Toot) = aoe Thea eter (9G) 


Pp —p p 


cae | 
p 
Expand - » and use the expanded form of J(ct); we get 
—1 
i p 
2 3 
Pear lict) = cE ete yee ae ) 
| pee De aD 
p 
ot? ot! o t6 
+3 + ap to. gg t ) 


This is integrable by sight, giving the following result: 


1 242 4y4 
=I (at) = af a = fs a eae ) 


a 
p 
ieee ott! 
242 = a ane 
$0) + oa t ais t ) 
il ot? Gale oe 
tesa tee aatee edt ) 


+ . : 
= u(t), say. (27) 
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The law of formation is obvious, every bracket term being the 
integral of the preceding bracket term. Hence, 


if —PtT (ot) = “ul. (28) 


Combining this with (25), we obtain the completely developed 
expression for the initial current 


eet s 
ee i 


which may be put in this form: 


Oy a (ro Iy(ot) + (1 Fa “ul : (29) 


For g = 0, ¢ = p = a, and the above reduces to (16). 

The Current at Any Point Distant x from the Transmitting 
End.—To obtain the expression for the current for any point on 
the line, it is necessary to operate on either of the expressions 


obtained above for the transmitting end current by the operator 
—Kz 


PI g(ot) + et 
Cp 


€ 
NON-DISSIPATING LINE 


We shall first take the case of a non-dissipating line free from 
energy losses; R = 0 and g = 0. The transmitting end current 
is given by 

w= jo 
If # is constant, the current jumps to its full value E/Lv and 
travels on with the velocity v reaching any point distance x at 


time t = a/v. The current jumps to its steady value at any 
point distance x at time ¢ = x/v. For EF = f(t); % = f(t), and 


ig = Ee-®*f(t). (30) 


Expanding e~**, remembering that for R = g = 0;K = pVLC = 
p/v. 


0 Lied ive a 
io (i-2 +p Pet). (31) 
If we compare this with Taylor’s series 


fet n= (1th + oe es Heer fie 
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we see that they are equivalent, and we may, therefore, write 


tz = s( - i) (32) 


Whatever the character of the current at the origin, it travels 
at a speed v without any change whatever. This is also true of 
the voltage on the line; if the voltage at the transmitting end is 
Ef(t), the voltage at x is Ef(t — x/v). 


DISTORTIONLESS LINE 


When the electrical constants of the line are all active but 
related to satisfy the condition 


the line is said to be distortionless; that is, any current or voltage 
wave impressed on the line is propagated at the speed v without 
any deformation; the waves suffer uniform attenuation in their 
transit on the line but no change in form. For this condition, 
' a = b, and by (9) and (5), 


te = Te 
and 
1 
K = me + 2a). 
If 
Ey = Ef(t), 
Hae Be >” f(t). 


Similarly, 
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Since the time of transit for the waves to reach any point x 
is t = x/v, we may write the above equations in this form: 


E % 

a crey(1 = ) 

ly | v (33) 
Bhee"y ( = a 


The current and voltage are in the same phase and in constant 
ratio for all points on the line. 

This may, perhaps, seem a little clearer by considering the 
steady-state condition for applied alternating voltages. If the 
frequency of the impressed voltage is w/2zr, the steady-state 
solution for the voltage and current on the line is readily obtained 
by replacing p by jw in equations (5), (6),and (9). We get, then, 


te 


i, 


i , 
KO = (jo + 2a)(jo + 26), 
V = He *7i, 
i= Lan r 2» ee. 
LuNjo + a° 
For the distortionless condition, a = b, the above reduce to 


K = "(jo + 2a), 
i Eat ée(7), 
pies a eae Gea) 


The real parts of the above equations give 
_ 24, x 
V = He >” cos of — a 
= — ate 
v 


The waves travel along the line at the speed v and are attenu- 


(34) 


ated at the rate «~*’", Since ¢ — : is the same for all points 


on the line, there is no change in phase in the waves as they travel 
on the line. The waves retain the same form though they are 
diminished in size as they travel along the line. The attenuation 
of the waves, in this case, is independent of frequency. Hence, 
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any complex wave will travel on the line without any deformation; 
the waves arriving at the receiving end may be reduced in size, 
but otherwise a fair copy of that sent out at the sending end. 

Inductance, Capacity, and Resistance Active.—The current at 
the transmitting end is, by (16), 


4 E 

iy = Tyo Lola); (16) 
and 

tn ered (al). (35) 


For this case, 
K = (Vp + 2ap = 1 V(p + a)? — @&. 


By shifting e«~* to the left, and compensating for it by changing 
p to p — a, a process explained before, equation (35) is trans- 
formed to this: 


; E Sy rar 
n= nee ae I,(at). (36) 


Expanding the exponential factor, and rearranging, the equation 
takes this form 


; lik er ipa). sp — ae) 4 

Cras iE i ve 2! hs vf A! ote (37) 
x? (p? — a?) xt (p? — a?) iz r: 

= E ak a He x B i Teo Be — a?) Hota 


We have, however, by (18), 


Ge ae! = (at), 


and — 
(p? — a?)*I,(at) = pl = 0, 


that is, operating by the second bracket term on J,(at) produces 
zero result and can be thereafter neglected. It remains therefore 
only to evaluate the first bracket term in (37) in its operation on 
T,(at), that is, 


: IH AGUS UE LN pele a oe 
z= 76 Z Peg aR ary aime? 0 a I,(at). (38) 


104 HEAVISIDE’S ELECTRICAL CIRCUIT THEORY 


The Bessel function satisfies the differential equation 


7 Sat + — ; tT To(at) — aI,(at) = 0, 


and (39) 


(p? — @)Io(at) = —B (at) = —aite), 


More generally, the function satisfies the relation 


T,,(at) _ (Qm +1) Tn(at) 
Wes atm a t Ptat)™ 


einer Dane (40%) 


Hence, the following: 


(p? — @)I(at) = —P (at) = — ga | 

(pt = at)PI(at) = —(@" = ata = gat a 
2 2)3 = Oye ey ie pe etey a gf (at) 

(p? — a?)®Io(at) = 3(p* — a?)a laiik 3° 5a Cat)? 
Si RO ed oe 8 gla(at) _ ce hate 

(p? — a®)4I(at) = — 3-5(p? — @)a (at)é =3-°5-7a (at)! 


Introducing these values in (88), we obtain 


Ge Uigah = ex) Taal) 


oe iat | Poa = 


v2 at bt 2-4 (at? 
PT dines 1 ‘I3(at) 
ye 2-4-6 (at)? re : (42) 


The bracket term of above equation is the equivalent of 


Te Swe — 37)* This identity can readily be established 
( 


by putting both expressions in series form, arranging and com- 
paring terms in the series of powers of 22. We thus arrive at the 
following expression for the current at any point on the line: 


. E 
te = —e “ly 


Lv 


1 See note on Bessel functions, in appendix. 


(ve? = ays], (43) 
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Comparing with (16), we see that vt is turned into (vt? — 2?) 
in the J) function. 

At the time ¢, the region occupied by the current extends 
only to the distance vt from origin. Beyond this distance, the 
current is zero. Equation (41), therefore, expresses a wave whose 
front travels at speed v. At the wave front, the current is 
E 
Lo 
the current curve along z, at successive moments of time, can be 
readily calculated. 

For a non-inductive cable, L = 0, the argument of the function 
in (43) approaches infinity as L approaches zero value, and the 
value of the J, function is given by 


e—, Since the function has been tabulated, the shape of 


14 
© (p22 — 22)? 
€? 


Io 
" am (ve — 7°) 


a (Can Sate a) % 
v 


Introducing this value in (43) gives 


ax? RCxz?2 
: Hoe 203 Hemet 
Vx = = _- = ee 
Lv V/ Qrat at 
G 


the same expression that was obtained for this case by direct 
development (see formula (15), Chap. IV). 


DERIVATION OF VOLTAGE WAVE 


The formula for the voltage can be obtained from the current 
formula by the relation between current and voltage, which is 

dV ; 

dz a (R ai Lp)i, 

= L(p + 2a)i, 


=~ + daeetle| H(t — 21)4|. (44) 
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Shifting «~* to the left, changing p to p — a, we have 


ae dV = Fe + a)Io 


7 ayy os AVL 


The voltage is the negative of the z integral of the right-hand 
member of the above equation. The integration can be readily 
effected, using the bracket term of (42), the equivalent of 


Ty pre — 2°)”!, which gives the following: 
ee p\ | ax _ faz\* Ii(at) 
y= —Beo(s +2)|81an — (22) Ee, + 


az\*® I2(at) ee 
Greer | x46) 


We must add to this a constant of integration to satisfy the 
condition V = H when z = 0. Since the above makes V vanish 
at the origin, the constant of integration is obviously #. Hence, 


at P\| (ax ax\* I,(at) 

iia (1 Bs a CG ota = CS) Sa 
ae\® — Fs(at} a: 
(leeetc: 


The complete development of the solution for V involves one time 
differentiation. We have, by (40), 


Tn(at) _ o,m+(at) 


We IBGE 


(47) 


(aty™ — " (at)mrt 


that is, 
pI,(at) = al, (at), 
Ii(at) _ 2h), 
at at 
I,(at) om I3(at) 
Paty? ~ © (at)? 


Introducing these values in (47), we obtain the complete 
developed solution for V as follows: 


— —at Ci D et ax ox pibras2 © =— 
Ute he ls TOOLS ler ican 


ax\?* I,+T; 
(=) G-actin 


The argument of all the J functions is af. 
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Another way of deriving the expression for the voltage is to 
work direct without any regard to the developed current formula 
from which (48) was derived. The voltage on the line is given by 


H, = He£2 
K = VLCp? + RCp = Vp + 2ap, 


1 
- 5V 10) = ar. 
We may put 
Hepat Di cae erg" 1) a0 


Shifting «” to the right by changing p to p —a, the above 
transforms to 


af Wana ae. 
Bes ea” eu. (49) 
: : — hipaa : : 
Expanding the exponential factor e in a series, we 
obtain 


B. = Be*}1 — 2p? — a) + Sp - 2) - 
; v 2!v? 


3 [pa (P” Set a?) + Selene |e, (50) 
operating by (p? — a?) or an integral power of (p? — a?) on 


e*t is obviously zero; hence, the integral powers of (p? — a?) in 
the above expression are to be disregarded, and it reduces to 


2 
E, = Bes] = "1 ae <0) 


ee Vor — a’) le (51) 


By (16), 
Ne = Gove 
(2 
and 
(p? — a’)4e" = (p + a)Io(at). 
Therefore, 


gan newer —2(14 Xe ee 


) 
ae JO Se aoe rota} (52) 
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The results produced by operating by (p? — a?) or integral powers 
of (p? — a?) on Iy (at) are given by (41), and, substituting these 
values in (52), we obtain 


E, = Ee* le = aa + PV(r(ap = aaa pestle 


at 
x‘yt — Is(at) _ 
2:4-5v* (at)? es )| Sa) 
which is the same as (46), and developing it fully will give again, 
of course, (48). 

When all the electrical constants, including leakage, are taken 
into consideration, and steady voltage applied at the origin, the 
problem is much more difficult. The final formulas assume a 
form which is too complex for interpretation, either physically or 
numerically, and the discussion of this problem is, therefore, 
omitted here. Heaviside gives the complete derivation begin- 
ning on page 312, vol. 2, Electromagnetic Theory. 

Applied Voltage Varying as «~?’.—A special case, however, of 

g 
some interest is when the applied voltage varies as care ean 
In this case, 


Eek? |p + p+ 2b. aete 
Lv pt 2a° 
omer FES + 2b _ ons 
pt a 
Shift «*% to the right and «~ to the left, which changes p to 
p — a — b, and the above transforms to 


le = 


=e ema etna 


. _ Be-et -2Ve@ta@-+) |p —o,, 
ee 9 P= Len, (54) 
By (16), we have 
| g < ert — 
t + s) Fee 
Hence, 
: Ee— et —=(p2—o2) 
ze = ae u I (ct). (55) 


This is exactly the same form as (36), for which the developed 
solution was obtained, and can be applied to this problem, which 
gives 


3 eet 
iz = “Ty 


© (v4? = x?) 74 


(56) 
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FINITE LINES 


The method of treatment adopted in Chap. IV in connection 
with the solution of problems relating to cables of finite length 
will be followed here, applying the expansion theorem, which is 
simpler than the direct operational process. We shall first con- 
sider the case of free lines, that is, either grounded or open at the 
ends, no terminal impedances, in which complete reflections at 
the ends with or without reversals are produced. The investiga- 
tion will be carried through for two cases: steady and alternating 
voltage applied at the origin. 

Line Open at Far End, and Steady Voltage Applied at Origin.— 
The voltage and current at any point on the line are given by (6) 
and (8), 


Ack* + Be-*e, 


ee pete 2b ce neh 
er (2) { —Aek* + Be-k*} (57) 


K = *(p + 2a)4(p + 20), 


The terminal conditions for this problem are 


x=0;V =#, 

z=l;7 =0. 
Hence, 

A+B=EH 


—Ack + Be*' = 0. 


From which the values of A and B are determined as follows: 


Ee-*®! 
Bee er 
one te (58) 
B= ex eK! 


Introducing these values in (57) gives 


V z E{ ke) -- eXt-®)} 


kl 4 @-Ki 
% (l—z) —. ¢-K(Il—z) 
we E(p + 26 & fi : (59) 
Lv\p + 2a Xl 1 ¢—Kl 
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Expressed in hyperbolic functions, 


1 cosh KD a2) 


y cosh Kl 
pee vi 
nl ae + a sinh K(l — 2) (60) 
Iv\p + 2a cosh Kl 


The expansion theorem can be applied to either V or 7; but it is 
sufficient to develop the solution for one, the other is obtained 
from the circuital relation between V and 7. We shall develop 
the solution for V. 

The determinantal equation is 


Z(p) = cosh Kl = 0. (61) 
Hence, 
Kl =j53n=1,8,5-°°°, 
and 
nny? 
KY = Ap 20)(p A- 2b ae 


From this, the values of p are readily determined. 


n? 242 
pa = —p £ alo? - a 


SS Palms (62) 


Also, 
OZ(p) _ an KPA 
Op j 
sinh Kl E 
ey jee as (63) 


Introducing the value of Kl = jnr/2, and the value of p as given 
by (62), the above transforms to 


sin — . 78 is sin 
OZ ; m2 2 
OO) ioe ees (64) 
Op nr v nT 
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For p = 0, 
ois 
cosh V ab, 


y cosh sl 
pers ae 2) aa 
v 


cosh 


where, for brevity, we put s = 2\/ab/v. 
Substituting in the expansion-theorem formula, we obtain 
the developed expression for the voltage, 


— pte +IBnt re 
y = Heosh lay 4 a e~ Me #16! cos SF (1 — 2) 
a cosh sl ee a, 
+(=pt i 
a ( Pes IBn) ohm sil ) 
Eat ne a oe 
. Ecosh s(l — 2) LR oS a nav SIN oy sae 
* cosh sl : fe 296 ate Jen) 


The complete solution must include all terms which come under 
the double sign in the summation term. Hence, 


pee cosh s— 2) | Mare _ -. We 
V= : a >on sin “97 


cosh sl or © ey, 
| cos Bat +7 sin Bat cos Bat — 7 sin Bat ‘ 
JEX—0 = JEs) jBn(o + jBn) 


which simplifies to 


, _ Heosh si — x) _ par? _,, 
Lee cosh sl Be > 


n sin Sa sin B,t + Bn cos Bnt) 
Br(p” =f Bn) 


(67) 


This is the complete solution for the voltage, the first right- 
hand term giving the steady-state condition, which theoretically 
obtains only after an infinite time, and the summation term giving 
the transient voltage component, the summation to be extended 
for all odd integrals from 1 to ». For x = 0, the summation 
term is zero, and the first right-hand term reduces to E, as it 
should. 
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The expression for the current can be derived from (67) by 
the circuital relation 


di 
= = (Cp+q)V. 

This gives 

di _,,coshs(l—z) _ ,rv’ _ MTX 
dx Hg cosh sl ey [2 > arias 

jenn Bet 0s os ES Cea 
Bn (p? ae Pz) Bn 

and : 


_ Egsinh s(l — z) Aree nmr 
“Ug cosh sl are Day ee 


ge SIN But + gBn COS Bnt _ sin Brit 
Bulp + Be) Bn 


» (68) 


which is the complete developed expression for the current; 
the first right-hand term is the steady-state component, and the 
summation term, the transient component. 

For g = 0, negligible leakage, equation (66) reduces to (p = a, 
in this case) 


ae ae But 


=== EO bay Sa (69) 


Line Grounded at the Far End.—The terminal conditions for 
this case are 


=0;V=H, 
keV =O 
which give, on substitution in (57), 
A+B=E, 


Ac + Bek! = Q, 
From these, A and B are determined, 


Pia aly 


B= 
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Substituting these values in (57), we obtain the following expres- 
sions for the véltage and current: 


¢K(—2) A=), | 


—e€ 
7 d 
eh — «= 


fe He ae a 
~ Lv\p + 2a 
Expressed in hyperbolic functions, 
Ve E sinh K(I — x) 
ip sinh Kl 


oe Ae + ar cosh K(J — v) 
os p + 2a sinh K] 


Kl 


es z| ae 


eK(l—2) 4 ¢-K(I-2) | 


eK! pees eK! 


(71) 


Applying the expansion theorem for the development of the 
solution from the above expressions, as in the preceding case, 
we have here the determinant 


Zip) = sinh Kl = 0, (72) 
Hence, 
Kl = jnr, 
and 


K 


SV @ + 2a)(p F 2) = 


from which p is determined as follows: 


De = pct Ee 2 ; 


12 cosh K! (74) 


sinh sl 


BODY = 2~/ ab au ~ sinh s(L — x) (75) 
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Substituting these values from (73), (74), and (75), in the 
expansion formula, we obtain the developed solution for the 
voltage, as follows: 


E sinh s(l j sin nal ~ eae 
ae Sal 24H Dae fp i 
+78n(—p as qBn) cos ee 


pst ea 


E sinh s(l — z) v? i 

= — Be re- : 6 

sinh sl Ep ra > + jBn(—p + jBn) ie 

The complete solution to include both terms which come under 
the double sign in the summation term is 


_ Esinhs@— 2) Env’ _y : ne eos Brt +7 sin Bat 
— ont Sr sin j8n(—p + JBn) 


cos B,t — j sin Bat 


sinh sl l2 


+ 58.0 + iBa) 
which simplifies to 
_°E sinh st — 2) 2H vr aby 
a sinh sl ea re 
n sin “*{p-sin Brat + Ba cos Bat} 
2 2 : 7 
= Bale? Ba) Se. 


This is the complete solution, the first term giving the steady- 
state voltage, and the summation term, the transient state, the 
summation to extend for all integral values of n from 0 to ~. 
When x = 0, the above reduces to V = E, the applied voltage 
at the origin. 

The expression for the current is obtained from the voltage 
formula by the circuital relation 


di 
me a (g + Cp)V, 


which gives 


di _ sinh s(l— x) 2H t - Ta 
eM ae op en 
| ge SiN Brt + gBn COS Bat 


L Bn (0? = Bn) Bn 
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and 


cosh s(T— x) _ ae, cos em 


sinh sl l l 
Jp SiN Brt + gBn COS Bat SIN Bat 
—C —t- =(78 
For g = 0, (78) and (79) reduce to the following: 
fe 2Hnv? 
V = ae ae s) = 7 €nt 
n sin ae sin Bnt + B, cos Bnt} 
- (79 
S Bala? + BE) 2 
lei Nx SiN But 
ade Ce 2S cos — Bn (80) 


If the line is infinitely long, either (69) or (80) should give the 
same result as (43), which was obtained by direct operational 
process for the same condition. The equivalence of the formulas 
(43) and (80) for! = © is not apparent, but it can be shown that 
(80) transforms into (43), which we shall now do. The sum- 
mation term goes by steps of 7/l, which become infinitely small 
when / is made infinitely great. We may put, therefore, 


Eom ene 
Le: m4 = dn: Bm = oe —@= V/ my? — a’, (81) 


and the summation term in (80) is converted into a definite 
integral, that is, 


je PHC | cos ma = ca (82) 
1 0 Bn 
Evaluation of the integral, 
= il cos mz Bml am, (83) 
0 Bm 
Ca = V/ mv? — a’. 


B is a function of a”; differentiating sin 6t/G with respect to t and 
a’, we obtain 
d sinpi_ d 


A P sin Bt 
dtda £6 da 


2 8B 


(84) 


5 cos Bt = 
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that is, 


d sin pt _ 1 (“tsin Bt), _ t sin Bt. 


5 
da 8B 2Jo 8B 2p p oe 


By Taylor’s theorem, 
4 
f(a’) = f(0) + af) + 57) + + - 


Considering sin 6t/8 as a function of a?, and remembering that for 
a=0, B =m, we get, by Taylor’s theorem, replacing the 
differentiations by successive operations of t/2p in accordance 
with (85), the following: 


sin Bt _ sin mvt , a’tsin mvt a@t\’ sin mvt 
B mv 7 2p mv zi (=) Dine 2 88) 


Substituting this in (83), we get 


ee) 

2p 2 sin mvt 

a nr, i cos oie (87) 
4 mv 


at 
oma ET i hn af 


The integral may be put in this form: 


{I cos mx sin muta al {sin m(vt + x)+ sin m(vt— X)\ am 
0 0 


mv os), mov 


We know, however, that 


“sin mz T 
|, idm = £5 (88) 


having positive value for x positive and a negative value for zx 
negative. Hence, 


“cos max sin mut rg 


T 


2 


when zx < vt, 
= ae — 1) =O whenz > vt, (89) 


that is, the operand which is the integral in (87) is zero when 


pees 
v 
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If we use the operand 1 only in (87), the limits of integration 
must be from 2/v to t instead of 0 to t. We finally arrive at the 
following: 


oe ae a Viral cc a De | 
w= Fit ye [et a | tae | edt + - - a 


8 


aes Lee 1 4a)? ey 
see ats oe ke 
2a) (ie ey, 
Trae ec ry we) 
a?\/# 
If we put al 2? — 4) = y, the above takes the form 
a a | fea y° 
Cap ee ee reat Drange 
ae spell | 2 _ vw a L 
= ay oly) = ayo a(t — >) (91) 


Substituting the value of wu from (91), which is the value of the 
definite integral, in (82), we arrive at the expression for the 
current, 


Tv 
He : x? “a 
= ay rae | 1G — ) ; (92) 


which is the same as (43), and thus establish the equivalence of 
the result obtained by two different processes—the direct 
operational process and the application of the expansion theorem. 

Terminal Impedances.—If the effects of the terminal apparatus 
at either end of the line are to be included in the general solution, 
the problem is by far more difficult. It is only in special cases 
that it is at all possible to obtain the completely developed solu- 
tion by the application of the expansion theorem. We shall 
give here only an outline of the problem and indicate some special 
cases where a solution is possible. 

Consider the most general case: a line of length /, impedances 
Z, and Z, connected at the transmitting and receiving ends, 
respectively. The solutions of the differential equations given 
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by (6) and (7) are perfectly general and also applicable for this 
case. We have 


V = Ack + Be-**, 
— Kk poss Ker —Kz (93) 
1.= Lp +R! Ae 4- Be K a 
where 
K? = (Lp + R)(Cp + 9g). (94) 


Ere. 22.- 


The constants A and B are to be determined from the terminal 
conditions, which, in this case, are 
when 


£=0; Vo = E — Zio. 


Dect Vee (95) 


Introducing the values of V.z.o and 7,0 from (93) into (95), we 
get the following two equations: 


ZK 
Ack! + Be-Ki = In pl Act + Be-*'), 


Rearranging, we get 
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Solving for A and B gives the following: 


A = ——- — = 
VE ; 
(: is gue eK) 4(Z, 42) KOBE 
2K (97) 
ad a 
Z,L, ar K eRe 
(+ apt Ep) + Et pp e gle ben) 


Introducing these values. into (93), we obtain the expressions 
for V and 7 as follows: 


( ZK ae xr) ( ZK Je (l—z) 
. E\\ 1 ip + R € — LotR cae 


(+; ZOiK Goes exk) (7, +2)F pp te*) 


(Lp+R)? 


api Z,K K(l—2) ( — Z,K ae 
egealaecso) As ese 


( LiF? JO + Cet LT 


(98) 


(eX! ¢-K?) 


(Lp+R)? 


Expressed in hyperbolic functions as follows: 


Z,K 
E Dp ee Ko-»)| 


ZL? 


_ EK 
Sp Hh 
Ay het ig) Gog : 

(: + UPR) sinh Kl AF 


To develop the complete solution from the above equations will 
necessitate obtaining the roots of the determinantal equation 


Z SAK ° 
(1 _ ce i) sinh Kl] + 


Be (701 7.) cosh Kl) = 0, (100) 


Lp+R 


d 


ae 


jaar 7. + Z,) cosh Kl 


ne Kl — a) + lp sinh K(J — x) 


(99) 


ia es + Z,) cosh Kl 


Lp +R 
Z(p) = 2a" 


Gtk 
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which for most combinations of Z, and Z, is nearly impossible 
to do. In some special cases, however, the roots of this equation 
can be determined, which makes it then possible to apply the 
expansion formula and develop the solution. 

The determinantal equation may be put in a more convenient 
form. We have, for x = 0, 


fos Bie | 
pER hhc Kl+ Tp ack sinh cane 


7m) sinh KI + (Z, + Z,) cosh Kl 


E 


ho 


ane ae 
(1+ a Lp +R 


ante (101) 


Lp +R 


(102) 


The above may be put in this form: 


PD ten KEE, 
Z(p) = Zi + KZ. =0. (103) 
L + jai: tanh Kl 
Consider now a special case: Z, = 0, and Z, = nl(Lp + R); 
n is a numerical factor, the receiving apparatus consisting of a 
coil whose time constant is L/R. For this condition, (103) 
reduces to the following: 


(R + i} tanh Al + nt)) 
AYO) 1 + nKl tanh Kl 


tanh Kl 
. (R + rp) Kl n) 


1+ nKl tanh Kl 


= 0, 


= 0. (104) 
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The roots of the equation Z(p) = 0, (104), are given by 
R+ Lp = 0,| 


tanh Kl + nKl = 0. mo) 
For the first equation (105), we have the solitary root 
R 
Dom ay: (106) 


The roots of the second equation (105) can be obtained 


graphically. Put 
Kl = 42. 


N= ng, 
tanh KI = tanh jx = —j tan g, 
and the second equation (105) transforms to 
nz = tan 2. (107) 


The intersecting points of the curves y; = nx and y2 = tan x 
give the values of x which are the roots of (107). From the 
relation Kl = jx, the values of p corresponding to root values of 
x are readily obtained. 

To complete the solution, it remains yet to obtain expressions 


for the oe terms of the expansion formula. 


For the root p = —R/L, we note that for this value of p, 
K =0. Hence, 
tanh Ki = 0, 
and 
tanh Kl _ 
kl : 


which gives 
Z(p) = (R + Lp)l(1 + n). 


Therefore, 
pew) Fri iid pres ial kc (108) 
For the other roots, we have 
OZ(p) _ aie ME 
Sts Ip)KE, : ik 
Se Rae ai {1 + n cosh Kl} A (109) 


This gives all the data necessary from which the complete solu- 
tion can be realized by the application of the expansion formula. 
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It must be remembered, however, that to get the solution of 
either the current or voltage at any point distance x from the 
transmitter, each of the summation terms in the expansion 
formula are to be multiplied by the numerators of either equation 
(101), using the proper root values of p. 

Another case which permits of a similar reduction is the 
following: 


Z, = n(R + Lp) + ng + Cp), (110) 
Z, = n3(R + Lp) + nag + Cp)-?. 
M1, Ne, M3, n4 being numerics. That is, apparatus at either end 
consisting of a coil and a condenser in series, the time constant 
of the coil being L/R and that of the condenser C/g. 
By (103), we may put the determinantal equation in this form: 


t+ (OR tanh Ki+ = 7“? tanh KU Z, 
Z(p) = e TK =0. (111) 
ih 45 Ip+R tanh Kl 
This gives the relation 
tanh Kl —(Zp +2) (112) 


~ KL Z,Z.4g + Cp)l + (R + Lp) 
Substituting the values of Z, and Z, from (111), we get 
tanh KI _ —(m + 3)(R + Lp) — (m2 + mag + Cpr) 


(g+Cp)l nyn3(R + Lp)? + (ning + nens) 


R+ 2) None | 
(fe + Cmmery + (R + Lp)l. 


af —(m1 + 3) K? — (nz +14) 
Ui ningK4 + (niN4 + Ngn3)K? + ngn4} + IK?’ 
ae { (ni + n3)K 3 + (m2 + ns)K} F 
ninsK* + (ning + nons + 1)K? + non, 
The roots of this equation are readily determined by graphic 
method, plotting the curves 


(113) 


tanh Kl = 


(114) 


y: = tanh Kl, 

y2 = f(Kl). 
f(KI) is the right-hand side of equation (114), the intersecting 
points of the two sets of curves giving the roots of the equation. 
Once it is possible to get the values of the roots of the determi- 
-nantal equation, the solution can be completely developed by 
the application of the expansion formula. 


(115) 


CHAPTER VI 
ARTIFICIAL LINES 


Artificial lines are of practical importance in their use for 
balancing telephone lines or ocean telegraph cables for duplexing 
or other purposes. They also find application in the study of 
electrical phenomena connected with transmission problems on 
long lines. An artificial line closely simulating the characteristics 
of an actual line permits experimental laboratory study of 
transmission-line phenomena which would be otherwise imprac- 
tical. It is, therefore, of some importance to investigate the 
problem of wave propagation on artificial lines and particularly 
in the matter of transients. It will also serve our purpose in 
further illustrating the utility of the operational method and the 
expansion theorem for the solution of a difficult and important 
problem. 

An artificial line is essentially a circuit structure consisting of a 
large number of recurrent sections of series and shunt impedances, 
the series impedance of each section representing the resistance or 
resistance and inductive reactance of unit length of line, while the 
shunt impedance represents the capacity reactance, or the 
capacity and leakage combined, of unit length of line. In 
structure, it is similar to the filter circuit discussed in Chap. III, 
and the basic formulas (4) and (11) derived there apply also to 
the artificial-line problem. The fundamental difference is in 
the number of sections. For an artificial line to simulate closely 
an actual line, theoretically the number of sections should be 
infinite in number, or, at least, a very large number. In the 
filter circuit, only a few sections are sufficient for many purposes. 

If the total number of sections is n, the last section being the 
nth section and the first section designated as the zero section, 
then for any intermediate section, say the qth section, we have, 
by (11) and (4), Chap. III, 


ee E cosh (n — qv (1) 
¢ 22 sinh y sinh ny, 
123 


124 HEAVISIDE’S ELECTRICAL CIRCUIT THEORY 


and 
cosh y = 1 t 9 a (2) 


These are the general formulas ee for any type of arti- 
ficial line. What we are generally concerned with most is to 
obtain an expression for either the input current or the output 
current, that is, 7) and 7,. 

Artificial Ocean Cable.—Consider the case of an artificial 
ocean telegraph cable, negligible inductance and leakage, L = 0, 
g = 0. In this case, 


a= i 
a Cp 


Fig. 23 


If we assume an infinite number of sections, n = «, then cosh ny 
= sinh ny, and (1) reduces to, for the input current, 


E 


fo = Z, sinh y (4) 
By (2), 
sinh y = +/cosh? y — 1 = Ne +34) | 
2 22 
= fa, 14 
V2 a 42 
and 
ge sinh y = V/2iz2 + 142°. (5) 
Introducing the values of z; and z2, we get 
Pri Ve a iB (6) 
Hence, 
E E 
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The complete developed solution is readily derived from (7) 
by direct operational process. Put, for brevity, h = 4/Rc and 
expand the denominator in a series of inverse powers of p, thus: 


ib fp 1-3 h? 
; =pespEr ts 
ve + Be an +h 
iho Bhs By 1 
22 A Oy" 


Remembering that operating on unity function gives 1/p, = 
t”/n\, we have, by substitution, 


242 
— Sh +5 Coe (hg i 


ee 1: 
Zak Dy 22 
ee 
p 


which may be put in this form: 
ioe ht 1-3/ntY AG 

yrake iG) + an) SOHC aes menos 
Le 


The series of (8) is the expanded form of the function e~*/? ‘To(ht/2), 
which can be readily verified by using the series for the two 
factors, multiplying and combining terms of the same powers 
of ht/2, and comparing with the corresponding terms in (8). 


Hence, 
h 
See eka anh) (9) 
h 2 


Introducing this value in (7) gives the complete solution for the 
current in the first section, 


hy (ht 
sr ot ne ) 
moot (et 
ieee or ee) oy 


If the number of sections is large, the resistance and capacity 
per section small, the product RC is a very small quantity. For 
values of ¢, therefore, differing appreciably from zero, the argu- 


ee 


73° 5 Wet eed a a 
4-6 3! 
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ment 2¢/RC of the Bessel function in (10) is very large, and, toa 


high degree of approximation,! 
2t 


I eee 7 (11) 

Nec) a 
RC 

Substituting the value from (11) in (10), it reduces to the 


follo wing : 
ae oy ae Re «/ xt ( ) 


This is precisely the same as formula (9), Chap. IV, derived 
for the input current in a cable of uniform distributed resistance 
and capacity. This condition of the identity of the character 
of the input currents in a uniform cable and an artificial cable 
is realized only theoretically when the number of sections in the 
artificial cable is infinitely large. Practically, this condition is 
approximated with a reasonably large number of sections. 

To obtain the solution for the output current, that is, the 
current in the nth section, it is best to apply the expansion 
theorem. 

For q = n, (1) reduces to 


. E 
Ui 2, sinh y sinh ny. (13) 


The determinantal equation is 


Z(p) = 2 sinh sinh ny = 0, (14) 
which gives 


«ST 


ny = jer, and y = J (sil 9 2250s (15) 


By 2, 
Sir 1 
cosh y = cos Ae 1 +- kt, 


which determines the values of p, thus: 


eels (16) 


1 See note on Bessel functions. 
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The expression for oe) 5 is readily obtained, 
OZ(p) _ : oy 
pee n sinh y cosh ny ap (17) 


By (15), however, we have 
; Oy 
sinh y Ip = =RC. (18) 


Introducing this value of dZ(p)/dp into (17), we get 


OZ il 
er gn@aRC cosh ry 
nk 
= oF cos (sz). (19) 


For p = 0, find the value of 2: sinh y sinh ny as p approaches 0. 
We have 
cosh y = 1 + 14RCp; 


as p approaches 0, cosh y approaches unity, that is, y is very 
small. Hence, we may write, 


Xe 1 
cosh y = 1 te = 1+ 5kCp, 


and 
7 = RCp. 


Also, for very small values of y, 
sinh y = y; sinh ny = ny 


We have, therefore, as p approaches zero, 


: : eS eas a ae 
2. sinh y sinh ny = Cries nR. 
Hence, 
Z(P)p—o = nh. (20) 


Introducing the values of p and 0Z(p)/dp from (16), (19), and 
(20) in the expansion formula, we obtain the complete developed 
solution for the current in the nth section, 


pein’ 3% 
RC a) 


21 
“ ++ DSS eee cos (sr) oe 


128 HEAVISIDE’S ELECTRICAL CIRCUIT THEORY 


When n is very large, sr/2n small, sin sr/2n = s/2n, and the 
above reduces to the following: 


82472 


i 2k Ron, 
—-+ nR Cos (sr) € 


(22) 


th = 


Comparing this with formula (27), Chap. IV, derived for the 
current in a uniform cable, which is 


np? 
= = +. = S' cos ares 

we note that, forz = 1, that is, at the end of the cable, the expression 
for the current is of the same form as (22), above. If, in one 
case, we consider RI and Cl as the total resistance and total 
capacity of the cable and, in the artificial cable, Rn and Cn as 
the total resistance and capacity, then the two formulas are 
identical. Formula (22), however, was arrived at on the assump- 
tion that the angle sr/2n is sufficiently small to permit replacing 
the sine of the angle by the angle itself; and this is what conditions 
the equivalence of the artificial cable to a uniform cable. The 
number of sections must be sufficiently large so that the sine of 
the angle sr/2n has practically the same numerical value as 
the angle in circular measure. 


NON-DISSIPATING ARTIFICIAL LINE 


For an artificial cable of negligible resistance and leakage, the 
series and shunt elements are 


1 
Shit On Se 
zy DP; 22 ap (23) 
If n is large, the expression for the input current is as before, 
j He 
e 22 sinh Vi 


Also, 

zo sinh y = VW 2122 +142. 
Introducing the values of z; and z2 given by (23), we have, for 
this case, 


ere = 
zesinh y = 4/7 + i Ltpt = “Ph + eH) 
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Hence, 
. 2H 
ie ae ie aay eo (25) 
, 
Lp = Pe 
y= a 
EC 


Expanding the denominator of (25), we obtain 

: aut 1 (40? 1-3/40°\? 1-3-5 /4*\3 } 

Ae es ee a | =. 

= Tal 3 (Ge) alae) — eae ae) + | 
Bearing in mind the operational relation 1/p” = t”/n!, the above 
transforms to 


ee 2H | 5 _ 1 40")? At 1:3 (40°)44 1.3.5 (402)848 saat 
Lp| 2 2! 2-4 4! DPIeG Gl | 
which may be put in this form: 
eres, HP tO 
rts 2? + 52g 2. 42-62 | (26) 


The bracket term is recognized as the series for the zero Bessel 
function and may, therefore, be put in the following convenient 


form: 
2H 


to = Lp J o(2vt) 
t 
a a J o(Qvt)dt. (27) 


0 


The expression for the current is given in the form of a definite 
integral of the Bessel function of the zero order. 

To develop the solution for the current in the nth section, it is 
best to apply the expansion formula. The expression for the 


current in the nth section is - 


22. sinh sinh ny 


tn = 


The determinantal equation 
Z(p) = 2 sinh y sinh ny = 0, 


gives 
sinh ny = 0, 


and 


‘ 8 
ny = jst; 7 = i (28) 
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The values of p corresponding to these values of y are obtained 
from the equation 


Lee 1 21 
cosh y = 1+ tes 
For this case, 
cosh y = cos @ = 1+- SLC p?. (29) 
Hence, 
sae Sir 
{Oy aE) LC 1 — cos ee 
ae ree 
+ j2v sin om (30) 


The expression for 0Z(p)/dp is readily obtained, 


OZ(p) _ E ; ay 
oe nz, sinh y cosh Was 
By (29) 
sinh oo =. LC p. 
Hence, 
OZ(p) _ 
ee nz2LC cosh ny 


nL cos (sr). (31) 


To determine the value of Z(p) for p = 0, take its value as p 
approaches zero. In this case, 


ue epee 
cosh y = 1 + Son 1+ oD’, 


and 
v= Lp 
also, 
sinh y sinh ny = ny? = nLCp’. 
Hence, 


Z(p)p—0 = 22.nLCp? = nLp = 0 for p = 0. 


This would indicate an infinite value of the steady-current 
component for applied steady voltage. This is a consequence 
of the fact that the resistances are entirely neglected. In actual 
practice, of course, the inevitable resistance of the coils would 
limit the magnitude of the steady-current component. We 
need not concern ourselves with it here. What we are interested 


4 
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in mainly is to obtain the solution for the transient-current 
component. 

Introducing the values of p and 0Z(p)/dp from (30) and (31) 
in the expansion formula, we get 


J _ti(2s sin a) 
i, = E> a (32) 
+j 2v sin onet nL cos (sr) 


Taking both terms under the double sign in the above equation 
to include all the roots of the determinantal equation, it reduces 
to the following: 


: 5. Sr 

sin | vt sin — 

Y 2n 
Un — E = 

nLv cos (sr) sin 


; (33) 
ST 


2n 


Equation (33) gives the complete solution to the problem of the 
current in the nth section, but it is in the form of a series which 
is not very convenient for numerical calculations. It is possible, 
however, to transform the solution to another form, that of a 
definite integral of a Bessel function, which is more suitable for 
computations. 

We note that in the expansion formula each term of the summa- 
tion may be replaced by a time definite integral as follows: 


t epst ePst 1 
ji TAC aa 10 10 a 
dZ(p) 
9 <= = Cap ae ata bf Sea et — Ts 
ap P = Ps Ds ap {Ph P dp | es & 
hence, 
eet t EPs 

mig), ~ JoSE@), LAO ae 

8 Op Pp s dp 8 8 ap 8 


Introducing the values of p, and dZ(p)/dp from (30) and (31) 
into the second right-hand term of (34), we get 


if 1 


pees sm 
Raper 


8 Ds +nL cos (sr)2jvt sin 


2n 


which equals to zero, taking into account the double sign. 
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Hence, 


ePst =|l5 eet i 
aZ Ip) Se 
Ps— (2), =e ep — ps 


and using this in the expansion formula, we get 
ePst 
sal Saag (35) 
= (p) 
ep = Ps 


Introducing the values of p, and dZ(p)/dp from (30) and (31) 
we get 


et J2ot sin 


e = | cos jae 
. sr 
cos { 2vt sin “) 


~ 1b cos (sz) 


the summation extending to all values of s from o to n, going by 
steps of z/n. If n is large, the summation term may be con- 
verted into a definite integral, thus: 

Put 
Sees 3, = 993 St = 2n6, 
and the summation term becomes a definite integral whose limits 
of integration are 0 and 7/2, also multiplying numerator and 
denominator by cos (sz) and remembering that cos? (sz) = 1 
for all values of s, equation (86) transforms into the following 
definite integral: 


tn : 
tn = el i cos (2n@) cos (2vt sin @)dédt. (37) 
8 ON Fy be 


The function cos (2vt sin 6), however, can be expanded into a 

Furier-Bessel series,! thus: 

cos (2vt sin 6) = Jo(Quvt) — 2J2(2vt) cos 26 + 2J4(2vt) cos 46 + 
(38) 

If we multiply both sides of (88) by cos 2n@ and integrate from 0 

to 7/2, each term on the right-hand side except the term 


1 See Gray and Maruews, ‘‘ Treatise on Bessel Functions,’ p. 18. 
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Jon(2vt) cos 2n@ will be separately zero. Multiplying the term 
Jon(2vt) cos 2n8@ by cos 2né, and integrating from 0 to 7/2, gives 


2 2 2 
Fag(20b) cos? 2nédé@ = Fan(2vt) [ee = : sin! né 
0 2n 2 4 0 


= 7 J on(20t). (39) 


We arrive, therefore, at the conclusion that 


2 
{ cos 2né@ cos (2vt sin @)d@ = 7 on(20t). 
0 


Hence, substituting in (37) gives 


E t 
in =F | Jen(2at (40) 
L Jo 


This is the equivalent of (33), expressed in an entirely different 
form, the integral of the Bessel function of the 2n order. The 
solution for the current in the zero section (27), arrived at by 
direct operational process, is of a similar form, the integral of the 
Bessel function of the zero order. 

These formulas are more convenient for numerical calculations. 
There are extensive tables of the Bessel functions obtainable, 
from which the values of either Jo(2vt) or J2n(2vt) for different 
values of ¢ are readily obtained. Hence, in using these formulas, 
curves can be plotted, the Bessel functions as the ordinates and 
t as the abscissas extending from 0 to any desired value of ¢, and 
the area enclosed by these curves giving the integrals of these 
functions. 


DISTORTIONLESS ARTIFICIAL LINE 


A brief discussion of the properties and characteristics of a 
distortionless line is given in Chap. V. It is shown that if the 
electrical constants of the line are so related as to satisfy the 
condition R/L = g/C, then the line is said to be distortionless, 
that is, one on which electromagnetic waves are propagated 
without deformation. We shall consider here the problem of an 
artificial line in which the distortionless condition obtains, that is, 
one in which each series element consists of an inductance and 
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resistance, and each shunt element, of a capacity and leakage in 
parallel, and so related as to satisfy the condition 

ls ¥ 

— =. 41 

yb, (Bi ( ) 
If we put, for brevity, 


R_g 
ee een ee 42 
7 = $= %, (42) 
then we have 
a=Lp+R = Lip + 2a), 
ey 1 a ii (43) 
Cp+g Cp + 2a) 
It may be observed in passing, that if the artificial cable closely 
simulates a uniform cable, which is realized when the number of 


£2 


Bren 242 


sections is large, then the propagation constant y should be the 
same as the propagation constant of a uniform line, which we 
designated by K in Chap. V, 


K =VJ/(Lp+ R)(Cp +49). 


When the number of sections is large, L, C, R, g per section are very 
small, and, in that case, from the relation 


1 1 
cosh y=1+5° = 1+ 5(Lp + R)(Cp +9), 


it is evident that cosh y differs from unity by a small quantity; 
that is, y is very small. We may, therefore, put 


cshy=1+Wy=14+40p4+ Rg t+ Cp), 
and 


y = V(Lp + R)(g + Cp), 
the same as K. 

We shall confine the discussion, in this case, to the determina- 
tion of the current in the mth section. We have here, as in the 
previous case, the general expression 

E 
z. sinh y sinh ny 


la = 


4 
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The determinantal equation is, in this case, also, 
sinh ny = 0, 
and 
: . St 
CG cle EY SID faa ie 
The values of p corresponding to the roots of the determinantal 
equation are determined from the relation 
1 21 
hiv = — 
cosh vy 1+ 9 2,’ 


which, for this case, gives 


cosh y = cos a = 1+ 5 LC(p + 2a)?. (44) 
Hence, peetenne Ce 
p = —2a + jon)2(1 — cos 2) 
n 
eam, : 
Se ae a < a : 
2a + j2v sin On (45) 


We also have here 


— = nz sinh y cosh ny = (46) 
By (44), 
Shi = eon 
Op 
Introducing this value in (46), we get 
ee = nzoLC(p + 2a) cosh ny 
= nL cos (sr). (47) 
For p = 0, Z(p)p=0 = 2 sinh y’ sinh ny’ where 2 = 1/2Ca; and 
y' is determined from the relation 
cosh 7’ = 1 + 2LCa?’. 
In the expansion formula, #/Z(p),-0 gives the permanent- 
current component. This, however, is not of any particular 


interest here. What concerns us chiefly is the transient-current 
component. This is given in accordance with the expansion 


theorem by 
: eP 
ae PRO) 
dp 
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Substituting the values of p and 0Z(p)/dp from (45) and (47) 
results in the following: 


s7 
+ j(20t sin — ) 
Hy € 2n 
| 


(48) 
nL cos (sn) + 72v sin sill as 
( 


2n 


Taking both terms under the double sign to include all the roots 
of the determinant, the above is expanded into 


af PSE 
_ Hens Se (20 sin ) + 7 sin (20 sin =) 
(2 sin a — 22) cos (sz) 


. ST sin ; S71 
cos (20 sin ) — 7 si (20 sin #4 


_ > ; »+ (49) 
(20 sin 3 sin - =f ) cos (sz) 


and this simplifies to 


eae \ ies ~ aR Pak 
. ren sin 7 sin (20 sin zs — 2a cos (20 sin ) 


(40 sin? = ou 10°) cos (sr) 


1 
. nL 


(50) 


This is the complete developed solution for the transient current 
in the nth section. The solution is in the form of a series which 
is rather complicated and would be very laborious to use in 
computations. It is possible, however, in this case, also, to 
obtain a solution in the form of a definite integral of a Bessel 
function, as in the case of the problem of the non-dissipative 


et 


04 (P) 0) P 
* Op 


artificial line. Asin the previous case, we can replace —_~— 


t 
[a@ ee and write 
dp 


E i Se” (51) 
Op 
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ODN ee the values of p and dZ(p)/dp, we get 
(20 sin =) 


=) > nL cos (sr) ae 
e«” cos (20 sin =) 
soe Ah » an ay (52) 


~ nb cos (sz) 


We have, however, already established that 


2 cos (20 sin a 
n Cos (sr) upon) 


(see equations (36) to (40)), therefore, 


E t 
i,== i «Jy, (Qut) dt, (53) 
L Jo 


the same as (40) modified by an attenuation factor e~%. 

Series Elements: Inductance and Resistance. Shunt Ele- 
ments: Capacities——The investigation of the properties and 
characteristics of an artificial line in which the leakage elements 
are not included is of considerable importance, inasmuch as this 
is the condition resembling more closely an ordinary transmission 
line in which leakage is generally a negligible factor. Also, this 
type of artificial line is very closely related in its properties and 
performance to the periodically loaded line. The solution of 
this problem, therefore, will serve also as an introduction to the 
study of the loaded-line problem. For a comprehensive mathe- 
matical discussion of the loaded line, reference must be had to 
Dr. Pupin’s celebrated papers on this subject.! 

For an artificial line of negligible leakage 


ait = Lp + 1% 
1 (54) 
Zo Cp’ 
and 
cosh y = 1+ : . =1+ sCpLp 4 P), (55) 


1 Trans. Am. Inst. Elec. Eng., Vol. XVI, 1899; Vol. XVII, 1900. 
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The determinantal equation is, of course, the same as in the 
previous cases. Hence, we have, as before, 


ea a 
We have, therefore, 


cos - == J] 5 (LCp? + RCp), 


from which the values of p are determined as follows: 


sr 
\iece—stc(1 — cos x) 


R 
5 a = 
f 2L 2LC 
= — } 2 ain2 or aa OL 
ac in| sin 5 (56) 
where, for brevity, 
Pd Bg / LC 


We also have, as before, 


OZ(p) _ Oy 
Saas nz sinh y cosh ops 


but, in this case, 


sinh 15. = LCp + BRC. 


Hence, 
Z 
mete es na LCp + aRC) cosh ny 
= m0 + a) cos (sr). (57) 


For p = 0, we take the value of Z(p) = 22 sinh y sinh ny as p 
approaches zero. For small values of p, cosh y differs from unity 
by a very small quantity only, and we may, therefore, put 


cosh y = 1+ 147? = 1+4 14(LCp? + RCp) 
Therefore, 
vy’ = LCp’ + RCp, 
and 
n(LCp* + RCp) 
mR. (58) 


ie 


Z(p)p=0 = 22 Sinh y sinh ny = cam 


ARTIFICIAL LINES 139 


That is, the steady-state component is simply E/nR. For 
the transient-current component, we have 


pt 
ee pt Noe 
: ue 

dp 


= ED ar + a) cos (sr) (59) 


Introducing the values of p from (56), we get 


. Ss i 
+j Vao2 ain? eget 
€ 2n 


in = Een#! (60) 


: meets 
+ ‘ 2 paras ae we 4 
+ jnL cos (ox) | sin’ 5 a 
Taking both terms under the double signs to include all the values 
of p corresponding to the roots of the determinantal equation, 
the above simplifies to the following: 


; ated ©). as 
fa Nae sin ([ sin Qn at) 
SG S 5 Sir , 

cos (sr) 4 /4v? sin? ae a 


The current is the sum of n oscillatory components, all having 
the same damping factor, but each of different frequency. The 
frequencies of oscillations are given by 


x= sh, Bieaeyay tee te 2 62 
i= 7/40 sin’ 5 a (62) 


If n is infinitely large, the condition of a uniform line is 
approached, and in that case the solution given by (61) should 
reduce to that of the uniform line. That this is actually the case 
can be readily verified by comparing the above solution with (80), 
Chap. V, the solution for the oscillations on a uniform line. 

If n is very large, we may put 


(61) 


Afi ee 

Qn  2n 
and (61) reduces to 
; ( vs? «) f 
sin ——-—a 
PAYER n? 
meer D S pane Co 

cos (sz) —@ 
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Referring back to equation (80), Chap. V, we note that, 
neglecting leakage and putting x = /, the expression for the 
current at the end of the line reduces to the following: 


; 2H? ; sin Bpt 
= — Le ——— ” 4 
11 i Cer 3. cos (nr) (64) 
For g = 0; 
pete 
p aL = ) 
ny? 
Br a Tes —~— 
ig aaa 
i ap Fo 
and 


a 


sin( posts - or 
2 
s a> Ewe cos (nz): (65) 


If we bear in mind that n, the total number of sections in (61), 
is the equivalent of 1, the length of the line, in (63), and that s 
in (61) and n in (63) designate the same thing—the summation 
steps—then it is evident that the two formulas are identical. 
The solution for the uniform line can be derived from that of the 
artificial line by making the number of sections infinitely large. 
The solution of the uniform line is a limiting case of that of the 
artificial line. 


ARTIFICIAL LINE-SERIES ELEMENTS INDUCTANCE AND RESIST- 
ANCE; SHUNT ELEMENTS CAPACITY AND LEAKAGE 


The basic formula (13) for the current in the nth section 
applies, also, of course, to this case; that is, 
E 


‘n= 2 sinh y sinh ny, CES 
and 
es 
Ves Las 


The values of p to be used in the expansion formula are here 
also determined from the relation 
1 21 


cosh y = 1 + 
2 22 
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In this case, however, 21 = (Lp + R) and z2 = pea hence, 


Co+g 
cosheyi = coe e =1+3Lp+R)(Cp+g), (66) 


which may be put in this form: 


1 
cos ) =1+ apt?” + 2(a + b)p + 4ab}- 


Solving for p, we obtain 


—2(a + b) + /4(@ + b)? — 16ab — sv(1 — cos ) 
[Os 2 ? 


= bau Se 7 2 ST 
(a + b) + /(a b) 4? sin On’ 


=-pt jaf Woe a (67) 
We may put, for brevity, 
where aeeane 
Bs = 4/40 sin — o* (68) 


The expression for 0Z(p)/dp to be used in the expansion formula 
is readily obtained, 


OZ(p) _ ; a(ny). 
= ih z. sinh y cosh ny ine (69) 
By (66), however, 
: Oy ek 
Sinh ane (2LCp + Lg + RC) 
1 1 
= a +a b) = Gp Fe): (70) 


Introducing this value in (69), we get 


22) = Nee B (p + p) cosh ny 
pap). 
= GES cos (sn) 


+ njB. +njB; cos (sm). 


= i {f = ‘ 1 
ee ee oo peas By 
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Substituting these values of p and 0Z(p)/dp in the expansion 
formula, we obtain the following: 

EB (— ot 9Bs)ei— 2 £780)t 
nL +78. (— p + jB;) cos (sr) 
Taking into account the double-sign terms to include the values 


of p corresponding to all the roots of the determinantal equation, 
we obtain the complete expression, as follows: 


gee eee > (=o tee eget 
; jBs(—p + 98s) cos (st) 7Bs(p + jBs) cos (sr) 


mF gn Sie tH — w= HBedePt = (© + 18) (0 = 5Be)e Ht 
nL j8.(0? + B3) cos (sr) 


_ 2He7 se a — p) cos Bt + (op + 6%) sin Bot 
Bs(p” + 82) cos (sm) 
This is the complete solution for the transient-current com- 
ponent in the last, the nth, section of an artificial line which is to 
simulate a line in which all the electrical constants, R, L, C, g, 
are active. 
Equation (73) may be put in a little more convenient form, 


(72) 


tn = 


(73) 


thus: 
; oe bes + 6? cos (Bt — Ys) 
tae -> eee + B2 B, cos (sr) (74) 
op Ff, 
tan ¢, = (oes ay (75) 


Formula (73) covers the most general case and should reduce 
to the special cases discussed before by the proper choice of 
values of R, C, L, g. 

For R = g = 0, the formula should reduce to (33), the expres- 
sion for the current in a non-dissipative line. For this condition, 


p = 0,0 = 0; B = 2 sin = and (73) reduces to 


: . sr 
‘ on sin (so sin -) 
Pen : 


. sr 
2v sin —— co 
an 8 (sr) 


which is formula (33). 
For the ba cas condition, = 0, p = 2a, 


GB ==02y sin * of ; and for these values (73) reduces to (50). 
For g = 0, p = 2a; 0 = 2a and (73) reduces to (61). 


CHAPTER VII 
HEAVISIDE’S DERIVATION OF EXPANSION FORMULA 


It will be of interest to give an outline here of the method by 
which Heaviside arrived at the expansion formula. Aside from 
the historical interest, the physical and mathematical ideas by 
which he was led to the formulation of this important theorem 
should be of considerable interest for the student of physics 
and engineering. Heaviside was profoundly interested in the 
study of the dynamics of an electromagnetic field; there are 
frequent recurrences to this subject all through his collected 
papers. Irrespective of what particular investigation in electrical 
theory he was working on, the subject of the dynamics of an 
electromagnetic field, and related matters, is touched upon from 
one angle or another. It was the study of energy distribution 
and subsidence in an electromagnetic system that suggested 
ideas which led to the formulation of the expansion theorem. 
Nowhere has he given a complete and formal demonstration 
of the derivation of the formula. But ideas and suggestions 
relating to this and connected problems run through the publica- 
tions of his collected papers. 

In his “Electromagnetic Theory,’! he gives the formula 
without proof. He contents himself with the following state- 
ment, and no reference is made to the process of derivation: 


Finally, there is a third and very general way of converting operational 
solutions to the form of the sum of normal solutions. It does not 
require special investigations of the properties of normal functions. It 
is very direct and uniform of application. It avoids, in general, a large 
amount of unnecessary work. The investigation of the conjugate 
property, and of the terminal apparatus in detail in order to apply it to 
the determination of the coefficients, is wholly avoided. It applied to 
all kinds of series of normal functions, as well as Fourier series. And it 
applies generally in electromagnetic problems, with a finite or infinite 
number of variables; or more generally, to the system of dynamical 
equations used by Lord Rayleigh in the first volume of his treatise on 


1 Vol. II, p. 127 
143 
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Sound, which covers the rest of the work, and upon which he bases his 
discussion of general properties. 

The method may be briefly (though imperfectly) stated as follows: 
Let e = ZC be the operational solutions of an electromagnetic problem; 
say, for definiteness, that C is the current at a certain place due to an 
impressed force e at the same or some other place. Let the form of Z 
be such as to indicate the existence of normal solutions for C. Then 
when e is steady, beginning at the moment ¢ = 0, the C due to e is 
expressed by 


aoe er 
Cee an 


This is recognized as the expansion formula which we have 
developed by algebraic process; the complete derivation is 
given in Chapter II. 

Because of the fact that in giving this formula in his electro- 
magnetic theory, no reference is made to his collected papers, 
nor is any hint given as to the how and when of its derivation, 
the impression was created that Heaviside gave this formula 
without proof. This, of course, is incorrect. Heaviside 
discussed the subject very extensively but in scattered form.! 

The basic ideas which led to the development of the expansion 
formula relate to questions of energy distribution and subsidence 
in a normal electromagnetic system. A normal system is defined 
as one which on subsiding remains similar to itself, the subsidence 
being represented by the time factor e?’. 

If we have any electric-circuit system, energy stored up 
magnetically and electrically in the inductances and condensers, 
and leave it to itself, removing the impressed force, it will subside 
to equilibrium in a manner determined by the distribution of the 
currents in the coils and the charges of the condensers. In 
general, the energy relations are as follows: 


W+SU+T) =0, (1) 


where U is the electric energy T the magnetic energy, and W the 
energy dissipated. This formula states that the rate of decrease 


1 Tt is partially dealt with in his papers, ‘Electromagnetic Induction and 
Its Propagation,’’ collected papers, Vol. I, pp. 429-560; ‘‘On Self Induction 
of Wires,”’ collected papers, Vol. II, pp. 168-323; and the final derivation of 
the formula is given in his paper “‘ Resistance and Conductance Operators,” 
collected papers, Vol. II, pp. 355-374. 


HEAVISIDE’S DERIVATION OF EXPANSION FORMULA 145 


of the electric and magnetic energies of the system is equal to 
the dissipativity. 

When the energy subsidence occurs in two normal systems, 
as in the case, for instance, of a double periodic oscillating circuit 
system, and if we designate the two systems by suffixes ; and ., 
an additional energy relation is given for this case by the equation 


Wie + Se: + Ti2) = 0, (2) 


an equation of mutual acticity. Ui, and T1. are the mutual 
electric and magnetic energies, and W,,. the mutual dissipativity. 
The mutual energies are the excess of the total energies when the 
two normal systems coexist over the sum of the separate energies 
if they would have existed independently. In addition, for a 
system in which no energy is communicated, leaving it only 
irreversibly through dissipativity, there is also the relation 
Ure — T 12. (3) 
The mutual electric and magnetic energies are equal. Heaviside 
established these relations from general considerations of a 
dynamical system, which will be given farther on. We shall 
first consider these questions from the electric-circuit standpoint, 
following the method of Wagner.! 
Take the simplest case, that of two circuite coupled through 
a common element, and consider the three possible cases— 
capacity coupling, inductance coupling, and resistance coupling. 
a. Capacity Coupling.—The circuit arrangement is shown in 
Fig. 25; the circuits are coupled through the condenser C4». 


BiGse2o: 


If the charge on condenser C, is Q:, and the charge on condenser 
Cz is Qe, the charge on the coupling condenser Cz is Qi + Qe. 


1Waaner, Karu Witty, “Der Satz von der Wechselseitingen Energie,’’ 
Elektrische Nachrichten-Technik, Bd. 2, 8. 376-392. 
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Zi, and Z2 are the impedances of the circuits 1 and 2 exclusive 
of capacities. We have the circuit equations 


a Qi + Qe 


le EVAR i SES = E,, 
0, Q +0 me 
a la to Ey 
Put 
1 1 
sa eer ee 
1 1 
Ko. — Cs + Ce (5) 
1 
= K = ——, 
Ki. 21 Cis 
Also, 
d d 
I, = a = pQi; In = oe = pQ:. 
Substituting these in equations (4), we get the following: 
(Ku aia pZir)Q1 ae K12Qs = E,, | (6) 


(Ko. + PZ22) Qo + KuQ = Bp. } 


b. Inductive Cowpling—Two circuite coupled through an 
inductance Li., the arrangement shown in Fig. 26. If we put 


Diu = Ty = D2 


Fie. 26. 


the total inductance of circuit 1, and 
Le. = Le = Tie 


the total inductance of circuit 2, we have the following equations 
for the two circuits: 


(pZ11 + p?Li1)Q: + p?Li2Q> Fi, (7) 
(pZ22 + p*L22)Q2 + p?LeiQi E2, 


Zi, and Zs» are the impedances of the circuits exclusive of 
inductances, 
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c. Resistance Coupling.—The circuit arrangement is shown in 
Fig. 27. Ris, is the coupling resistance, Z1,; and Zy». are the 
impedances of the circuits exclusive of resistances. If we put 
Ri = Ri + Riz, the total resistance of circuit 1, and Ra = 
R2 + Ro the total resistance of circuit 2, the circuit equations 
will be as follows: 


p(Zi1 —- Rii)Qi at PRi2Q2 = fii, (8) 
p(Zoe = E2)Q2 te pi = EH. | 


RiGe 2. 


In the general case for any number of circuits and all three 

types of coupling simultaneously effective between the circuits, 

the circuit equations in accordance with the preceding will be 

as follows: 

(Ki + pit + p11) Q1 + (Ki. =f pRiz ae p?L12)Qo+ Sere” 
ss Ey, 


(Kee = DR» + p?Li22) Qe aa (Kor Se pRo1 + p?L1)Q1 See 
— Eo. 


(9) 


We also have the following expressions for the energies in the 
electric circuit system. 
The electric energy in the condensers 


U= WKuQ IPs K12Q1Q2 =F 16K 05 amet Pe (10) 


The magnetic energy in the inductances 


T= LyT? + Lyell, + L6LoI? esis (11) 
The energy dissipated by the resistances 
Wea Ryul? + 2Rio1l WwW. + Roel? to-+-- (12) 


The power, rate of energy, supplied to the circuits by the 
impressed electromotive forces, 
S = Hil, + Helo + ea eine (13) 
By the aid of the above formule, we can establish the relation 
between mutual electric and magnetic energies, 
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Assume two systems of electromotive forces acting on a 
circuit network. One system we shall distinguish by ’ thus: 
Ei, Ei, ..., the corresponding charges by Q), Q3, . . 
and the corresponding currents by I/, Jj,... We shall 
also designate by U’ the electric energy; T’ the magnetic energy; 
W’ the dissipativity; and S’ the power supplied. 

A second system of electromotive forces, different from the 


first, we shall distinguish by ” thus: H’{, HZ .. . ; the charges 
corresponding to these electromotive forces are designated by 
NoiQas. si. y the: currents by slaty. se Oo aeeiy = 


nated by U”, the electric energy, 7T’’ the magnetic energy, W’’ 
the dissipativity, and S’’ the power supplied. 

When both systems of electromotive forces act simultaneously 
on the circuit network, the charges and the currents are additive 
so the effective charges are Qi + Q7, Q; + Q2 ... and the 
effective currents are [i +I, I,+ 7, ... The energies, 
however, of the two systems are not additive, because the 
energies are proportional to the squares of the charges and the 
currents. The resultant energies-are as follows: 


U =U +" + Uys, | 
THT 4T' 47, | 
W=W LW" + Wr, 
S = 8S’ + 8” + Syo. 


(14) 


U2 designates the mutual electric energy between systems | and 2, 


Ure = KiQ,Q07 + Ki2(QiQ5 + Q’Q)) + Ko2QiQ7 + --- (15) 

Ti. is the mutual magnetic energy between systems 1 and 2. 

Pre = Lull) + Lal” + 5) + Deli’ +: (16) 

The corresponding mutual dissipativity is 

Wi = 2{RuTiT” + Rel’, + I) + Reali} +--- 7 

The corresponding power, rate of energy, supply is given by 
Sv= (AI, + BD + (A + BUD +... (18) 


Of the mutual power supply, the part contributed by system 
1 is 

Sp = ET! + BI, +--: (19) 
We will now assume that in system 1, all amplitudes vary as the 
exponential function ¢? and that in system 2, all amplitudes 


tf 
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vary as the exponential function ¢?:*. We shall have then, by 
(9), 
kyQ) + Rul, + pila, + k12Q) + Rooly + piliel + 


I’; = peQy’, ete. 
Introducing these values in (19), we obtain 
Sip = pUw + mle + Wx. (21) 
In a similar way, from the relation 
eee lag tees tetas 
we get 
Si= piU 12 ae p2T 12 ata 16 Wr. (22) 
These are the equations of mutual activity. Subtracting (22) 
from (21), we get 
Sy — Su = (pe — pi)(U2 — Tr). (23) 
Now let the H’s vanish, so that no energy can be communicated 
to the system, then S/, = 0, and Si} = 0, and 
(p2 = pi)(Ui2 — Ti) = 0, 
giving 
Uie — T 12. (24) 
By the aid of this relation between the mutual energies, the 
amplitudes of the normal functions of a system of any number 
of degrees of freedom can be determined. The expansion 
formula was derived by the application of this principle to 
electric-circuit systems. 

Heaviside established this relation between the mutual energies 
from the general considerations of a dynamical system. He 
states the whole problem very clearly and concisely in his paper 
as follows: 


THE CONJUGATE PROPERTY Uy. = 7): IN A DYNAMICAL SYSTEM 
WITH LINEAR CONNECTIONS 


Considering only a dynamical system in which the forces of reaction 
are proportional to displacements, and the forces of resistance to 
velocities, there are three important quantities—the potential energy, 
the kinetic energy, and the dissipativity, say U, T, and Q which are 
quadratic functions of the variables or their velocities. When there is 
no kinetic energy, the conjugate properties of normal systems are 
Ui. = 0 and Qi2 = 0, these standing for the mutual potential energy 


1“On the Self Induction of Wires,”’ collected papers Vol. II, p. 202, 
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and the mutual dissipativity of a pair of normal systems. When there 
is no potential energy, we have 712 = 0 and Qi2 = 0. When there is 
no dissipation of energy, Ui2=0 and 71. = 0. And, in general, 
Ui. = T12, which covers all cases and has two equivalents, 1Qi2 + 
Uy. = 0, and Qi. + Ti. = 0; for, as the mutual potential and kinetic 
energies are equal, the mutual dissipativity is derived half from each. 

Let the variables be 71, x2 - - * , their velocities 7: = 4, - - - and the 
equations of motion 


F, = (Au + Bup + Cup?) a1 + (Arie + Bip + Crp?) v2 + 


@} 
88 
= (Aga + Bop ae Cop?) ayt (A o2 = Boop + C oop?) totee: | ) 
where F,,F: - - + are impressed forces and p stands for d/dt. Form- 
ing the equation of total activity, we obtain 
ZFu=Q+U+4+T (89) 
where 
VUn— Axi + 2A 10% 1X2 — Art? a sit Bry 
Q = Bui + 2B yo 4 Bow? — Pee (90) 
Did ee Cy} + 2C 2012 + C203 +++. 


So far will define in the briefest manner, U, 7, Q and activity. 

Now let the F’s vanish, so that no energy can be communicated to the 
system, whilst it can only leave it irreversibly, through Q. Then let 
Pi, P2 be any two values of p satisfying (88) regarded as algebraic. Let 
Q:, U:, T; belong to the system p; existing alone; then, by (89) and (90), 


0=Q+ Ui + 11, or 0 = Qi + 2p1(Ui + 73); 
0O=Q.+ U2 + Te, or 0 = Qo + 2p2(U2 + T°). 


But when existing simultaneously, so that 
Q=Q+0.:4+ Qn, U = Ui+ U24+ Uy T= T+ 724+ Tr, 
where U2, T12, Qi2 depend upon products from both systems, thus: 


Qie = 21Bivw, + Bow, + Bio(ow, + vw), + - °° 3 
Oy. = Aiwa, + Aster, + Ar(aia, + wet,) +s: 
Tie = Civ, + Covyoy + Cro(vivy + vv) + > 


the accents distinguishing one system from the other, we shall find, by 
forming the equations of mutual activity, 2F,;=---, and [Fi = 
, that is, with the F’s of one system, and the v’s of the other, in 
turn, 
= 4002 + poU 2 + pili, 
= Qi + piU ie + pol 9; 


adding which, there results the equation of mutual activity, 
0 = Qi + (pr + pe) (Uie + Ts) or O = Qire + Ure + Pro; 
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and, on subtraction, there results 
= (pi — p2)(Ui2 — Ti), (91) 


giving Uj. = Ty. if the p’s are unequal. But this property is true 
whether the p’s be equal or not; that is, Ui: = 71, when pi: isarepeated 
root. I have before discussed various cases of the above, with social 
reference to the dynamical system expressed by Maxwell’s electro- 
magnetic equations. 


The above quotation gives Heaviside’s approach to the problem. 
He arrives at the relation of the mutual energies given by (24), 
Heaviside’s equation (91), from the general energy considerations 
of a dynamical system. 


DETERMINATION OF THE AMPLITUDES OF NORMAL DISTRI- 
BUTIONS OF VOLTAGES AND CURRENTS IN A CIRCUIT 
NETWORK 


Assume a network consisting of n circuits, each comprising 
an inductance and a condenser; energy stored up electrically 
and magnetically in the circuit system. On removal of the 
impressed electromotive force, the system will subside to equilib- 
rium in normal current and voltage distributions, the subsidence 
being represented by the time factor e?’.. The voltage and current 
in any circuit is given by 

V = 2Aue’, 
I = DAwe*’, 


V being the real voltage at a place where the corresponding 
normal voltage is wu and J the real current where the correspond- 
ing normal current is w. 

For any particular normal distribution, the amplitudes of the 
voltage across the n condensers are 


Ayu, A,yu®, Beis A,u™) 


where A, is the coefficient of the normal distribution correspond- 
ing to the vth distribution. 
The amplitudes of the currents in the inductance coils are 
A,w?, A,w®, ~~» A,w®. 
The effective voltage on any one condenser, of the pth circuit, 
say, is given by 
Vu = Ayu yt + Anu ot +--+: (25) 
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The current in the coil of the uv” circuit is given by 

I, = Aywt eit + Apwh) 2 + --- (26) 

Initially, when ¢ = 0, the voltage and current in the uw” circuit 
are given by 

Vio = Au? + Agu +--+ | 


2 
Tyo = Aw + Aw) +--- (27) 


Now calculate the mutual electric and magnetic energies of 
the given state with respect to, say, the v’* normal distribution. 
To do this, multiply (27) by C“u% and sum up over all the n 
circuits. We obtain the following, leaving off the designation 
u, and using the summation sign: 


Uo = LCuoVo = TC{Auitly + Aotiatty ++ -- } 
= Aa Ui» + AUey + ft Poe (28) 


The terms U;,, on the right-hand side of the above equation are 
the mutual electrical energies of any K‘’ distribution and the 
vy’ normal distribution. 
In a similar way, we obtain the following expression for the 
magnetic energies: 
To = ArT + AcT2, + * °° (29) 


Subtracting (29) from (28), all terms Ax(U;z, — Tx.) on the 
right-hand side cancel by the relation (24) except the term 
A,(U,, — T,,). The energies U,, and T,, are not mutual 
energies. U,,, is double the self energy of the v* normal distribu- 
tion, and so is T,, = 2T,. 

We obtain, therefore, the relation 


Uo, ae Tow = ZAG, a Pe), 
and 
= Cea De 
~ 20, =T.) 
This relation between the amplitudes of the normal distributions 
and the energies in the initial state was developed by Heaviside 
in his paper, ‘‘ Electromagnetic Induction and its Propagation, ” 
collected papers, Vol. I, p. 523. 

From this and another relation that he establishes between the 
energies of the system and the electrical characteristics of the 
circuit system, the energy terms are eliminated and an expression 
is obtained for the amplitudes of the normal distributions in 
terms of the circuit characteristics and the initial states. 


A, (30) 


fe 
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The derivation of formula (30) is not set forth in Heaviside’s 
papers connectedly in a way that could be followed easily; 
hence, the demonstration given above. From this point on, 
however, the rest of the work in the process of deriving the 
expansion formula is very clearly given in his papers, and it 
will be best to quote him directly.! The two sections referred 
to, together with formula (30), give the complete derivation of 
the expansion formula. Parts of these sections which have 
specific reference to the problem are given below. 


Suppose, for example, we have two fine-wire terminals, a and b, that 
are joined through any electromagnetic and electrostatic combination 
which does not contain impressed forces, nor receives energy from with- 
out, except by means of the current, say C, entering it at a and leaving 
it at b. Let also V be the excess of the potential of a over that of b. 
Then VC is the energy-current or the amount of energy added per second 
to the combination through the terminal connections with, necessarily, 
some other combination. . . . The combination need not be of mere 
linear circuits, in which differences of current-density are insensible; 
there may, for example, be induction of currents in a mass of matter 
either connected conductively or not with a and b; but in any case it is 
necessary that the arrangement should terminate in fine wires at a and 
b, in order that the two quantities V and C may suffice to specify, by 
their product, the energy current at the terminals. Even in this we 
completely ignore the dialectric currents and also the displacement, in 
the neighborhood of the terminals, 7.e., we assume c = 0 to stop dis- 
placement. This is, of course, what is always done, unless specially 
allowed for. 

Now, supposing the structure of the combination to be given, we can 
always, by writing out the equations of its different parts, arrive at a 
characteristic equation connecting the terminal V and C. For instance, 


V =Z0 (98) 


where Z is a function of d/dt. In the simplest case Z is a mere resist- 
ance. A common form of this equation is 


foV + AV +f2V + owes = mC + ac +90. + 506 


where the f’s and g’s are constants. But there is no restriction to such 
simple forms. All that is necessary is that the equation should be 


1 It is necessary to refer to only two sections, one given on p. 204, Vol. II, 
collected papers, under the heading “‘ Applications to Any Electromagnetic 
Arrangement Subject to V = ZC,” and the other section on p. 371, Vol. II, 
collected papers, under the heading “‘The Use of the Resistance-operator in 
Normal Solutions.”’ 
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linear, so that Z may be a function of p. If, for example, (dC/dt)? 
occurred, we could not do it. 

Now this combination must necessarily be joined on to another, how- 
ever elementary, to make a complete system, unless V is to be zero 
always. The complete system, without impressed forces in it, has its 
proper normal modes of subsidence, corresponding to definite values of 
p. Consequently 


U2 —T12 = (U2C1 —U C2) + (pi— pa)! (99) 


if Vi,C, belong to p; and V2, C2 to p», whilst the left member refers to 
the combination given by V = ZC. Or, 


Vi } i741 — Le 
Uye — Tre. = CiC2 — — —) + — =) Cp — LOO 
12 12 1U2 an C, (po Pi) 1 Oe a, ( ) 
and the value of 2(U — T) is a single normal system is 
a Ede Ar Ve OV Se eagle 
2(U — T) Gr oop = —C Fr a 6 a6 (101) 


In a similar manner we can write down the energy differences for the 
complementary combination, whose equation is, say, V = yC; remem- 
bering that — VC is the energy entering it per second, we get 


¥: dY 


SESE rear and ae respectively. 
By addition, the complete Ui2 — T12 is 
CC a i ee See (102) 
Disa JIN UO 
and the complete 2(U — T) is 
d dg 
C2 (Veo ==, 
FP ) or de (103) 


where g = 0, or Y — Z = 0, is the determinantal equation of the com- 
plete system (both combinations which join on at a@ and b where V 
and C are reckoned), expressed in such a form-that every term in is of 
the dimensions of a resistance.” 


THE USE OF THE RESISTANCE-OPERATOR IN NORMAL 
SOLUTIONS 


In conclusion, consider the application of tne resistance-operator to 
normal solutions. If we leave a combination to itself without impressed 
force, it will subside to equilibrium (when there is resistance) in a 
manner determined by the normal distributions of electric and magnetic 

1 See equation (23), p. 149. 

* Collected papers, Vol. II, pp. 204-206 
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force, or of charges of condensers and currents in coils; a normal system 
being, in the most extended sense, a system that, in subsiding, remains 
similar, the subsidence being represented by the time factor e?t, where 
p is a root of the equation Z = 0. It is true that each part of the com- 
bination will usually have a distinct resistance-operator; but the 
resistance-operators of all parts involve, and are contained in, the same 
characteristic function, which is merely the Z of any part cleared of 
fractions. It is sometimes useful to remember that we should clear of 
fractions, for the omission to do so may lead to the neglect of the whole 
series of roots; but such cases are exceptional and may be foreseen; 
while the employment of a resistance-operator rather than the char- 
acteristic function is of far greater general utility, both for ease of 
manipulation and for physical interpretation. 

Given a combination containing energy and left to itself, it is upon 
the distribution of the energy that the manner of subsidence depends, 
or upon the distribution of the electric and magnetic forces in those parts 
of the system where the permittivity and the inductivity are finite, or are 
reckoned finite for the purpose of calculation. Thus conductors, if 
they be not also dielectrics, have only to be considered as regards the 
magnetic force, whilst in a dielectric, we must consider both the electric 
and the magnetic force. Now the internal connexions of the system 
determine what ratios the variables chosen should bear to one another 
in passing from place to place in order that the resulting system should 
be known; and a constant multiplier will fix the size of the normal 
system. Thus, supposing wu and w are the normal functions of voltage 
and current, which are in most problems the most practical variables, 
the state of the whole system at time ¢ will be represented by 


V = LAuert, C = DAwer; (46) 


V being the real voltage at a place where the corresponding normal 
voltage is u, and C the real current where the normal current is w, the 
summation extending over all the p-roots of the characteristic equa- 
tion. The size of the systems, settled by the A’s (one for each p) are 
to be found by the conjugate property of the vanishing of the mutual 
energy difference of any pair of p-systems, applied to the initial dis- 
tributions of V and C. 

To find the effect of impressed force is a frequently recurring problem 
in practical applications; and here the resistance-operator is specially 
useful, giving a general solution of great simplicity. Thus, suppose we 
insert a steady impressed force e at a place where the resistance-operator 
is Z, producing e = ZC thereafter. Find C in terms of e and Z. The 
following demonstration appears quite comprehensive. Convert the 
problem into a case of subsidence first, by substituting a condenser of 
permittance S, and initial charge Se, for the impressed force. By 
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making S infinite later we arrive at the effect of the steady e. In 
getting the subsidence solution we have only to deal with the energy of 
the condenser, so that a knowledge of the internal connexions of the 
system is quite superfluous. 

The resistance-operator of the condenser being (Sp)-!, that of the 
combination, when we use the condenser, is 7:1, where 


Li Sp) 2. (47) 


Let V and C be the voltage and the current respectively at time ¢ after 
insertion of the condenser, and due entirely to its initial charge. Equa- 
tions (46) above express them, if w and w have the special ratio proper 
at the condenser, given by 


w= —Spu (48) 


because the current equals the rate of decrease of its charge. Initially, 
we have e = 2Au and ZAw = 0. So, making use of the conjugate 
property, we have 


Seu = 2(U, — T,)Av3 (49) 


if U, be the electric and 7, the magnetic energy in the normal system. 
But thefollowing property of the resistance-operator is also true. 


AT a Us) = ee (50) 


That is, dZ,/dp is the impulsive inductance in the p-system at a place 
where the resistance-operator is Z;, p being a root of Z; = 0; just as 
dZ,/dp with p = 0 is the impulsive inductance (complete) at the same 
place. Using (50) and (49) gives 


dp 


Now use (48) in (51) and insert the resulting A in the second of (46), 


and there results 
(ee >s ept 52 
pZ\ eh 


where the accent means differentiation to p. This is the complete 
subsidence solution. Now increase S infinitely, keeping e constant. 

Z, ultimately becomes Z; but in doing so one root of Z; = 0 becomes 
zero. We have, by (47), and remembering that Z,; = 0 


pZ, = —(Sp) + pZ' = Z + pl’; (53) 


1See equation (30), p. 152; the initial energy uo, condenser charge no 
initial magnetic energy. 
2 See equation (103) p. 154. 


Ae iseaee (w (51) 


4 


HEAVISIDE’S DERIVATION OF EXPANSION FORMULA 157 


so, when S = » and Z = 0, we have pZ; = pZ’ for all roots except the 
one Just mentioned, in which case p tends to zero and Z’ is finite, making 
in the limit pZ, = Z» by (53), where Zp is the p = O value of Z, or the 
steady resistance. Therefore, finally, 


oS ee (54) 
where the summation extends over the roots of Z = 0, shows the manner 
of establishment of the current by the impressed force e. The use of 
this equation (54), even in comparatively elementary problems, leads to 
a considerable saving of labor whilst in cases involving partial differen- 
tial equations, it is invaluable. 


It is interesting to compare the methods of deriving the 
expansion formula as originally given by Heaviside, and the 
method developed in Chap. II. It is observed that Heaviside 
starts from certain fundamental considerations of energy distribu- 
tions in an electric-circuit system, first arriving at an expression 
for the energy subsidence when the system is left to itself, no 
external electromotive forces acting on it, and then developing 
the more general expression for the current rise in an electric- 
circuit system under the action of an impressed electromotive 
force. In the method given in Chap. II, the subject is discussed 
from the standpoint of current distributions in a circuit system, 
arriving at the general-expression formula for the current rise 
in a circuit system under the application of an electromotive 
force, and from this the formula for current subsidence in a 
circuit system is deduced. 


1 Collected papers, Vol. II, pp. 37-373. 


APPENDIX 
NOTE ON BESSEL FUNCTIONS 
The Bessel functions of the first kind J,(z) and I,(z) are defined, 


when n is zero or a positive integer, by the absolutely convergent 
series: 


4 
J,(z) = a te Wie + 2) T 9-4@n + 2)(Q2n+4) 


78 
Q-4- 6Gn + 2)n+4)Qn £6) * oe | 
Be Hd x 
oS ane + 3@n $2) 24m + 2Qn +H 7 


7? 


2-4-6(2n + 2)(2n 4+ 4)(2n + 6) 


. 


The function J,,(x) satisfies the differential equation 


ad (x) 4 Lad, (2) Sis (: ae m2) J) aa 0, 


dx? aR ahs 


and the function J,(z) satisfies the differential equation 


ges 


dx? an as 


The connection between the J and J functions is given by the 
following equations: 


I(x) = 7-"J n(x), 
(j= V-}) 
Log(e) la): 
The series given above J,(x) and I,(x) become practically 


useless for numerical computations when the argument z is even 


moderately large. These functions, however, are expressible in 
158 
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other series forms which are well adapted for numerical calcula- 
tions for large values of the argument. Thus: 


ae An? —1 | (4n? — 1)(4n? — 9) 
le A/a 0 118 21(82)? 
(4n? — 1)(4n? — 9)(4n? — 25) ; 
31(82)3 

| 

ane Ja IP, cos (x es cian ',) = Qusin(z = m1, 

where ; 

oat (4n? — 1)(4n?— 9) 
ee 21(82)? af 
(4n? — 1)(4n? — 9)(4n? — 25)(4n? — 49) _ 
41(8x)4 
Piste Ant I) An 9) ant = 25) 
Be) a Sr 31(8x)3 


For very large values of z, the functions approximate the values 
given by 


I(t) = Fo 


JaA(z) = ie cos (2 -*nt1,), 


The following formulas are useful giving relations between the 
functions and their differentials, 7; and J/, indicating the differen- 
tials of the functions: 


Ft = Jn — Ines 
Ji = ete + inet 
x 
2n 
Int =, mak + Jnr1 = 0 


Jo = —di 
Similarly, for the J functions, 


I= In te ate 


U —an piesa, 


ll 


Tae sole + In = 0 
Le, 
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Derivation of formula, 


p In(at) 
t Gly 


@ — py — (am + 1)? 


Tin1(at 
= (2m + 1)a? mt 


This formula is used in Chap. V in the solution of wave propaga- 
tion problems (see p. 104). 
Disregard the constant factor a for the moment. We have, 
by direct differentiation, 
In(t) _ tpIn(t) — mt™In(t) _ pIm(t) _ mI m(t) | 


{™ {2m t™ tmtl 


By the relation given above, we have 


pIm(t) = FIm(t) + Imsilt). 


Hence, 


I PA: Lie 2 | (: 
a = —— ) a Lat ad 


Differentiating again, 


pS felt) p AG Eyer eat 


t” { = 
ne ae, Tal) 
“a t Pm ae 
pal Bale Tate) 
t Pr(g)m 


“e ie = 


fee a) nO Gerla 


2m Be ljm—i 
P (t)™ + r i Denth) ae In(d) aa 


m Lote 
a 


a 


_ =—2m I(t) 
See: 
a EL 
meee ae ae" 


Imlay Visa) 
ie t ike + 


2mIn(t) 1 Im(t) _ 2m In In(t) 


P tm iP t” t? tm ah {™ 


a — 


+ 


and 


sIn(t) _ = (2m + 1) Init) 
t™ 


Pa pes, 
Vee i tm 


(2m + vf gue 


ll 
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If we replace ¢ by at, we have 
me Ld 
P* d(at) ~ a dé 
ou L @, 
a’ dt? 


and the above transforms to 


a (at)™ 


or 


(pt — 22 = 2m + yale) 
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p TrAat) Talat); 
(1) age — Om + viens 


(a ). 
(a peas 


MATHEMATICAL FORMULAS 


For convenience, the more frequently 
formulas in the text are assembled here. 


e* = cosh x + sinh x 
e* = cosh x — sinh x 
@* = cosx+jsinz 
e7* = cosx —jsinz 
cos x = cosh jz; cosh x = cos jx 


used mathematical 


sin z = —j sinh jz; sinh x = —j7 sin jz 
tan ¢ = —j tanh x; tanh x = —j tan jr 


cos? z + sin? z = 1 


cosh? « — sinh? x 


sin 5 =4 5a — cos x); cos 


sinh = = [a(co (cosh x ey: cosh 5 


sin (x + y) ap ee 


cos (x + y) = cosxcosy + sinz sin y 


bN!s 


= ac + Gee 


- /3 (cosh x + 1). 


sinh (2 + y) = sinh xz cosh y + sinh y cosh x 


cosh (x + y) = cosh x cosh y + sinhzs 
sinh (mjr) = O (m is an integer) 
cosh (mjr) = (—1)” 
tanh (mjr) = 0 
sinh (x + mjr) = (—1)” sinh x 
cosh (x + mjr) = (—1)” cosh x 
sinh 2x = 2 sinh x cosh x 


inh y 


sinh 3z = 4 sinh? x + 3 sinh z = sinh z (4 cosh? x — 1) 
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sinh (n + 1)z = 2 cosh x sinh nx — sinh (n — 1)ax 


sinh nz = n cosh”! 2 sinh 2 + 


cosh 27 = 
Coshro7a— 
cosh (n + 


cosh nz = 


sin x 


COS & 


n(n — ae — 2) tduaoee 


sinh'za+:-:-: 


cosh? x + sinh? x = 2 cosh? z —1=1+4 2 sinh? x 
4 cosh® x — 3 cosh x = cosh z(4 sinh? x + 1) 
1)x = 2 cosh zx cosh nz — cosh (n — 1)z 


cosh” 2 -+- ——— es 9 rad) cosh” z sinh? x + - 


2 23 at 


x x 
ppeess haze esi iigane ans 


a3 7? al 


pairs (Ran ieee VA a aie 


Lig a4 76 


BE eee eee 


x a? at 


Seta ee 
eg at Cee 
mir otic ole 
2 2 
Shae 2 + ee 
OMTIOMTC 
r ze x 
=1+ : 4s om LE ee 


sin x 


cos 7 = 


7 
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(—1)""1(2n — 1)r 
sec 7 = > ae i = 
( 2 )e-= 


Pst! (—1)*-122 
cose x = | + > Boas 


nx? ay x 


1 
= 2 the Ut 
(Lt o)™ = 1 4 ma + MO De, mim — Mm — 2)? | 


Co a eas el i ec 
a A ee 


Ie dl JL oil oR Pod oiors 
+ 4 = As ee) = 54. 
(1 z) Se = f:4 § 9.4.67 2-4-6-8" 
ioe: aye 5 1983 5) 7/ 
oe Se c= ey 2 3 5 
os 2) : BF th. Ae ede a Gore Gas 
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TaBLE I.—BrsseL FUNCTIONS 


x | J (x) | J\(x) | I(x) Ii(x) 
0.0 1.00000 0.00000 1.00000 0.00000 
0.2 0.99002 0.09950 1.01002 0.10050 
0.4 0.96040 0.19603 1.04040 0.20403 
0.6 0.91200 0.28670 1.09205 0.31370 
0.8 0.84629 0.36884 1.16651 0.43286 
1.0 0.76520 0.44005 1.26607 0.56516 
1.2 0.67118 0.49829 1.39373 0.71468 
1.4 0.56685 0.54195 1.55340 0.88609 
L6G 0.45540 0.56990 1.74998 1.08481 
ies 0.33999 0.58152 1.98956 1.31717 
250 0.22389 0.57672 2.27959 1.59064 
2.2 0.11036 0.55596 2.62914 1.91409 
2.4 0.00251 0.52019 3.04926 2.29812 
2.6 —0.09680 0.47082 3.55327 2.195388 
2.8 —0.18504 0.40970 4.15730 3.30106 
3.0 —0.26005 0.33906 4.88079 3.95337 
3.2 —0.32019 0.26134 5.74721 4.73425 
3.4 —0.36430 0.17923 6.78481 5.67010 
3.6 —0.39177 0.09547 8.02768 6.79271 
3.8 —0.40256 0.01282 9.51689 8.14042 
4.0 —0.39715 —0.06604 11.30192 9.75947 
4.2 —0.37656 —0.13865 13 .44246 11.70562 
4.4 —0.34226 —0.20278 16.01043 14.04622 
4.6 —0.29614 —0.25655 19 .09262 16. 86256 
4.8 —0.24043 —0.29850 22 .79368 20 .25283 
5.0 —0.17760 —0.32758 27 .23987 24 .33564 
i) —0.11029 —(0.34322 32 .58359 29 .25431 
5.4 —0.04121 —0.34534 39 .00879 35. 18206 
5.6 +0 .02697 —0.33433 46 .73755 42 .32829 
5.8 +0.09170 —0.31103 56 03810 50 .94618 
6.0 +0.15064 —0.27668 67 .23441 61.34194 
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Ks Ce cae Cosh x Sinh x 
0.0 1.0000 1.0000 1.0000 0.0000 
OnzZ 1.2214 0.8187 1.02007 0.20134 
0.4 1.4918 0.6703 1.08107 0.41075 
0.6 1.8221 0.5488 1.18547 0.63665 
0.8 22225 0.4493 1.33743 0.88811 
1.0 2.7183 0.3679 1.54308 17520 
2 3.3201 0.3012 1.81066 1.50946 
1.4 4.0552 0.2466 2.15090 1.90430 
1.6 4.9530 0.2019 2.57746 DS lool 
1.8 6.0496 0.1653 3.10747 2.94217 
2.0 7.3891 Onsos 3.76220 3.62686 
Dee, 9.0250 0.1108 4.56791 4.45711 
2.4 O232 0.0907 5.55695 5.46623 
2:6 13.4637 0.07427 6.76900 6.69473 
2.8 16.4446 0.06081 8.25273 8.19192 
oO) 20.0855 0.04979 10.0677 10.0179 
me 2475325 0.04076 12.2866 12.2459 
3.4 29.9641 0.03337 14.9987 14.9654 
20) 36.5982 0.02732 18.3128 18.2855 
3.8 44.7012 0.02237 22.3618 22 .3394 
4.0 54.5982 0.01832 27 .3082 27 .2899 
4.2 66.6863 0.01500 33.3507 33}. SS A// 
4.4 81.4509 0.01228 40.7316 40.7193 
4.6 99.4843 0.01005 49.7472 49.7371 
4.8 PT 510 0.00823 60.7593 60.7511 
Sa, 148.413 0.00674 74.2099 74.2032 
Dee 181.272 0.00552 90.6388 90.6333 
71 221.406 0.00452 110.705 110.701 
5.6 270.426 0.00370 185215 15-211 
5.8 330.300 0.00303 165.151 165.148 
6.0 403 .429 0.00248 201.716 201.713 


166 


NON KE oO 
Oven: 


DPOo>»r WwW WwW 
aAononaod ao 


NIN 
ano 


80 


Jo (x) and J, (x) 


HEAVISIDE’S ELECTRICAL CIRCUIT THEORY 


2 


2 Vie 
Tasie III.—erfy = Vil, ela 
TJO 


erfy A y erfy 
0.05637 0.05637 1.05 0.86243 0 
0.11246 0.05609 1.10 0.88020 0 
0.16799 0.05553 1.15 0.89612 0 
0.22270 0.05471 1.20 0.91031 0 
0.27632 0.05362 1.25 0.92290 0 
0.32862 0.05230 1.30 0.93400 0 
0.37938 0.05076 1.35 0.94376 0 
0.42839 0.04901 1.40 0.95228 0 
0.47548 0.04709 1.45 0.95969 0 
0.52049 0.04501 1.50 0.96610 0 
0.56332 0.04283 1.55 0.97162 0 
0.60385 0.04053 1.60 0.97634 0 
0.64202 0.03817 1.65 0.98037 0 
0.67780 0.03578 1.70 0.98379 0 
0.71115 0.03335 1.75 0.98667 0 
0.74210 0.03095 1.80 0.98909 0 
0.77066 0.02856 1.85 0.99111 0 
0.79690 0.02624 1.90 0.99279 0 
0.82089 0.02399 1.95 0.99417 0 
0.84270 0.02181 2.00 0.99532 0 


INDEX 


A 


Activity, mutual, 145 
Alternating voltage, 6 
applied to filter circuits, 43 
Applied voltage at some phase 
angle, 20 
varying as e?', 97 
Artificial lines, 123 
distortionless, 133 
non-dissipating, 128 
series elements inductance and 
resistance, shunt capacity, 
and leakage, 140 
ocean cables, 124 


B 


Band tuning, 43 

Bessel functions, 158 
approximate formule, 159 
table, 164 

Bessel-Fourier series, 132 


C 


Cable, 65 
artificial, 124 
condenser at receiving end, 90 
transmitting end, 91 
finite length, 69 
with terminal impedances, 80, 
85 
grounded at both ends, 69 
infinite length, 66 
open at one end, 73 
resistance at receiving end; trans- 
mitting end grounded, 88 
voltage applied at intermediate 
point, 76 
Cables, ocean, 65 
artificial ocean, 124 


Calculus, operational, 1 
Capacity circuit, applied alternating 
voltage, 12 
and resistance circuits, 5 
Circuit network, 138, 15 


Circuits distributed inductance, 
capacity, resistance and leak- 
age, 94 

divided, 28 


electric filter, 38 
magnetically coupled, 31, 33 
four-coupled, 60 
three-coupled, 59 
of distributed capacity and resist- 
ance, 65 
Cohen, Louis, 15 
Conjugate property, 149 
Coupled, circuits dampings and 
frequencies, 58, 59 
negligible resistance, 34 
inductive circuits, 31 
oscillatory circuits, 33 
Critical case, oscillatory circuit, 25 
Current distribution in divided 
circuits, 30 


D 


Derivation of expansion formula, 15 
Dissipativity, 145 

mutual, 148 
Distortionless artificial line, 133 
Divided circuits, 28 
Dynamicial system, 149 


1D) 


Electric filter circuits, 38 
Energy, magnetic, mutual, 148 
mutual, 145 
subsidence, 145 


167 


168 


Error function, 68 
table, 166 
Expansion theorem, 13 
Exponential functions table, 165 
Extension of expansion theorem for 
alternation voltage, 18 


r 


Filter circuits, 38 
applied alternating voltage, 43 
electric, 38 
series elements inductances and 
capacities in series, shunt 
elements inductances and 
capacities in parallel, 53 
capacities, shunt elements 
inductances and capacities 
in parallel, 51 
inductances, shunt elements 
inductances and capacities 
in parallel, 50 
magnetically coupled, 54 
low-pass, 43 
high-pass, 47 
Four-coupled circuits, 60 
Fourier-Bessel series, 132 
Fractional derivatives, 67 
Free cables, ends grounded or 
open, 69 
vibrating circuit system, 38 


G 


Generalized impedance, 13 


H 


Heaviside’s derivation of expansion 
formula, 1438 

High-pass filter, 47 

Hyperbolic functions, table, 165 


I 


Impedance at cable terminal, 83 
generalized, 13 

Indices, negative, 10 

Inductive circuit, 2, 6, 10 
coupling, 146 

Infinite line, 95 
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L 


Line, finite, 103 
terminal impedances, 117 

grounded at far end, 112 
infinite, 95 
loaded, 137 
periodically loaded, 137 

Lines, artificial (see Artificial lines). 
non-dissipating, 100 
transmission, 94 

Low-pass filter, 43 


M 


Machinari Goto, 18 
Magnetically coupled circuits, 31, 33 
filter circuits, 54 
Mathematical formule, 161 
formulation of expansion theorem, 
13, 17 
Multiperiodic circuit system, 38 
Mutual electric energy, 148 
dissipativity, 148 
magnetic energy, 148 


N 


Negative indices, 10 
Non-dissipating lines, 100 
Normal distributions, amplitudes 
of, 152 
of voltage and current in 
circuit network, 151 
systems, 144 


O 


Ocean cables, 65 
Operand shifting, 9 
Operating on exponential function, 7 
on unity function, 3, 4 
Operational calculus, 1 
method, inductive circuit solution, 
2 
Oscillatory circuit, 25 
circuits coupled, 33 


le 


Partial fractions, 15 
Periodically loaded line, 137 


INDEX 


Propagation equations of voltage and. 
current, 65 
Pupin, M. I., 137 


R 


Resistance at cable terminal, 81 
coupling, 147 
operator, 154 
Resonance effect for frequencies of 
free oscillations of filter circuit, 
45 
of multiperiodic circuit system, 60 


Ss) 


Shift principle, 9 
Shifting of operand, 9 
Subsidence of current in circuits, 17, 
35 
of energy, 145 
Symbolic solutions, 3 


4 


169 


a 


Taylor’s theorem, 116 

Terminal impedances, cable, 80 
transmission line, 117 

Three-coupled circuits, 59 

Transmission lines, 94 

Transmitting end current, cable, 67 
infinite line, 96 

Two-coupled circuits, 59 


U 
Unity function, 3 
V 


Voltage, alternating, 6 
applied, at some phase angle, 20 
to filter circuits, 43 
applied varying as «~24, 108 
wave derivation of, 105 


WwW 
Wagner, Karl Willy, 145 
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